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We review applications of the sine-Gordon model, the 0(3) non-linear sigma model, the 
U(l) Thirring model, and the 0(N) Gross-Neveu model to quasi one-dimensional quan- 
tum magnets, Mott insulators, and carbon nanotubes. We focus upon the determination 
of dynamical response functions for these problems. These quantities are computed by 
means of form factor expansions of quantum correlation functions in integrable quantum 
field theories. This approach is reviewed here in some detail. 



1 



2 F.H.L. Essler and R.M. Konik 



Table of Contents 



1 Introduction 5 

2 Correlation Functions in Integrable, Massive Quantum 
Field Theories 7 

2.1 Computing Correlation Functions with Form Factors: Exact 
Results at Zero Temperature 7 

2.2 Spectrum and Scattering in an Integrable Model 9 

2.3 Computation of Form Factors: Form Factor Axioms .... 11 

2.3.1 Scattering Axiom 11 

2.3.2 Periodicity Axiom 12 

2.3.3 Annihilation Pole Axiom 12 

2.3.4 Lorentz Invariance 13 

2.3.5 Form Factor Normalization 13 

2.3.6 Minimality Principle 14 

2.4 Simple Example: Form Factors for the Sine-Gordon Model in 

the Repulsive Regime 15 

2.4.1 Soliton-Antisoliton Form Factor for the Current Op- 
erator 15 

2.4.2 Soliton-Antisoliton Form Factor for Non-Local Oper- 
ators 17 

2.5 Relevance of Higher Particle Form Factors: Quantum Ising 
Model as an Example 19 

2.6 Form Factors for Integrable Models with Bound States ... 21 

2.7 Form Factors in Finite Temperature Correlators 23 

2.7.1 Regularization of Form Factors 24 

3 Sine-Gordon Model and Spin-| Quantum Magnets 26 

3.1 The Spin-^ Heisenberg Chain in a Uniform Magnetic Field 26 

3.1.1 Excitation Spectrum of the Lattice Model 27 

3.1.2 Continuum Limit 28 

3.1.3 Operators and their Normalizations 35 

3.1.4 Dynamical Spin Correlation Functions 37 

3.2 Massive Perturbations and Window of Applicability of 
MIQFT 39 

3.3 Field-Induced Gap Problem 39 

3.3.1 Spectrum of the SGM and Quantum Numbers ... 40 

3.3.2 Thermodynamics 42 



Integrable QFTs in condensed matter physics 3 

3.3.3 Normalization of Form Factors in the Sine-Gordon 
Model 44 

3.3.4 Dynamical Structure Factor 45 

3.4 Transverse-Field Model 52 

3.5 Dimerized Chain 55 

3.6 Quasi-ID Spin-^ Antiferromagnets in the Ordered Phase . . 58 

3.6.1 Mean-Field Theory 59 

3.6.2 Random Phase Approximation for the Interchain Cou- 
pling 62 

4 0(3) Non- linear Sigma Model and Gapped, Integer Spin 
Chains 64 

4.1 From Integer Spin Heisenberg Chains to the 0(3) Non- linear 
Sigma Model 66 

4.2 Basic Description of the 0(3) Non- linear Sigma Model ... 68 

4.2.1 Form Factors for the Staggered Magnetization Opera- 
tor, n 69 

4.2.2 Form Factors for the Magnetization Operator, M . . 70 

4.3 Zero Temperature Correlators 71 

4.3.1 k ^ Magnetization Correlations 71 

4.3.2 A; ~ vr Magnetization Correlations 72 

4.4 Low Temperature Properties of Correlation Functions ... 75 

4.4.1 Magnetization 76 

4.4.2 NMR Relaxation Rate 80 

4.4.3 Drude Weight of Spin Conductance at Finite Field . 83 

4.5 Transport: Ballistic or Diffusive 85 

5 U(l) Thirring Model and Quasi One Dimensional Mott 
Insulators 90 

5.1 Lattice Models of Correlated Electrons 91 

5.2 Field Theory Description of the Low Energy Limit 92 

5.2.1 Half Filled Band 92 

5.2.2 Quarter Filled Band 96 

5.3 Correlation Functions 97 

5.3.1 Correlation Functions in the Charge Sector 98 

5.4 Optical Conductivity 98 

5.5 Spectral Function 101 

5.5.1 Half-Filled Mott Insulator 101 

5.5.2 Zero Temperature 102 

5.5.3 Finite Temperature T < A 104 



4 F.H.L. Essler and R.M. Konik 



5.5.4 Quarter-Filled CDW Insulator 104 

5.6 A Remark on Luttinger's Theorem 107 

5.7 Interchain Tunneling 108 

5.7.1 Expansion Around Uncoupled Chains 108 

5.7.2 Formation of a Fermi surface 109 

6 Hubbard Ladders and Carbon Nanotubes: SO (8) Gross— 

Neveu model 111 

6.1 Armchair Carbon Nanotubes and Hubbard Ladders: Identical 
Low Energy Behavior 112 

6.2 Weak Coupling Flow Onto S0(8) Gross-Neveu 116 

6.3 Limitations of the RG Analysis 120 

6.4 Excitations and Physical Fields in S0(8) Gross-Neveu . . . 124 

6.4.1 Excitation Spectrum 124 

6.4.2 Relationship between Lattice Operators and Gross- 
Neveu Fields 125 

6.5 Form Factors in S0(8) Gross-Neveu 127 

6.5.1 One and Two Particle Current Form-Factors .... 127 

6.5.2 One and Two Particle Kink Form-Factors 128 

6.6 Exact Low Energy Correlation Functions in S0(8) Gross- 
Neveu 129 

6.6.1 Behavior of Optical Conductivity in a Hubbard Ladder 129 

6.6.2 Single Particle Spectral Function 133 

6.6.3 STM Tunneling Current 135 

6.7 Effect of Integrable Breaking Perturbations 138 

References 142 



Integrable QFTs in condensed matter physics 5 

In memory of Ian Kogan. 

1. Introduction 

The study of strongly correlated electrons in low dimensional systems lies 
at the heart of much of modern condensed matter physics. Interest in these 
systems arises as the behavior of strongly correlated systems is not in gen- 
eral adequately captured in approximations based upon 'free-particle' non- 
interacting models. In the presence of generically strong interactions, the 
physics is typically much richer, exhibiting qualitatively new features. 

Accessing this physics presents an imposing challenge. In one route to 
comprehension, quantum field theories can be employed. Field theories are 
typically able to describe the low energy behavior of strongly correlated 
systems and so are able to extract universal characteristics. It is these char- 
acteristics that are of greatest interest precisely because they provide the 
most robust experimental signatures and do not, in general, depend upon 
particular experimental details. 

In low dimensions, the study of numerous quantum field theories is aided 
by their integrability. Integrable theories are characterized by an infinite 
number of non-trivial conserved charges. The existence of these charges al- 
lows for an exact characterization of many features of these models. Both 
the spectrum and the scattering matrices of an integrable quantum field 
theory can be explicitly written down. On the other hand, correlation func- 
tions in an integrable quantum field theory cannot, in general, be exactly 
computed. However there exist approaches, based on exploiting the inte- 
grability of the model, that can be used to obtain some information on 
correlation functions. And while incomplete, this information does reflect 
the non-perturbative structure of the theory. 

Of these approaches, we describe in this review one based upon form- 
factors. Form factors are matrix elements of quantum fields with exact 
eigenstates. Under a spectral decomposition, all correlation functions can 
be written in terms of form factors. While this decomposition is exact, its 
practical manipulation requires truncation of the spectral sum. In this review 
we will explore how this truncation may be done and under what conditions 
exact information on the theory remains available. As the reader will see, the 
behavior of correlation functions at low energy scales fortuitously remains 
exactly computable. 

The review is organized as follows. In the first substantive section. Section 

2, we give in detail the form factor programme to computing correlation 
functions in integrable quantum field theories. Our purpose here is to be 
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pedagogical. We begin by giving an overview of how form factors are used 
to calculate zero temperature correlators. In Section 2.2 we then turn to 

describing the key features of an integrable model for the programme, the 
spectrum and its associated scattering. Form factors are computable in an 
integrable model because they must satisfy a series of constraints arising 
from both consistency with the underlying scattering and various assumed 
analyticitics in the energy-momentum. We list these constraints in Section 
2.3. In Section 2.4 we give examples of the application of these constraints to 
an archetypal integrable theory, the sine-Gordon model. While form factors 
do not generically give the behavior of correlation functions at all energy 
scales, they do provide exact information at low energies. And at higher 
energy scales, corrective terms, as a rule, are extremely small. We illustrate 
this principle in Section 2.6 with the specific example of the Ising model. In 
the final two subsections we turn to form-factors in more involved settings: 
form-factors in models with bound states and the use of form-factors to 
compute correlation functions at finite temperature. 

Having given an overview of the form-factor programme, we turn to spe- 
cific applications. As our first example, we consider applications of the 
sine-Gordon model to half-integer spin chains. An anisotropic Heisenberg 
spin-1/2 chain (in a magnetic field) has gapless excitations. In this section 
we first detail the bosonization of the spin chain and so exhibit the chain's 
Luttinger liquid phase. We then consider a number of physical perturbations 
under which the chain becomes gapped and is described by the massive sine- 
Gordon model. These perturbations include a staggered magnetization, a 
transverse magnetic field and dimerization. The dynamical structure fac- 
tor, which is measured in inelastic ncTitron scattering experiments, can be 
determined using sine-Gordon form factors. 

In Section 4, we study integer spin chains using the same techniques. 
Unlike their half-integer counter parts, an integer spin Heisenberg chain is 
completely gapped and its low energy behavior is believed to be described by 
the 0(3) non-linear sigma model without topological term. For pedagogical 
reasons we begin the section by giving the map between the two, pointing 
out what approximations are used. We then turn to a basic description of 
the 0(3) non-linear sigma model including its spectrum, scattering matrices, 
and form factors of various physical fields. With these form factors in hand, 
we next consider the computation of various zero temperature spin-spin cor- 
relation functions that would be measured in neutron scattering experiments. 
Having considered the correlation functions at T = 0, we turn to their finite 
temperature counterparts. This will be our primary example demonstrating 
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the possibility of a form-factor computation at T 7^ and it is done in some 
detail. Using these techniques we will address the issue of whether transport 
is ballistic or diffusive at finite temperatures in integer spin chains. 

In the next section, Section 5, we turn to applications of the U(l) Thirring 
model to quasi-one-dimensional Mott insulating chains. Such insulators are 
modeled by various extended Hubbard models which in turn are related to 
the U(l) Thirring model. The latter can be bosonized in terms of a free 
boson and a sine-Gordon model. Using the form factors of various operators 
in the sine-Gordon model we compute both the optical response as well 
as the single-particle spectral function of half-filled and quarter filled Mott 
insulating chains. While this treatment is appropriate to one-dimensional 
materials, we also show how the case of weakly coupled chains can be treated. 

In the final section. Section 6, we consider our last application of form 
factors, that to Hubbard ladders and armchair carbon nanotubes. In the 
weak interaction limit, both materials can be described in an RG sense by 
the S0(8) Gross-Ncvcu model, an integrable theory of four interacting Dirac 
fermions. We outline this mapping in some detail and point out its limita- 
tions. We then use the form factors of SO (8) Gross-Neveu to compute a 
number of physical quantities including the optical response and single par- 
ticle spectral function. We end by pointing out how the results are changed 
if small, integrable-breaking perturbations are introduced. 

2. Correlation Functions in Integrable, Massive Quantum 
Field Theories 

In this section we outline the form-factor programme by which correlation 
functions in massive integrable field theories can be computed. 

2.1. Computing Correlation Functions with Form Factors: 
Exact Results at Zero Temperature 

In general computing correlation functions in integrable field theories is an 
open problem. There exists no general technique that is able to access a 
generic correlation function at all energy scales. Now it is certainly true 
that progress can be made in specific instances. For some integrable models, 
such as Ising model variants, the sine-Gordon model at is free fermion point, 
and the Bose gas with delta-function interactions, some correlations can be 
determined with the aid of Fredholm determinants, see e.g. [36, 58, 59, 90- 
92,124-126,152,162,174-177,183,191] and references therein. This method, 
however, is technically involved, and at present only works in a handful 
of cases. More promising is the form- factor programme for the calculation 
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of correlation functions [167, 280]. Form factors do not allow in practice 
the complete determination of a correlation function. However in a gapped 

(massive) field theory, they do permit the exact determination of the low 
energy properties of the corresponding spectral function. 

This arises as the form-factor representation of any zero tempera- 
ture correlation function is obtained by inserting a resolution of the 
identity corresponding to the basis of eigenstates. (We will con- 
sider finite temperature correlation functions later in this section.) 
Thus, for an operator, 0(x,t), we write the spectral decomposi- 
tion schematically (r denotes imaginary time, and T time ordering) ^ 

G?(x,r) = -(0|r(O(a;,r)Ot(0,0))|0) 

= - J] J]e--^-(0P(x,0)|n;s„)(n;5„pt(0,0)|0), (r > 0), 

n=0 Sn 

where Es„ is the energy of an eigenstate, |n; Sn), with n particles described 
by quantum numbers, {s„}. By inserting a resolution of the identity, we 
have reduced the problem to one of computing individual matrix elements. 
In an integrable model the matrix elements of a physical operator between 
the vacuum and the exact eigenstates can in principle be computed exactly 
from the two-body S-matrix. However the calculation of these matrix el- 
ements, as well as the evaluation of the sums/integrals. V . becomes in- 
creasingly cumbersome as the particle number n becomes large, so that the 
full expression for the correlation function cannot be evaluated in closed 
form. Often, however, a truncation of the sum at the level of two or three 
particle states already provides a good approximation to the full correlation 
function [54, 72, 73, 197, 201]. This may be understood in terms of phase 
space arguments [54,230]. On the other hand, this truncation is no longer 
necessary in a massive theory, if one considers the corresponding spectral 
function. Only eigenstates with a fixed energy, a;, contribute to the spectral 



^ The sum over s„ is meant to include integrals over the momenta of all particles and sums over 
particle/excitation types. 



Integrable QFTs in condensed matter physics 



9 



function: 

--ImG^(x, -iLV + 6) 

TT 

oo f 

= EE 1 (0P(x,0)|n;s„)(n;s„|Ot(0,0)|0)<5(a; -E,J 

n=0 ^ 

-e(0|Ot(0, 0)|n; s„)(n; 0)|0)(5(u; + E^J^ , (2.2) 

where e = zt for fields, O, that are bosonic/fermionic. 

In a massive theory the creation of an extra particle in the intermediate 
exact eigenstate costs a finite amount of energy, and so the sum in Eq. (2.2) 
is finite. For example, when uj is smaller than the energy of all three-particle 
states (i.e. when to is below the three-particle threshold), then only the form 
factors with one and two particles (n = 1,2) have to be determined in order 
to obtain an exact result. All of the spectral functions we have computed in 
this review are exact at sufficiently low energies for this reason. 

As knowledge of the exact eigenfunctions in an integrable model is central 
to computing correlations functions, we consider this in more detail in the 
next section. 



2.2. Spectrum and Scattering in an Integrable Model 

The key feature of an integrable system is the exact knowledge of a basis of 
eigenstates of the fully interacting Hamiltonian. At the root of integrability 
is a well defined notion of "particles", or "elementary excitations" in the 
fully interacting system. These particles scatter off each other according to 
two-body S'-matrices, that is, all particle production processes are absent 
and particle number is conserved. This is due to special conservation laws 
which exist in an integrable model, preventing the decay of these particles. In 
this sense, an integrable system is similar to a Fermi liquid. An additional 
feature is that new particles can arise as bound states of already existing 
ones. However the total number of different types of particles is finite which 
makes the system analytically tractable. 

Formally the elementary excitations arc created and destroyed through 
the Faddeev-Zamolodchikov operators, denoted by Aa{6). 9 is termed the 
rapidity and it encodes the energy-momentum carried by the excitation via 

P = Asinh((9); E = AcoshiO). (2.3) 

Two excitations will in general scatter according to a non-trivial two-body 
S-matrix, S. S gives the amplitude of the process by which two particles 
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{a, b} scatter into two potentially different particles, {a', b'}. In terms of the 
Faddeev-Zamolodchikov operators, this scattering determines the commu- 
tation relationship between operators 

Aai9i)Abie2) = Si^\9i - e2)A'(02)^a'(^l); 

A\{eMl{02) = si"' (9, - 02)4'(^2)4'(^i); 

^a(^l)4(^2) = 27T6abS{9i - 62) + SlUOi - ^2)4' (^l) A' (^1 ) ■ (2-4) 

It is solely a function of 9i — 92 = ^12 by Lorentz invariance. 



time 




Figure 1. A graphical representation of the Yang-Baxter equation. The left-hand figure repre- 
sents a three body process by which three particles a,b, and c scatter into a",b", and c". Via 
integrability, this three-body process can be factorized into two different sets of three two-body 
processes pictured in the central and leftmost figures. Under the Yang-Baxter relation, these 
different sets of two-body scattering processes are equivalent. 

These two-body S-matrices encode all scattering information in the the- 
ory. Scattering processes involving higher number of particles can always be 
expressed in terms of two-body scattering matrices. Suppose we consider a 
three particle scattering process, scattering particles {a,b,c} into {a',b',c'} 
given by the S-matrix 5"^^ . This three-body S-matrix factorizes into a set 
of two-body S-matrices 

ca'h'c' (a a n \ _ oa'b'/n a \Qa"c' (a a \cb"c" /n n \ 

= St'{92,9,)S'if{9u9,)Si!^"{9,,92). (2.5) 

From the above equation we see that we can factorize the three-body S- 
matrix in two different ways. This is illustrated graphically in Fig. 1. Because 
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the theory is integrable, the different ways of factorizing are equivalent. This 
equivalence is known as the Yang-Baxter equation. The ability to factorize 

the higher-body S-matrices results from the existence in the integrable theory 
of conserved charges with are non-trivial powers of energy and momentum 
[194,299,329-331]. 

As two-body scattering provides complete information on an integrable 
theory, the form factors (even those involving large number of particles) are 
determined by this S-matrix alone. We will see this in operation in the next 
section. 

Using the Faddeev-Zamolodchikov operators, a Fock space of states can 
be constructed as follows. The vacuum is defined by 

Aa{9)\0)=0. (2.6) 

Multiparticle states are then obtained by acting with strings of creation 
operators aI{9) on the vacuum 

|^n...ei)a„...a, =4„(^n)...4i(^l)|0) • (2.7) 

In terms of this basis the resolution of the identity is given by 

1 = |0)(0| + EE/ ^^7^^ (2-8) 

2.3. Computation of Form Factors: Form Factor Axioms 

The form factors of a field O are defined as the matrix elements of the field 
with some number of particles, Aa{0): 

C-aJdl, ■■■,dn) = {O{0, 0)AaM ■ ■ ■ Aa^m. (2.9) 

These matrix elements are constrained by a variety of requirements arising 
from the scattering relations just discussed, Lorentz invariance, hermiticity, 
analyticity, and the locality of the fields. We consider each in turn. 

2.3.1. Scattering Axiom 

For the form factor to be consistent with two body scattering we must have 

Sai(H+-^^{6i-0i+i)f2^^...^a',,a',+ i,---,aS^^^''' i ^« ' 1 ' ' ' '^n)- (2.10) 

This relation is arrived at by commuting the i-th and i + 1-th particle and 
using the Faddeev-Zamolodchikov algebra in Eq. (2.4). 
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2.3.2. Periodicity Axiom 

A second constraint upon the form factor can be thought of as a periodicity 
axiom. In continuing the rapidity, 9, of a particle to 9 — 2Tri, the parti- 
cle's energy-momentum is unchanged. However the form-factor is not so 
invariant. We instead have 

/a^,...,a„(^l>--- ,9n) = fZ,a,,-,ar.-A^n-2TTi,9,,--- ,9n-l). (2.11) 

This constraint is derived from crossing symmetry [36] . It implicitly assumes 
that the field O is local: if O is non-local additional braiding phases appear 
in the above relation [37, 185,279]. 

These braiding phases arise when two fields are interchanged: 

^{x, t)0{y, t) = R^oO{y, t)i;ix, t); x<y. (2.12) 

(On occasion, wc must deal with braiding matrices, not merely phases. But 
we will not consider such a situation in this work.) If R^o is non-trivial we 
must alter the above periodicity axiom to read 

fau-,aJOl,--- ,0n) = RrPr.ofZ,a,,-,a„.A0n-2Tri,9,,--- ,^n-l). (2-13) 

Here ipn can be thought of as the field which creates the excitation An- 

In order to employ the periodicity axiom (Eq. (2.13)), we then need to 
know how to specify the braiding of the fields. In one approach, we identify 
both fields with their corresponding excitations, Aq and An- If both fields 
are right-moving, the braiding of the fields is then encoded in the asymptotic 
limits of the corresponding S-matrix (see [279]): 

Ron = {SSni+<^))- (2.14) 

If on the other hand the fields are left-moving, we find instead Ron = 

2.3.3. Annihilation Pole Axiom 

Another condition related to analyticity that a form factor must satisfy is 
the annihilation pole axiom. This condition arises in form factors involving 
a particle and its anti-particle. Under the appropriate analytical continua- 
tion, such a combination of particles are able to annihilate one another. As 
such this condition relates form factors with n particles to those with n — 2 
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particles, 

{&n-ln-2)^, (2.15) 

where C is the charge conjugation matrix. This relation as written assumes 
that we are normalizing our particle states as {6\9') = 27r6{6 — 9'). Note 
that the braiding phase plays a role in this axiom (see r.h.s. of the above 
equation). 

2.3.4. Lorentz Invariance 

The form factor miist also satisfy constraints coming from Lorentz covari- 
ancc. In general, the form factor of a field, O, carrying Lorentz spin, s, must 
transform under a Lorentz boost, 6i ^ 6i + a, via 

+ «,••• ,0„ + a) = e-/„%jei,--- ,en). (2.16) 

Often in this work we will consider the correlation functions of (topological) 
current operators. Let us thus consider jo{x, t), a charge density, and ji (x, t), 
its corresponding conserved current. Together they form a Lorentz two- 
current. (Here 0, 1 are Lorentz indices.) The form factors for jo and ji 
appear as 

/ar--a„(6'l, • • • ,Gn) = efxuP''{9i)fai--an{9l, - ■ ■ ,9n), (2-17) 

where = ^•Acosh(^?j) and = ^^Asinh(^j). The function, 
/ai---an(^i) ■ ■ ■ ) 0n), on the r.h.s. of Eq. (2.17) is solely a function of 6i — 6j. 

2.3.5. Form Factor Normalization 

It is sometimes possible to determine the absolute normalization of a set of 
form-factors. (Of course, normalization of form-factors with different particle 
numbers but of the same field can be fixed, among other ways, by the use 
of the annihilation pole axiom.) In the case of current operator we can rely 
upon the action of the conserved charge 



Q = dxjo{x,0), 
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upon the single particle states with charge q. We expect (with the normal- 
ization of Section 2.3.3) 

{e,q\Q\e',q)=27rq6i9-9'). (2.18) 

Using crossing, we can relate this matrix element to a two-particle form 
factor: 

{e,q\Q\e\q) = {Q\e',q;9-iTr,q), (2.19) 

where q denotes the charge conjugate of q. Thus we are able to fix the 
normalization of the two particle form factor in a natural fashion (and so all 
other higher particle form factors through the annihilation pole axiom). 

With operators other than currents, we can still fix the phase of the 
normalization using hermiticity. For this purpose it is sufficient to consider 
2-particle form factors. Hermiticity then gives us 

= (Ot(o, 0)^„-, (01 - iTr)Aa-,{92 - in)) , (2.20) 
where the last line follows from crossing and so 

fa^a,iOi,02)* = fgU92 - in,9i - in). (2.21) 



2.3.6. Minimality Principle 

These conditions do not uniquely specify the form factors. It is easily seen 
that if f{9i, • • • , 9n)ai,-- ,an satisfies these axioms then so does 

f(ft ft ) Pn(cosh(%)) 

where P„ and Qn are symmetric polynomials in cosh(0jj), 1 < i,j < n, and 
are such that 

Pn\en=0n-i+m = Pn-2 Qn|6»„=6»„_i+7ri = Qn-2 ■ (2.23) 

To deal with this ambiguity, we employ a minimalist axiom. We choose Pn 
and Qn such that Pn/Qn has the minimal number of poles and zeros in the 
physical strip, Re(0) = 0, < Im^? < 2tt. Additional poles are only added in 
accordance with the theory's bound state structure. Using this minimalist 
ansatz, one can determine Pn/Qn up to a constant. 
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2.4. Simple Example: Form Factors for the Sine- Gordon 
Model in the Repulsive Regime 

In this section we consider form factors of fields in the sine-Gordon model 
in its repulsive regime. The sine-Gordon model is described by the following 
action 

S = — [ dxdr [da^d^"^ - /xcos(/3$)] . (2.24) 
iDvr J 

For (3 > l/\/2 this theory describes repulsively interacting solitons alone, i.e. 
there are no bound states. Classically these solitons arise as interpolations 
of the field $(x, t) between minima of the cosine potential. The two solitons 
are characterized by a topological U{1) charge, +/-. 

As it is integrable, the model is characterized solely by a two particle 
S-matrix. Its nonzero, U (l)-conserving elements are given by 



Slli9) = SZZ{9) = Soi9) = -exp 



dx , /■x6\ sinh((^ — |)a;) 

— sinh — — 7- 

L^o X Vztt/ cosh(|)sinh(^) 



sinh(f) 

+-V V -+y J sinh(^) ^ ^ 

zsin(f) 

St^iO) = S^l{9) = - . ;X Soi9), (2.25) 

smh(2-^) 

where 

^ = T^- (2.26) 

With this brief description in hand we go on to compute two particle form- 
factors in this theory. 



2.4.1. Soliton-Antisoliton Form Factor for the Current Operator 

We will first consider the form factor of the current correlator with a two 
soliton state. The (topological) current operator in sine-Gordon is given by 

f = e^^d^^ , (2.27) 

(where e*^^ = 1). As the operator itself carries no U{1) charge, it couples to 
a soliton-anti-soliton pair. The matrix element has the general form 

fe',e,i0uO2) = (01/(0,0)14^(^2)4.(^1)) • (2-28) 
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To determine /Ilea (^i; ^2), we note that the current couples anti- 
symmetrically to the soliton-anti-soliton pair (as is clear if one examines 
the limit = 1/2 where the model reduces to free massive fermions). Thus 
/eie2(^i)^2) takes the form, 

/fi«(^i,^2) = e.i«r(^i,^2), (2.29) 

where e is the anti-symmetric tensor. We can scalarize the above by explicitly 
exhibiting the piece of the form factor satisfying Lorentz covariance 

^{01,02) = (^e(^i+^^)/2 - (-l)/^e-(^^+''^)/2) f{eu) . (2.30) 

Having so constrained the form of /ile2(^i) ^2), we now apply the scattering 
axiom. Using the anti-symmetry of /Ilea (^1)^2) in the indices €1,62, we find 

/e'^.ei (^2,^1) =5(^12)/e'^,ei (^1,^2), (2-31) 

where 

sinh(l) - isin(f ) / f°° dr /tO\ \ 

^(^) = ie-.^, So{e) = - exp / ^ sinh (^) Ge(x) , (2.32) 

smh(2-5^) \Jo X ^*7r/ / 



Gc{x) = Gciix) + Gc2{xy, 

Q . . _ sinh((| - l)x) cosh((| - l)x) _ ^ sinh((^- |)x) 

' sinh(f) cosh(^) ' ^ cosh(f)sinh(f ) ■ 

We note that S{6) tends to —1 as ^ goes to zero. This implies the form factor 
will vanish in the same low-energy limit. As the current operator is local 
and bosonic, the braiding here is trivial. The periodicity axiom reading, 

/e'^,e,(^2,^l) = /e1e.(^l-27r^0l) (2.34) 

thus reduces to 

fi-e) = fi9-2ni). (2.35) 
A minimal solution satisfying these constraints on the form factor is 

ft,,,iei, 92) = aee„e,(e(^^+^^)/2 - (-l)'^e-(^i+^^)/2)s(0i2/2) 



X exp 



r--?#^Sin^(;^(m + 012)) 
.^0 



X sinh(x) 27r 



, (2.36) 
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where lA is some normalization with mass dimension [m]. The phase of A 
is determined through the hermiticity condition 



(2.37) 



This impHes that A is real. 
We note that 



exp 



cosh(^: 



(2.38) 



The form factor can then be written in the same form appearing in Ref . [280] : 



(6^,92) = iAe,, ,,(e('^^+^^)/2 - (-i)Me-('^i+^2)/2^ K^i2/2) 

Je,e^\ 1) 61,62 V V ) ^cOsh(^i|tHL) 



X exp 



-JO 



X smh(x) 27r 



, (2.39) 



2.4.2. Soliton-Antisoliton Form Factor for Non-Local Operators 

We now consider the two soliton form factor of the non-local operator, e^*2 
This operator's non-locality can be seen if we relate the field $ to the current 
operator 

g±if * ^ g±i| Sl^ dx'f{x,t)_ ^2.40) 

This non-locality will lead to a non-trivial braiding relation for the field. 

Again these operators carry no i7(l) charge and so will couple to a soliton- 
anti-soliton pair. It will prove to be convenient to consider form-factors 
involving symmetric and anti-symmetric combinations of the solitons: 

/|(ei,^2) = (O|e±^t*|^L(02)4(^i)) + (0|e±^f*|4(^2)^L(ei)); 

f^{ei,02) = (0|e±'t*|Al(e2)4(^i)) - (0|e±^i*|^i(^2)^L(ei)) • (2.41) 

As the operators e^*2* are Lorentz scalars, these form factors are solely a 
function of 612 = 0\ — 92- Applying the scattering axiom to the form factors 
we find 

ff{-e) = Ss{e)ff{e)- 

fi{-e) = SA{e)f^{e)- (2.42) 
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where 



sinh(^) / foc dx xO \ 



where Gc2{x) is as for the current form- factor. The non-locahty of the fields, 
gitij*^ relative to the solitons imphcs that R should be taken to be —1 in 
the periodicity axiom (see Eq. (2.13)). Thus we have 

/^_(-f?) = -/^+(e-27ri); 

f^4-e) = -f^_{e-2T:i). (2.44) 
In terms of the symmetric and anti-symmetric combinations we obtain 

fi{-e) = fi{e-2T,i). (2.45) 

Having expressed the constraints on the form factors in this way, we can 
readily write down a minimal solution 

. sinh(^i2/2) r dx Gc2(x) . o/^ /■ n \< 

f^{0i,02)=At \ . exp / T^^^sm-^i—ii-K + 612)) 

^ Ssinh(toi) ^L7o X sinh(x) ^tt^ 

^ (2.46) 

^ , sinli(0i2/2) r dx Gc2(x) . 9, x ,. ^ ,,1 

f^{eue2) = A^ v.X+^m "^p / -:^7e77IXSln2(-(^7^ + M) 



X sinhs(x) 27r 



where A^^^ are normalization constants. We now turn to their determina- 
tion. 

Prom the hermiticity condition, 

/|/^(^i, ^2)* = /|/^(^2 - ^7^, 91 - iTT), (2.47) 

we see that Ag* = A^ and A^* = —A"^. The operators e^*2* have a nonzero 
vacuum expectation value. Using the annihilation pole axiom it is possible 
to relate the value of A"^ to this expectation value. In this case it reads 

iTes0,=e,-nifi (01,02) = 4(e±4*). (2.48) 

We thus obtain = |(e^*^*). Using the transformation properties under 
charge conjugation we conclude that Ag = Ag. To relate A^ and A^ we 
appeal to the explicit calculation in Ref. [280] which shows that Ag = . 



Integrable QFTs in condensed matter physics 



19 



2.5. Relevance of Higher Particle Form Factors: Quantum 
Ising Model as an Example 

We have argued in Section 2.2 that to obtain exact information on spectral 
functions at low energies it is enough to include form factors involving only 
a few particles or excitations. However as a rule of thumb, form factor cal- 
culations fare even better. In practice, form factor sums have been found to 
be strongly convergent for operators in massive theories [54,72,73]. To ob- 
tain a good approximation to correlators involving such fields at all energies^ 
only the first few terms need to be kept. Even in massless theories where 
there are no explicit thresholds, convergence is good provided the engineer- 
ing dimension of the operator matches its anomalous dimension [197,201]. 
At high energies where higher particle form factors do begin to contribute, 
their contributions are progressively (often exponentially) smaller. We will 
illustrate this with a simple example, the spectral function of the spin-spin 
correlator in the Ising model. 

The quantum Ising model is a model of one dimensional spins governed 
by the Hamiltonian 



H = -J^(afaf+i+5<), 

i 

where and are Pauli matrices and J is assumed to be positive. At 
zero temperature, this model undergoes a Z2 phase transition as a function 
of g between an ordered and a paramagnetic state. This transition is in the 
same universality class as a classical two dimensional Ising model. In the 
continuum limit, it is thus described by a massive Majorana (real) fermion: 

S = ^ J dxdr (^^dg^ + tpd^ip + 2imip^^ , (2.49) 

where z/z = {T±ix)/2. ip and are the left and right moving self-conjugate 

components of the Majorana spinor. The mass of the fermion corresponds 
to the distance from criticality m ~ {gc — g) which we will assume to be 
positive (the model is in its ordered phase). The fundamental excitation, 
Aa{9), of the model arises from the mode expansions of the Fermi fields: 

^ = ^ r _^e-^/2fA(e)e-'"(^«"'+^"«') - ^t(e)e'"(^e-«+5e«)\ 

= -i^ r -^e^/2fA(^)e-'"(^^"'+^"^') + ^t(^)e-(^e-«+.V)Y (2.50) 
J -00 \/2m \ J 

The defining commutation relations of A and are {A{9),A\d')} = 
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2TT8{e - e'). 

While the form factors of the Fermi fields are trivial (only the one particle 

form factors arc non-zero), the form factors of the order parameter field, 
az{x,t), arc more complicated (ultimately a result of the spin and Fermi 
fields being mutually non-local). From Refs. [40,53,183,321], we have 



(0|a^(0, 0)|^(^2n) • • • ^(^l)) = fi0l, ■■■ , d2n) 

2n 

e 

i<j 



ntanh(^). (2.51) 



Only even numbers of fermions couple to the spin field (in contrast, odd 
numbers of fermions couple to the dual disorder field). The fcrmion S- 
matrix of this model is S = —1. The mutual non- locality of the spin, az, 
and Fermi, i/;, fields is encoded in the braiding phase, Ra^ip = — 1- Knowing 
S = Ra^ip = —1, it is easy to check that the form factors in Eq. (2.51) satisfy 
the axioms of Section 2.3. 

The spectral function of the spin-spin correlator is defined as 



S{u;, k) = — Im 



1 poo poo 

- dx dTe''^^-'''%-{T{azix,T)aziO,0))))l^-^u.^ 

J —oo J —oo 

(2.52) 

Here {T{az{x,T)az{0,0))) is the time ordered correlator in imaginary time. 
If we insert a resolution of the identity in between the two fields, do the 
Fourier transforms, and take the appropriate analytic continuation, we find 
at positive frequencies: 



2n „(VD 2n 



S(a; > 0, A;) = 27r 2^ — JJ / d9i (5(A; - m ^ sinh(0i)) 

n=i i=i ''-'^ i=l 

2n 

xS{u-mJ2cosh{ei)) 1/(^1,- •• ,92n\\ (2.53) 



i=l 



where dO = d6/{2iT). Integrating over 6i and 62 we obtain 

2 00 -. 2n 

sico>o,k) = -Y — l[ de. 

n=l ■ i=3 " 



' —00 



X 



-Vk^+^)\f{9ir-- ,e2nf, (2.54) 
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where 6i and 62 are given by 




i=3 
2n 



k' 



k — my ^ smh(gj); 

i=3 




(2.55) 



We plot as a function of energy, (j, the first three non-zero contributions 
(the 2-particle, the 4-particle, and the 6-particle) to the spectral function 
in Fig. 2. We first note that the 2n-particle form-factor only contributes for 
energies u > 2nm. Thus the plot of these contributions is exact for uj < 16m. 
We also see that for energies beyond 4m the higher particle contributions 
barely register. The 4-particle contribution, at its peak, is still 1/3000 of 
the 2-particle contribution. And the 6-particle contribution is considerably 
smaller yet. To be seen, it must be plotted on a scaled magnified by 10^. 
We thus learn that to compute the spectral function to energies far above 
the mass gap, it would seem to be sufficient to consider only the two-particle 
contribution. While we have only shown this to be the case in the off-critical 
Ising model, it is a generic feature of all massive integrable models. We will 
see other examples in Sections 4 and 5 with the spin spectral function of 
the 0(3) non-linear sigma model and the optical conductivity in the U(l) 
Thirring model. 

2.6. Form Factors for Integrable Models with Bound States 

When there are bound states in an integrable model, additional analytical 
requirements are imposed upon the form factors. Suppose we have an exci- 
tation, Ac, which is a bound state of excitations, Aa and A^. The presence 
of this bound state is indicated by the analytic structure of the two-body 
S-matrix between particles Aa and Ai,. The S-matrix must have a pole at 
iu':^ of the form 




(2.56) 



e 



— lU, 



,c 



'ab 
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CO/m 

Figure 2. Plots of the response function derived from the spin-spin correlator in the Ising model. 
In the first plot, the contributions from the two through six particle form factors are included. In 
addition, we plot the individual contributions of the four and six particle form factors. We see 
that the four and six particle form factors, to be seen, require a magnified scale. 

The coefficients forming the residue of the pole, 5^^, effectively measure the 
strength with which the two particles bind. 

We can formally indicate that Aa and Aj, have a bound state Ac via 

igabM^) = res5=o Aa{d + 6 + - iuts) ■ (2-57) 

Here u^^ = vr — u'^- and is the location of the pole in the two-body S- 
matrix between A), and Ac indicative of the formation of the bound state, 
Aa- The locations, u, of the various poles satisfy a simple relation 

27T = ul, + utc + ul-c. (2.58) 

which can be proven on the basis of kinematical considerations. If nia and 
mf, are the masses of the particles Aa and Af,, then the mass of the bound 
state, nic, satisfies = + + 2mamf, cos(u^^). 

The form of Eq. (2.57) indicates that we are able to relate n — 1 particle 
to n-particle form factors in the presence of bound states. In particular we 
have 

VeSS=of{9l, • • • , On-2, On-1 " iu^c, On-l + 5 + iM^g)ai,... ,a„_2,a,fe • (2.59) 

Thus with theories with bound states, consistency demands that form- factors 
come equipped with poles indicative of the bound states together with con- 
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straints upon the corresponding residues. 



2.7. Form Factors in Finite Temperature Correlators 

Up to this point we have only considered correlation functions at zero tem- 
perature. However the form factor approach can also be fruitful at finite 
temperatures. A finite temperature expansion of correlators is given in terms 
of a trace over the Boltzmann density matrix, 

^ E^,^e-f^^'r.(n,Sn\O{x,t)O{0,0)\n,Sn) 

Here the state, |n, s„), denotes a set of n-particles carrying quantum num- 
bers, {sn}- Inserting a resolution of the identity between the two fields as 
we did in the zero temperature case then leads us to a double sum, 

e~f^^'"{n,Sn\O{x,t)\m,Sm){m,Sni\O{0,0)\n,Sn) 

(2.61) 

We thus again have reduced the evaluation of the finite temperature corre- 
lator to the evaluation of form factors. 

To evaluate this expression for the correlator, we again focus on the as- 
sociated spectral function, G'^{k,(jj). In computing G'-^{k,uj), only terms in 
the form factor sum with a given energy, uj, and momentum, k, contribute 
to the sum, i.e. 



nsn 



msm 



xe-/5-^»n ^n, Sn\O{0, 0)|m, s™)(m, Sm\O{0, 0)|n, (2.62) 

as enforced by the presence of the two delta functions. For any uj,k, this 
dramatically reduces the number of matrix elements one must compute.'^ 
This reduction nonetheless usually leaves a difficult computation. However 
we can exploit the gapped nature of the spin chain to make the problem 
more tractable. Because the theory is gapped or massive (with gap. A), the 



''Here is simply the Fourier transform of G^{x,t), but similar considerations also apply to 
the corresponding retarded correlator. 
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correlator admits a low temperature expansion of the form, 

G'^{k,io) = ^a„(A:,w)e-"'^^. (2.63) 

n 

For the particular correlators of concern in this review and for the range of 
CO and k in which we are interested, each a„ is determined by a single matrix 
element. Because we can compute these matrix elements, we obtain an exact 
low temperature expansion. 

The ability to compute such an expansion should be compared with the 
approach taken by LeClair and Mussardo [199] (following Ref. [198]). These 
authors argued that it was possible to use the same form-factors we employ 
here to compute finite temperature correlators. However rather than directly 
evaluate individual terms in the sum appearing in Eq. (2.61), they first 
conjectured an ansatz involving a resummation of terms in the sum. This 
procedure was criticized in Ref. [265] (see also Ref. [52]). There it was 
argued that while this worked for the computation of one-point functions, 
it was problematic for two-point functions.^ Rather Ref [265] put forth 
the view that such problems should be attacked through the use of form 
factors computed against a thermalized vacuum [113,114,191]. However the 
counterexample cited in [265] , a computation involving interacting quantum 
Hall edge states, involved a gapless theory, and so is in a different class than 
the models considered in this review. (Without a gap, the low temperature 
expansion we consider above ceases to make sense.) The results of Section 
4 show that it is indeed possible, at least in certain cases, to make sense 
of the form-factor expansion of two point functions at finite temperature. 
But while one can make sense of this expansion, it will not be possible to 
directly compare computations in this review to the ansatz posited in [199]. 
Their ansatz as is applies only to diagonal theories where scattering does 
not permute internal quantum numbers, contrary, for example, to the case 
of the 0(3) non-linear sigma model considered in Section 4. 

2.7.1. Regularization of Form Factors 

Unlike in the zero temperature case, at finite temperature we must compute 
form factors with particles both to the right and to the left of the field, O, 

{A,^ {9m) ■■■A,, (0~i)O(O, 0)^„„ {On) ■■■Aa, (^l)). 



° That such a procedure is valid for one point functions is also supported by the study in Ref. [252] 
where form-factors are used to compute CFT correlators. 
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To understand such an object we must contend with the possibiUty that 
9i = 9j, tti = bj for some Prom the algebra of the Faddeev- 

Zamolodchikov operators (Eq. (2.4)), we know the commutation relations 
involve 5-functions, i.e. 

AaM)Al.i9,) = 27r6i9i - 9j)5a,b, + ■■■■ (2.64) 

It is crucial to include the contributions of the (5-functions to the correlators. 
To do so, we must understand the above form factor to equal 

{A,^ {9m) • • • {0i)O{0, 0) A„ {9n) ■■■Aa, (9,)) 

= Yl SA,AiSB,BABl\Al){B2\O{0,0)\A2)connected- (2-65) 

{aJ=AiUA2 
{6j}=SiUS2 

The sum in the above is over all possible subsets of {a^} and {bi}. 
The S-matrix Sa,Ai arises from the commutations necessary to rewrite 
^an(^Ti) ■ ■ ■ ^ai(^i)|0) ^ ^2^1 10) and similarly for Sb,Bi- The matrix el- 
ement, {Bi\Ai), is evaluated using the Faddeev-Zamolodchikov algebra. In 
this way (ill-defined) terms proportional to 5(0) are produced but which 
cancel similarly ill-defined terms arising from the evaluation of the partition 
function. 

The 'connected' form factor appearing in the above expression is to be 
understood as follows. Using crossing symmetry, the form factor can be 
rewritten as 

(I?2|O(0, 0)1^2) connected 

= {A,, i9,,)---A,, (4)0(0, 0)A„. {9jJ---A^. (^,J)connected 

= (0(0, 0)A^r {9jJ ■ ■ ■ A^r i9j,)A-^, {9i^ -m)--- Ay, {9i, - ztt)) connected 

Jq Jl t/j 11 

~ -^fe' ■■■b' a' ■■■a' ( '1 ~ ' ' ' ' ^ik ~ %l ' ' ' ' ' %q)connected) (2.66) 

n »fe ii jq 

where the last relation holds provided we do not have 9i = 9j, Qi = bj for 
any If this does occur we see from the annihilation pole axiom that the 
form factor is not well defined, having a pole at 9i = 9j. In such cases the 
form factor requires regulation. 

To regulate the form factor, we employ a scheme suggested by Balog [30] 
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and used by LeClair and Mussardo [199]. We define 

•^6' ■ b' a' ■■■a' ^ 'i ' ' ' ' H ^Ji ' ' ' ' ' ^jq) connected 
finite 



(2.67) 



'of l™n-^5' -6' a' -a'. (Oh-iT^ + im, ' ' ' , -ITT + irik, Oj,, ■ ' ' 

■ " T?i— »U XI Ik Jl Jq 

In taking the finite piece of f^, we discard terms proportional to rj^^ as 
well as terms proportional to rji/rjj. In this way the connected piece is in- 
dependent of the way the various limits ry.j are taken. Balog [30] has 
already used this prescription to compute one point functions and success- 
fully compare them to TBA calculations. In Ref. [30] it was argued that the 
delta functions leading to such terms arise from the use of infinite volume 
wavefunctions. If such wavefunctions are replaced instead with finite volume 
counterparts, the delta functions are regulated. For example, a pole in rj is 
changed as follows 



1 

IT] 



J e + iT] J e + iri' ^ ' 

where f{9) is some sharply peaked function about ^ = which in the infinite 
volume limit evolves into a (5-function. However the principal value of this 
regularized integral is zero. Thus discarding the pole terms is justified in this 
sense. For terms that are ratios of infinitesimals, Balog also demonstrates 
that such terms, once regularized, disappear in the infinite volume limit. 



3. Sine-Gordon Model and Spin-^ Quantum Magnets 

The first set of applications of massive integral QFTs we will discuss involves 
the quantum sine-Gordon model (SGM). It describes the low energy limit of 
a variety of gapped "deformations" of the spin-i Heisenberg chain. In order 
to fix notations and set the stage we first review the low-energy limit of the 
anisotropic spin-^ Heisenberg model. 



3.1. The Spin-^ Heisenberg Chain in a Uniform Magnetic 
Field 

The anisotropic spin-| Heisenberg chain in a uniform "longitudinal" mag- 
netic field H is defined by the Hamiltonian 

Wxxz = JY1 S^S^^, + + -HY^S^, (3.1) 

j j 
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and is exactly solvable by Bethe's ansatz [38]. In what follows we will con- 
sider only the region 



which corresponds to an "XY"-like exchange anisotropy. Ground state prop- 
erties, the excitation spectrum in the thermodynamic limit, thermodynamic 
properties and the large-distance asymptotics of correlation functions have 
been determined from the Bethe ansatz solution [150,164,179,286-288,318— 
320]. Dynamical correlation functions of the lattice model at zero tempera- 
ture and zero magnetic field have been evaluated in form of a spectral sum 
using the results of [163] in Refs. [1,46,47,165]. 

3.1.1. Excitation Spectrum of the Lattice Model 

The spectrum of low-lying excitations for the isotropic Heisenberg chain in 

zero field was derived from Bcthc's equations [38] in Rcfs. [104, 105]. It was 
shown that for even chain lengths the spectrum is given in terms of scattering 
states of an even number of gapless, elementary excitations carrying spin i, 
called "spinous" . There are no single-particle excitations: spinous occur only 
in pairs. The two-spinon continuum is shown in Fig. 3. 



-1< (5< 1, 



(3.2) 



4 



3 



■-2 
3 








7C 



Figure 3. Two-spinon continuum in the isotropic spin-^ Heisenberg chain {S = 1) in zero field. 
The triplet 5 = 1 and singlet 5 = states are degenerate. 



Application of magnetic field leads to a change in the ground state prop- 
erties. For any H > there is a finite magnetization per site (see below). 
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At a critical field 

He = J(l + S) , (3.3) 

the ground state becomes fully polarized. At He a phase transition be- 
tween a gapped, commensurate {H > He) phase and a gapless, incommen- 
surate {H < He) phase occurs, which is in the universality class of the 
commensurate-incommensurate transition [83,136,253,269]. In the follow- 
ing we will consider only the gapless regime H < He- 

In the presence of a magnetic field the spectrum is more complicated [166, 
229,294]. Now single-particle excitations exist,*^ but they have spectral gaps. 
The most relevant excitations at low energies are two-fold parametric and are 
similar in nature to the two-spinon excitations in zero field. We therefore call 
the elementary excitations giving rise to the gapless two-particle continua 
for H > spinous. ° The two-spinon excitations with 6S^ = and 6S^ = 1 
are shown in Fig. 4. The SS^ = 1 excitation remains gapless at P = tt/oo, 
whereas the SS^ = excitation becomes gapless at an incommensurate wave 
number for H > 0. Several variants of (3.1) are of considerable interest as 
they describe situations of direct experimental relevance. However, unlike 
the Heisenberg model itself these lattice models are not solvable by Bethe's 
ansatz. Interestingly, the low energy limits of these models are described by 
integrable massive quantum field theories, which allows for the calculation 
of the low-energy behavior of dynamical correlation functions. The starting 
point for all these considerations is the low-energy limit of the Heisenberg 
model (3.1), which we derive next following a well established procedure 
[5,133,134,215]. 



3.1.2. Continuum Limit 



The Heisenberg chain (3.1) is equivalent to a model of spinless fermions as 
can be shown by means of the Jordan-Wigner transformation 

1 



which maps (3.1) onto 



+5 



(3.4) 



1 

(3.5) 



They correspond to "string solutions" of the Bethe ansatz equations. 

In the Uterature various other terminologies have been used to denote these excitations. 
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Figure 4. Two-particle continua in the isotropic spin-i Heisenberg chain (5 = 1) in an applied 
field H = 0.5J. (top figure) States with SS^ = 1, which are probed by the transverse correlation 
functions (S^ (t)S-j~ (0)). (bottom figure) States with 55^ = relevant to longitudinal correlation 
functions (S^{t)Sf{0)). The boxes indicate the low-energy regime, which is described by the field 
theory approximation. 

It is instructive to analyze (3.5) for small values of 6. Denoting normal 
ordering by 

■.c]cj:= c]cj - {c]cj), (3.6) 

we rewrite (3.5) as 

- + J6{1 - 2{c]cj))^^c] Cj + const . (3.7) 
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Neglecting the interaction term, we may diagonalize the quadratic part 
of (3.7) by Fourier transform and determine the fermion number self- 
consistcntly. For small S the low-energy degrees of freedom are then found 
to occur in the vicinity of ztfe^, where 



2(5 



+ 



+ 0{6^). (3i 



Hence the continuum limit of (3.7) is obtained by keeping only Fourier modes 
with k ^kp, i.e. 



rv 1-1. 



(3.9) 



Substituting (3.9) into (3.7) we arrive at the low-energy fermion Hamiltonian 



n 



dx 



= v,J 

+JaoS J dx {R^R + L^L){R^R + L^L) 

+JaoS j dx {R^Le^'^> + h.c.)(i2^Le2<(^+"o) + h.c). (3.10) 

Here vp = Jaosin(fe^ao) and point splitting and normal ordering are im- 
plicit. Finally we may bosonize (3.10) using 



R\x) 



1 



exp 



, Lt(a 



exp ( ^f'i^) ) ' (3-11) 



where we choose the chiral bosons (p and to be normalized such that 

1 



(exp (iQ;(/7'(r, x)) exp (— ia(/?'(0))) 
(exp (iQ;(^'(r, x)) exp (— iQ;(^'(0))) 



{vT — ix)'^"^ 
1 



{vT + ix) 



2a2 



(3.12) 



The bosonized form of (3.10) is 

VF 



n = 

+ 



J 6 sin(2A;pao 



j dx\{i + ^^^^){d^^'f + {d.Q') 

^ ^ , dx d^^' -g f dx cos fV2$' - 4A;^[; 



Oq. 

2 



(3.13) 
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where g = J5ao/27r^ and the scalar field and its dual 6' are defined as 

= ip'{x) + (p'{x) , G'(a;) = ip'{x) - ^'{x). (3.14) 

We emphasize that the normalization condition (3.12) implies that the oper- 
ator cos(\/2*&') has dimension length"^ and hence the units in Eq. (3.13) 
are correct. The third term in (3.13) is important only in the vicinity 
of the isotropic limit (5 — ^ 1 in zero magnetic field. In this limit it is 
marginally irrelevant and gives rise to logarithmic corrections. These have 
been calculated by renormalization group improved perturbation theory in 
Refs. [14,15,34,35,212,214]. The second term in (3.13) is removed by shifting 
the field 



4V^()cos(A'2.ao) 



X . 



(3.15) 



It follows from (3.11) and (3.9) that the effect of this shift is to change the 
Fermi wave number to 

SH 



The canonical transformation 

1 

45 sin(A;^ao) 



irJao 



(3.16) 



1 + 



TT 



1 + 



46 sin(A:^ao) 



TT 



brings the Hamiltonian to the standard Gaussian form 



H = 



IGtt 



where 



vf(1 + 



2(5sin(A;FOo)\ 



TT 



(3.17) 



(3.18) 



(3.19) 



Here we have made use of the fact that we arc working in the small— 5 limit 
and all formulas are valid only to lowest order in 5. Importantly, the field 
G' and hence also G is compactified, i.e. 



e{x) ^ Q{x) + J . 



(3.20) 



This is a consequence of the U(l) symmetry cj - 
Hamiltonian (3.5), under which is invariant whereas G' transforms as 



e*"cj of the fermion 



32 F.H.L. Essler and R.M. Konik 



G' — > &' + 2\/2a- The identification (3.20) follows by setting a = 27r and 
using that the corresponding U{1) transformation is the identity. Similarly, 
the field and hence $ is compactified 

<^ = ^ + 87r/3. (3.21) 



It turns out that the above considerations can be generalized to the entire 
range (3.2) of 6. This may be done by comparing the 0{L~^) spectrum of 
the Heisenberg Hamiltonian calculated directly from the Bethe ansatz with 
the spectrum of the compactified boson (3.18). They are found to agree if 
the compactification radii fulfill (3.21), (3.20) and (3.29) and the velocity v 
and Fermi wave number kp arc chosen as follows [191]. Let us introduce 
a dressed energy e(A), dressed momentum p{X), dressed density p{X) and 
dressed charge Z{X) through the Predholm integral equations 



p{X) = — [ dn p{p) , 
«o Jo 

"A 



cosh 2A — cos 7 ' 

A 



P(A)- J ^^K{X-p) p{p) 



2 sin 7 



27r[cosh 2 A — cos 7] ' 
Z{X) - ^ K{X - p) Z(m) = 1 , (3.22) 

where 5 = cos (7) and the integral kernel is given by 

K{X) = -2 sin 27/(cosh 2A - cos 27) . (3.23) 
The integration boundary A is fixed by the condition 



£{±A) = . (3.24) 

The physical meaning of the various quantities is as follows: e{X) and p{X) 
are the energy and momentum of an elementary "spinon" excitation carrying 
spin 5^ = ib^. We note that spinous can only be excited in pairs. The 
magnetization per site in the ground state is given in terms of the ground 
state root density p{X) as 



(3.25) 
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The Fermi momentum is equal to 



27r IT 
kp = p{A) = — d\ p{X) = — 
«o Jo 0,0 



2 - <si: 



(3.26) 



where we have used that p{ — X) = p(A). The spin velocity is equal to the 
derivative of the spinon energy with respect to the momentum at the Fermi 
points 



V = 



56(A) 



dp{X) 



de{\)/d\ 



\=A 



2ttp{X) 



ao 



(3.27) 



X=A 



Finally, the dressed charge is related to /3 by 



V8Z{A) ' 



(3.28) 



In order to determine v and we solve (3.22) numerically, which is easily 
done to very high precision as the equations are linear. The results are shown 
in Fig. 5, 6 and 7. 




Figure 5. Spin velocity as a function of magnetic field for different values of S. Note that the 
spin velocity goes to zero as the magnetic field approaches He = J(l + 5). 



For zero magnetic field we have kp = 7r/2ao and 

-9 1 , „ Jan sin47r/3^ 

0^ = — arccos(-(5) , v = ^ . (3.29) 

47r ^ ' ' 2 1 - 4/32 ^ ' 
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Figure 6. "Fermi momentum" kp as a function of magnetic field for different values of S. We 
note that kp goes to zero for H —* He- 




Figure 7. Parameter /3 as a function of magnetic field for different values of S. We note that 
P^-^iorH-* He. 



For later convenience we define a dimensionless spin velocity by 

V 



V 



Jao 



(3.30) 



In what range of energies do we expect the continuum limit to provide a 
good description of the physics of the underlying lattice model? In zero field 
the answer is simply that the energy at which we probe the system must 
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be small compared to the cutoff of the continuum theory, which is equal to 
the exchange integral J. In the presence of a magnetic field, the cutoff is 
actually provided by H rather than J. 



3.1.3. Operators and their Norm,alizations 

For small values of S the lattice spin operators are related to the scalar fields 
Q through formulas (3.4), (3.9), (3.11) and (3.17). A straightforward 
calculation gives 



5-? 



d^^ - a{H)af sin - Ikpx ]+... , 
87r/3 \4/3 



Slr^A{H)a, « (-1) 



(3.31) 



The expressions (3.31) turn out to hold for —1<6<1 and the dots indicate 
operators with larger scaling dimensions, which contribute to the subleading 
terms in the large distance asymptotics of the spin correlation functions. In 
order to obtain the bosonized expression for S'^ it is necessary to symmetrize 
the Jordan-Wigner string operator as [5] 



kpap 



(3.32) 



and use that (c-c/,) = ^ . 

Interestingly, the coefficients a{H = 0), c{H = 0), A{H = 0) are known 
exactly [212,214]: 



I(a, b, c) 

a(0) 
c(0) 

^(0) 



.2-27]) 



2x/^r(^) 



exp 



sinh(6t) 



^dt 



2t [sinh(rf) cosh([l — ri]t) c 



2(T^^(i''''-0 • 



-2t 

(3.33) 
(3.34) 

(3.35) 



where rj = 4/3^. We note that the limit /3 — > ^, which corresponds to the 
isotropic spin-^ chain J — 1 is singular. This is because in this limit a 
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marginally irrelevant interaction of spin currents is present in the Hamilto- 
nian, which gives rise to logarithmic corrections in spin correlation functions. 
The isotropic limit 5 ^ 1 is discussed in detail in Refs. [14,212]. 

For finite magnetic fields the amplitudes have been determined numer- 
ically in [101, 140]. In Table 3.1.3 we list the results for the case of the 
isotropic spin-^ chain ((5 = 1 in (3.1)). 

Table 1. Amplitudes a and c, dimcnsionlcss spin velocity v, coupling /3 and 
magnetic field H as functions of the magnetization m for the isotropic Heisen- 
berg chain {S = 1). The amplitudes are determined numerically except for 
m = 0.5 where exact values arc shown. The figures in parentheses for a and c 
indicate the error on the last quoted digits. 
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3.1.4. Dynamical Spin Correlation Functions 

It is now a straightforward exercise ^ to determine the dynamical spin corre- 
lation functions at low energies, i.e. in the regions of the (a;, /c)-plane marked 
in Fig. 4. Using the bosonization dictionary (3.31) we find in imaginary time 
(r = it, X = jao) 



+e- 



2ik^x(j^ ,-<^0-^*) (^'^),<^0-^*) (0.0) 



A\H){-iy 
+c\H){-iy 



rj.2 _|_ y'ij-2 



il3e{T,x) ^2/36(0,0) 





2/32+1/8/32 1 


J.2 _j_ y2^2 ^ 





VT — IX 



2/32 



(3.36) 



{Tr 5|+i(r) 5f (0)) 



n 2 



+ a^(i7)ao'' (T^ sin 



(T, a,$(T,x)a,$(o,o)) 

<I>(t, x) 



1 



2a2 



_^a2(i7) cos{2kFx) 



4/3 
f r + ix 
VT — ix 



2k,x)sin{^)) 



4/3 



+ h.c. 



2-2 _|_ ^2^2 



(3.37) 



We see that the slowest decay of correlations occurs with an oscillating fac- 
tor (—1)-' in the transverse correlations and with factors exp{ziz2ikpx) in the 
longitudinal ones. Concomitantly their Fourier transforms will be most sin- 
gular at these wave numbers, which in turn makes k ^ ±2^;^ and ^ ~ ^ the 
most interesting regions in the Brillouin zone from an experimental point 
of view. The inelastic neutron scattering intensity (at zero temperature) is 



^We exclude the limit i? — > 0, <5 — > 1 from our discussion, see Refs. [14,212] for a discussion of the 
latter case. 
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proportional to 

I{iv,k) oc ^ (Sa-y - ^] y^' (c.',fc) , (3.38) 
where a, 7 = x, y, z and the dynamical structure factor (DSF) 8°'"' is defined 

by 

5"" {uj, k) = -- Im x"^(w, A;) , (3.39) 

TT 

/•oo 

x"^(u;,A;) = -z^ / e-^'=^«°+-*([5f+i(0, ^7(0)]) . (3.40) 



Here A; denotes the component of k along the chain direction and x°'"'{'^y k) 
are the components of the dynamical susceptibility. We note that in (3.38) 
only the symmetric combinations S^°''^\lv, k) = S"'^{u, k) + S"'"{u!, k) enter. 
Using a spectral representation one may show that for positive frequencies 
a; > 

/OO ^J. 

(3.41) 

We evaluate the retarded correlator in (3.40) by Fourier transforming the 

corresponding time-ordered correlation function in imaginary time and then 
analytically continuing to real frequencies. We arrive at the following results 
for the DSF in the vicinity of 2kF and ir/ao respectively 

5xx(^, ^ + g) = Syy{u, - + q) = ^ F{u;' - v\\ 2^^) ^ 
ao ao / 

S-\u:,2kF + q) = ^^ fL^-v\\^ 
^ ^ 4 V ^ ' 8/32 



T{a) 2vJ 



{2vJf 



(3.42) 



We see that the components of the DSF exhibit power-law singularities. 
The exponents vary with the applied field and also depend on the exchange 
anisotropy. However, the transverse components are always more singular 
than the longitudinal one. 
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3.2. Massive Perturbations and Window of Applicability of 
MIQFT 

As discussed above, the spin-i Heisenberg chain is quantum critical in the 
sense that the elementary spinon excitations are gapless and spin correlation 
functions decay as power laws. In what follows we consider several pertur- 
bations of the Heisenberg chain, in which excitations have a gap and spin 
correlation functions decay exponentially with distance. In particular, we 
will determine how the DSF is changed compared to the critical form (3.42). 
One issue to keep in mind that the field theory approach we employ has a 
limited window of applicability. Field theory becomes exact in particular 
scaling limits of the underlying lattice models. However, in experimentally 
relevant situations one usually is at some distance in parameter space from 
the scaling limit and a practical criterion is needed to judge, whether field 
theory will provide good approximations to the results for the underlying 
lattice model. A simple such criterion is that the spectral gap should be 
small compared to the cutoff set by the lattice model. The latter is simply 
min(H, J) for the Heisenberg chain in a magnetic field. Another possible 
criterion is to demand that the low energy excitations are relativistic up to 
some energy scale beyond that of concern (see, for example, Ref. [186]). 

3.3. Field-Induced Gap Problem 

In many materials with low-symmetry crystal structure the coupling of a 
uniform magnetic field to spin degrees of freedom is of a tensorial nature. In 
the case where there are two spins per unit cell the ^-tensor has generally 
both uniform and staggered components (see e.g. [245,246]) 

Wmag = -fiBHa J^idl + (-^Y 9lb]S'j ■ (3.43) 

j 

A low crystal symmetry also frequently results in the presence of a 
Dzyaloshinskii-Moriya (DM) [84,226] interaction 

^^DM = Y.^j (Si-1 X S,) . (3.44) 

i 

The direction of the DM vector Dj is constrained by the crystal symmetries. 
Of particular interest is the case when Dj has a staggered component, i.e. 
Dj = (— 1)''D. It was shown by Oshikawa and Affleck in [242] that in this 
case the application of a uniform magnetic field H leads to the generation of 
a staggered field perpendicular to the direction of H. The effects of having 
a staggered component of the ^-tensor and a staggered DM-interaction are 
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captured by the following model [242] 

H = Y^ JSj ■ Sj+i - HS^ - h{-iyS^ . (3.45) 
i 

Here the staggered field is proportional to the uniform field 

h = jH , (3.46) 

where 7 <C 1 depends on the details of the material under investigation. The 
model (3.45) has been applied successfully to experiments on several quasi- 
ID spin-| Heisenberg magnets: Copper Benzoate [Cu(C6D5COO)2 • 3D2O] 
[16,19,20,74,75], CDC [CuCls • 2((CD3)2SO)] [172], Copper-Pyrimidine 
[(PM • Cu(N03)2 • (H20)2)„ (PM=pyrimidine)] [110, 315, 316] and Yb4As3 
[181,243]. Theoretical calculations have been carried out for the excitation 
spectrum [94,207,242], the dynamical structure factor [94,101], the specific 
heat [95], the magnetic susceptibility [15,316] and the electron-spin reso- 
nance lineshape [244]. 

Bosonizing the Hamiltonian (3.45) by means of the identities (3.31) yields 
a sine-Gordon model [242] 

7^ = -|- J dx {d^^f + {d^ef -nih) J dxcos{P@). (3.47) 

where in our normalizations /x(/i) is a dimensionful quantity 

H{h) ^ h c{H) af^~\ (3.48) 

In our approach the velocity v and the parameter /? in (3.47) are determined 
for the lattice Hamiltonian (3.45) with h = 0. This is a good approximation 
as long as h is small compared to H and J. For given J and H , [3 \s 
then calculated from the Bcthc ansatz solution of the Heisenberg model in 
a uniform field, see Eq. (3.28) and Fig. 7, and 

/3 = /3 . (3.49) 

Similarly, the spin velocity is given by (3.27) and Fig. 5. The parameters 
entering the effective SGM are summarized in Table 3.1.3. 



3.3.1. Spectrum of the SGM and Quantum Numbers 
It is useful to define a parameter 
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The spectrum of the sine-Gordon model (3.47) in the relevant range of (3 
consists of a soliton-antisoliton doublet and several soliton-antisoliton bound 

states called "breathers." There are altogether [1/^] breathers, where [x] 
denotes the integer part of x. The breather gaps [66,103] can be determined 
for example from the Bethe ansatz solution of the massive Thirring model 
[39, 190] and are given by 



A„ = 2Asin(^), n = l,.... 



1 



(3.51) 



Here A is the gap for soliton and antisoliton. The parameter ^ depends on 
the applied uniform field H through (3. We introduce labels for the n*^ 
breather, s for the soliton and s for the antisoliton. A convenient basis for the 
SGM (3.47) can be constructed in terms of scattering states of (anti)solitons 
and breathers by means of the Faddeev-Zamolodchikov algebra (see Section 
2.2). The action of the charge conjugation operator C on these basis states 
states follows from 

C|0) = |0), 

cAl{e)c-' = AUe), 

c4^(^)C-i = (-l)"4^(e). (3.52) 

The topological charge 



Q 



B f°° 

^ / dx d^Q , (3.53) 



is a conserved quantity in the sine-Gordon model. We will use conventions in 
which soliton and antisoliton have topological charges —1 and 1 respectively. 
Breathers have topological charge zero and so are neutral. 

The soliton gap as a function of the parameters H and h of the underlying 
spin chain was determined in Ref. [15] in the regime A <^ H, where 

-{1+0/2 

(3.54) 



J 



with B = 0.422169. Equation (3.54) is applicable as long as H is sufficiently 
smaller than J or more precisely as long as the magnetization is small. In 
the derivation of (3.54) both the magnetic field dependences of spin velocity 
and the normalization constants a and c of the spin operators in (3.31) have 
been neglected, which is justified in weak fields. An additional complication 
in weak fields is the presence of a only slightly irrelevant perturbation to the 
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free bosonic effective Hamiltonian describing the low energy physics of the 
Heisenberg chain (see Eq. (3.13)). By taking this perturbation into account, 
renormahzation group improved perturbation theory essentially leads to the 
result (3.54). 




For magnetic fields comparable to J it is necessary to take into account 
the magnetic field dependences of the spin velocity and the normalization 
c{H) of the spin operator, whereas the effects of the irrelevant perturbations 
to the free boson Hamiltonian may be neglected. Using the results of Ref. 
[334] we obtain the following expression for the gap in the regime of H 
comparable to J (but still h <^ J) 

(1+0/2 



A 

7 



r(|) 



2v{H) 



c{h)ti r(^) h 
2v{H) r(4^) J 



(3.55) 



Here v = v/{Jao) is the "dimensionless spin velocity". 



3.3.2. Thermodynamics 

The thermodynamics of the SGM is most efficiently studied [77] via the 
recently developed Thermal Bethe Ansatz approach [76,179,182,284,289], 
which circumvents problems associated with solving the infinite number of 
coupled nonlinear integral equations that emerge in the standard approach 
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based on the string hypothesis [111, 112] (note that the coupUng constant 
(3 in our problem is a continuously varying quantity and no truncation to a 

finite number of coupled equations is possible). It was shown in [77] that the 
free energy of the SGM can be expressed in terms of the solution of a single 
nonlinear integral equation for the complex quantity e(9) (we set ks = 1) 

s{e) = -i- smh{e + ir]') - / dO'CoiO - 9') In (l + exp [-e{9')] ) 

J — oo 

/oo 
de'Goie -e' + 2ir/') In (1 + exp [-e{e')] ) , (3.56) 
-oo 



-oo 

where A is the soliton mass and 



dw cos(2a;g/7r)sinh(a;(^-l)) 
^'^^>-Jo ^ sinh(a;Ocosh(a;) " ^^■^^> 

The free energy density is given by 

/(T) = Im / de smh{e + ir]') In 1 + e'^^^ . (3.58) 

As we are interested in the attractive regime of the SGM we have 

< V < . (3.59) 

Note that the free energy does not depend on the value of rj' as long as 
it is chosen in the interval (3.59). The set (3.56) of two coupled nonlinear 
integral equations is solved by iteration. For T — > the first iterations can 
be calculated analytically and the corresponding contributions to the free 
energy density are seen to be of the form 



^ _2TA g (-iTi^, _ (^) ^ . . . , (3,0) 

TTV ^-^ n \ T J TTV \ T J 

n=l ^ ^ ^ ^ 

where Ai = 2A sin ^ is the mass of the first breather and Ki is a modified 
Bessel function. The first term is the contribution of soliton-antisoliton scat- 
tering states to the free energy, whereas the second term is the contribution 
of the first breather. Both terms have the form characteristic of massive 
relativistic bosons. The contributions of the heavier breathers are found in 
higher orders of the iterative procedure employed in solving (3.56). The 
specific heat is obtained from the free energy 

C = T^^. (3.61) 
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At low temperatures it is found to be of the form 



C 



[i/C] 

a=0 



1 + 5ao)A 



2-KV 



T 



3 / r 

4 1 a: 



A, 
T 



exp 



A. 



(3.62) 



where Aq = A. In order to compare theoretical predictions based on the 
SGM with the specific heat data of [75] we need the free energy at "inter- 
mediate" temperatures and thus have to resort to a numerical solution of 
(3.56) by iteration. A detailed comparison of the specific heat calculated 
in the framework of the low-energy effective sine-Gordon theory to experi- 
mental data on Copper Benzoate was carried out in Ref. [95] . In Fig. 9 we 
show a fit of the specific heat calculated from (3.61) [95] to data taken for a 
particular orientation of the magnetic field on Copper Benzoate [75]. 




0.0 



0.0 




Figure 9. Specific heat as a function of temperature for fields of H=1T and H=2.5T applied along 
the c" axis. 



3.3.3. Normalization of Form Factors in the Sine-Gordon Model 

The short-distance behavior of correlation function in the sine-Gordon model 
is governed by the operator product expansion (in imaginary time) 

hm e^"*(^'^) e^"'*(o.o) = + t;VY""' e^("+"')*(0'0). (3.63) 
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Given the normalization implied by (3.63) the large-distance asymptotics is § 
lim (e^"*(^'^) e^"'*(°'°)> = <e'"*) (e^"'*) . (3.64) 

where the one-point functions have been determined in [211] 



Ga = a^" (e 



2a2 , jQ,e(0,0)\ 



2a^ 



X exp 



Jo T 



A v? r((i + o/2) \ 
Jv 2 r(e/2) ; 

sinh^(a/3t) 
^2sinh(f t) sinh(l) cosh[i^t] 



2a2e-* 



(3.65) 



Equation (3.64) can be used to fix the overall constant factor in the form 
factor expansion of correlation functions. 



3.3.4. Dynamical Structure Factor 

Let us now turn to the issue how the staggered field modifies the dynamical 
structure factor. Without a staggered field one simply has the result (3.42). 
At low energies the field theory description can be used to determine the 
dynamical spin-spin correlation functions. We substitute the decomposition 
(3.31) into the expression (3.40) for the dynamical susceptibility. Neglecting 
the contributions of rapidly oscillating terms to the integrals, we find that 



X 



/"oo poo 

{io, — + g) ~ -ic^{H)af "W dx dt e^'^*-^«^ 

O-O J-oo Jo 

X ([cos(/3e(i,x)),cos(/3e(0)]), 

/-oo /-oo 

X^^ico, — + g) ~ -ic^{H)af dx dt e^'^*-^«^ 

«0 J-oo Jo 

X ([sin(/3e(t,x)),sin(/3G(0)]), 

J-oo Jo 

exp — — — ,sm ' 



(3.66) 



(3.67) 



4/3 



4/3 



(3.68) 



8 Here we do not consider the connected correlation function. 
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Using a spectral representation in terms of scattering states of solitons, an- 
tisolitons and breathers (2.2) we have e.g. 

27r 2/rr\ 4/3^-1 v^v^ [ dOi . . . dOn 



X-(a;,^+g) = 2.c^(iI)af-^5:/ 



, (27r)"n! 

n=l £i " ^ ' 

x|(O|cos/3e(O)|0„,...,9i),„,...,,j2 

Y.j -r sinh Oj) 5{q + J2j sinh ( 



UJ 



- X) ■ Aej cosh 6*^ + ie uj + Y^^ cosh 6*^ + ie 



(3.69) 



Here the indices £j G {s, s, Bi,B2, . . . , and = = A, Ab„ = A„. 

Many of the form factors actually vanish for the operators we are interested 
in here. Using the transformation property of the scalar field Q under charge 
conjugation 

CeC-^ = -e, (3.70) 
we observe see that by virtue of (3.52) 

(o|cos(/?e(o))|0)B,„_i=o, 

(0|sin(/3e(0))|^)B,„=0, 
(O|sin(/3e(O))|02,0i)B„,B„=O, (3.71) 

and so on. On the other hand, the ±2A;ir components of carry topological 
charge ±1. This follows from their commutators with the operator (3.53) 



Ti 

Q,exp I ±i— $ 



As a result the only non-vanishing form factors are of the form 

(O|e-**('^)|0),, (O|e-^*(°>i,02).,B^, {0\e-i^^'^\ei,92,93)s,s,s, (3.73) 

and so on. Combining (3.71), (3.73) and (3.68) we can identify which excited 
states inelastic neutron scattering experiments probe at different momentum 
transfers. This is summarized in Fig. 10. 

The form factors required to evaluate the first few terms in the spec- 
tral representation (3.69) have been determined in [24 26,210,213,278,280]. 
Determining the breather form factors via the bootstrap axiom (2.59) one 
obtains [94] 
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CO ii 




2kFao 71 2n-2kFao ^^o 

Figure 10. Schematic structure of the lowest-energy excited states relevant to neutron scattering 
experiments. Soliton and antisoliton occur in the vicinity of the incommensurate wave numbers 
2kF, — kp) and are seen in S^^, whereas the first breather Bi occurs in the vicinity of Tr/ao 
and contributes to S^^ . At higher energies further breather bound states as well as multi-particle 
scattering continua are present. 



- + q)= C,{H) Rei 27r ^ Sis'" - Ai J 



|^cos(/3e)[g(A„Ai,.)]P jF-(^Q)[g(A,A,.)]P 1 

s - 4Af - 4A2 ■■■]' ^ ■ ' 



where 



6{mi,m2,s) = arccosh I ) , (3.75) 

\ 2mim2 / 

and the overall normalization is given by 

C^(iJ) = - c\H)v{H)J gl . (3.76) 

Here is given by (3.65). As wc arc dealing with a two point function of a 
Lorentz scalar, the susceptibility depends only on the Mandelstam variable 

s = Va;2 - i;2g2 (3 77) 

The first terms in (3.74) correspond to single-particle breather states and 
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are given by [94] 



2(sin27re)^ 
Z2 = ■ — - — exp 



COtTT^ 



dx (sinh 2^x)^ sinh a:(l — ^) 
X cosh a; sinh 2a; sinh^x 



(3.78) 



^4 = 



2(sin 47r^)^ 



(cot 7r02(cot 37r^/2)2 cot(27rO 



X exp 



, r°° dx 
Jo X 



dx (sinh4^x)'^ sinhx(l — ^) 
cosh a; sinh2x sinh^a; 



(3.79) 



The sohton-antisohton contribution is 

(2cot7r^/2 sinh6')2 cosh e"/^ + cos tt/^ 



F 



cos(/3e) 



E{e) = exp 



^2 cosh2^/^-cos27r/e^ 
°° dx [cosh 2x cos 2a;0/7r — 1] sinha;(^ — 1) 



[I 



x 



coshx sinh2x sinh^a; 
and the BiBi breather-breather contribution is 

4 



(3.80) 
(3.81) 



(sinh 0)2 



2cosyW2sin — 



exp 



dx 



27r sin x 



(sinh 6') 2 + (sin7r^)2 



exp 



Jo 



dx cos ^ sinh ^a: sinh a;(l + ^ 



sinh 2x cosh x 



. (3.82) 



The next most important contribution (i.e. the term with the lowest thresh- 
old) in the expansion (3.74) comes from B^B^ breather-breather states. It 
will contribute at energies larger than Ai + A3, where Ai^3 are given by 
(3.51). In Fig. 11 wc plot S^^{uj, ^ + ?) as a function of s = y^cj^ _|_ ^2^2 
for two values of the magnetization m. We see that the coherent particle 
peak due to the second breather carries most of the spectral weight. For 
m = 0.1 the contribution due to the soliton-antisoliton continuum is consid- 
erably sharper above the threshold than for m = 0.03. This is because the ss 
continuum has just given birth to the fourth breather. The latter has a very 
small binding energy and carries very little spectral weight. In general the 
ss continuum becomes singular at the threshold, whenever a new breather 
peak splits off from it. 

Let us now turn to the yy component of the dynamical susceptibility. As 
the U(l) spin rotational symmetry around the axis of the uniform field H is 
broken by the staggered field h, must be different from x™(a;,g). The 
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Figure 11. xx component of the dynamical structure factor for magnetizations (a) m = 0.03 and 
(b) m = 0.1. The delta-functions corresponding to the second breather have been broadened by a 
Lorentzian to make them visible. 



form factor expansion yields 



s^^uJ,— 



+ q) = C,{H) Re| 27r ^ Z2n-i 6{s'' - A^.^^ 

I n=l 



sVs'^ - 4A^ v/(s2-A^-A^)2-4A^A^ ■■■]' ^ ' ' 
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Here the breather contributions are 



8(cos^)4 



cot 



exp 



^ °° dx sinhc^a; sinhx(l — ^) 
X cosh a; sinh2a; 



(3.84) 



Z3 



4sin(37rO(sin2|S)2 



dx (sinhS^x)^ sinhx(l — ^) 



(cotTrO^ 

r .c 

X exp -2 / , . , „ . , ^ 

I Jo X coshx sinh2x smh^x 

The sohton-antisohton contribution is 



(3.85) 



in(/3G)^^^|2 _ (2cot7r^/2 sinhg)^ cosh6'/g - cos7r/g 



cosh 26*/^ -cos 27r/^ 



E{e), (3.86) 



where E{6) is given by (3.81). Finally we take into account the B1B2 
breather-breather state which contributes 



^sin(/3e)^^-j|2 ^ taUTT^ 

dx X 

X exp —6 / - — : 

ZTT Sin X 



sinh2(0-i^) 



sinh2(0-z^) + sin2 7re 
1 

1 + 



2cosyy2siny 



X exp 



4cos^(cosh(9 + cos^) 



dx [cos ^ sinh^x + coshx^ sinh(^)] sinh(a;^^) 
X sinhxcosh(l) 



. (3.87) 



The next most important contribution to (3.83) is due to BiBiBi three 
breather states with a threshold at 3Ai. In Fig. 12 we plot S^^^lo, ^ + g) as 

a function of s = y^oj'^ + v'^q^ for two values of the magnetization m. The 
two coherent single-particle peaks due to the first and third breathers carry 
most of the spectral weight. The incoherent scattering continua are always 

weak and become less pronounced as the magnetization increases. 

The leading contributions to the longitudinal structure factor at the in- 
commensurate wave numbers 2kp, — — 2kp are due to soliton and anti- 

soliton single-particle states. The form factor (0|e ^/s -"l^)^ ig a rapidity- 
independent constant (as the operator is a Lorentz scalar), which has been 
determined by Lukyanov and Zamolodchikov Ref. [213]. A simple calculation 
then gives [101] 



S^^ (w, 2kF + q)= C,(H) 27T Zi 6 ((j2 _ (vqf _ A^) , 



(3.^ 
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Figure 12. yy component of the dynamical structure factor for magnetizations (a) m = 0.03 
and (b) m = 0.1. The delta- functions corresponding to the first and third breathers have been 
broadened by a Lorentzian to make them visible. 



where C,{H) = ^ a^H) 



2n 




X exp 



exp[-(l + e)nt] - 1 



°° dt 

t \ 2 smh(^t) sinh[(l + ^t] cosh(t) 



+ 



n 



2smh(tO 4/32 



(3.89) 
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Here the constants Ci,2 are given by [213] 

= exp (- r - 1^'^ - I'' I . (3.90) 



\ JO t sinh(2t) sinh(^t) cosh(t)^ 

\ Jo t sinh(2t) sinh(^t) y ^ ^ 

As was pointed out in Ref. [15], at if = the low-energy effective theory of 
(3.45) is SU(2) symmetric. In our notations this imphes that 

1™ 2/m = 2 Hm ^ ^ . 3.92 
Equation (3.92) is easily verified numerically. 



3.4. Transverse-Field Model 

Some quasi-lD spin-| antiferromagnets are described by Heisenberg models 
with exchange anisotropy. One example is CS2C0CI4 [171], where the Co 
spin- 1 multiplets are split by a strong single- ion anisotropy. Projecting to the 
low-lying spin-| doublet leads to a Hamiltonian of the form (3.93) with ^ = j 
and H = 0. Anisotropic Heisenberg models can also be used to describe 
spin chains with Dzyaloshinskii-Moriya [84, 226] interaction. The bulk DM 
interaction can be removed by a gauge transformation (local rotation of the 
coordinate system in spin space) ; the resulting Hamiltonian has an exchange 
anisotropy and twisted boundary conditions (in a ring geometry) [17,216]. 

When applying a magnetic field to an anisotropic Heisenberg model we 
have to distinguish between two cases: 

(i) the magnetic field is along the direction singled out by the anisotropy 
and leaves the spin rotational U{1) symmetry of the Hamiltonian 
unchanged; 

(ii) the magnetic field is at an angle to the direction singled out by 
the anisotropy and breaks the spin rotational U{1) symmetry of the 
Hamiltonian. 

In case (i) the model remains exactly solvable and has been discussed in 
section 3.1. The simplest realization of case (ii) is when the field is perpen- 
dicular to the anisotropy 

^ = E -^i^jSj+i + ^J^J+i + ^SjS^+i] - , (3.93) 
j 
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where S < 1. The model (3.93) has been studied in Refs. [51,55,78,79]. 
At low fields there are both uniform and staggered ordered moments and 

excitations are gapped. At a critical field Hc{5) there is quantum phase 
transition in the universality class of the 2D Ising model, at which the stag- 
gered magnetization vanishes. In high fields only the uniform magnetization 
remains and excitations are again gapped. The simplest case to deal with by 
field theory methods is the one where the magnetic field H is much stronger 
than the anisotropy J(l — 6), but is smaller than Hc{5). The field theory 
limit in this case has been studied in Refs. [55,95,244]. It is useful to change 
coordinate system and rewrite the Hamiltonian (3.93) as 

:^zxx,H = -^^i ■ - ^^1 + -^(^ - 1) •^J^l+i = Wo + . (3.94) 

3 3 

Our strategy in the regime H S> J(l — 5) is to bosonize at the point (5 = 1 at 
a finite value for the magnetization per site and then to treat the exchange 
anisotropy Tii as a perturbation. The continuum limit of TIq is constructed 
in section (3.1.2) and is given by a compactified boson (3.18). The per- 
turbing Hamiltonian 7ii can then be bosonized using (3.31). Fusion of the 
y-components of the staggered magnetizations gives a contribution 

Oj = S]S]^-^ — > C af cos(2/3e) + . . . . (3.95) 

In addition there is a small marginal contribution that shifts the compactifi- 
cation radius. For simplicity we neglect it here. Putting everything together 
we conclude that at low energies compared to the scale set by the applied 
field H < 2J, the effective Hamiltonian is given by a SGM 

H = [(^x$)' + {dccQf] - m(<5) cos(/3e) , (3.96) 

where 

I3 = 2P . (3.97) 

The cosine term in the SGM is relevant and generates a spectral gap. As 
P > (see Fig. 7) the SGM is in the repulsive regime and the the spectrum 
consists of soliton and antisoliton only. No breather bound states exist. 
The magnetic field dependence enters both via the prefactor and via the 
i?-dependence of /?. In order to calculate the prefactor as well as quantities 
like the magnetization we need to know the normalization C in (3.95) in 
the Heiscnberg chain in a field, i.e. the Hamiltonian (3.94) with 5=1. In 
Ref. [141] C was estimated numerically from the large distance asymptotics 
of an appropriately chosen four-point function and found to be very small. 
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An independent method for determining C from the staggered magnetization 
was suggested in [55] . The gap is given by 



A 

7 




(3.98) 



where v is the dimensionless spin velocity v = -j^. As we have mentioned 
before, in the low field phase the staggered magnetization in x-direction is 
nonzero in the presence of a transverse field. One has 

ii-irS^) = c{H) <cos(^G)) = c{H)g,_ , (3.99) 

where Qis is given by (3.65) (which in turn depends upon -y). Knowledge 

of the staggered magnetization as a function of H (say from numerics) then 
will yield C. 

Next we turn to the low-energy behavior of dynamical correlation func- 
tions. We start with the transverse correlations and concentrate on momenta 
close to vr/aQ. The staggered components of the spin operators are given by 
(3.31) and calculating their correlation functions reduces to the calculation of 
appropriate correlators in the SGM (3.96). The leading contribution to the 
yy-component of the dynamical structure factor is due to soliton-antisoliton 
two-particle intermediate states. One finds [55] 

^yy. ^ , _ vJc\H)glf^ _ ^(g(A, A, s))eH (i - 2) 
^"^'ao^'^ TT^^A^ . cosh(^(^)+cos(f) 

-I- contributions from 4, 6, . . . particles. (3.100) 

where 9h{x) is the Heaviside function, ^(A,A,s), Qa and E{9) are given 
by(3.75), (3.65) and (3.81) respectively and s = \J uP- — v'^q^. The analogous 
result for S^{ijj, ^ + q) \s 

^xx, ^, X ^J<^\H)Gl^^ - 4A2 E(0(A, A, s))eH (i - 2) 
^ {LO, h q) 



ao ^' vre2A2 s cosh ( ^(^) - cos (f) 

-I- contributions from 4, 6, . . . particles. (3.101) 

Finally we determine the longitudinal structure factor S^^{uj, k) in the vicin- 
ity of k = 2kF and k = ^ — 2kF- Recalling that the sine-Gordon coupling 
constant is given by (3.97), we see from the bosonized expressions for the 
lattice spin operators (3.31) that we need to consider two-point functions of 
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exp in the SGM (3.96). By (3.72) these operators have topological 

charge =f2 respectively. Hence the leading contribution to the zz-componcnt 
of the dynamical structure factor is due to intermediate states with two soli- 
tons or two antisolitons. The corresponding form factors are readily calcu- 
lated from the axioms (2.9)-(2.16). The normalization has been determined 
in Ref. [213]. A short calculation then gives 



2., + ,) = "'■;^f^y? ^IBE Ele(A. A, (i - 2) 

-|- contributions from 4, 6, . . . particles, (3.102) 

where Z2, E{9), Ci are given by (3.89), (3.81) and (3.90) respectively. We see 
that the leading contributions to the various components of the dynamical 
structure factor are always due to two particles. The structure factor is 
entirely incoherent and always vanishes at the threshold (as (3 is always 
strictly larger than -^), which occurs at a; = 2A. 

3.5. Dimerized Chain 

In actual materials the spin degrees of freedom arc always coupled to 
phonons. In many cases this coupling is small and can be ignored, par- 
ticularly at low temperatures. However, in quasi-lD compounds such as 
CuGeOs the magneto-elastic coupling leads to a phase transition, known as 
the spin-Peierls transition. In the low-temperature phase the lattice is then 
dimerized. This is turn leads to a modulation of the Heisenberg exchange 
and an appropriate model Hamiltonian is 

7i = J^[l-(-l)"5]Sn-S„+i . (3.103) 

n 

For 5^1 one may study (3.103) as a perturbation of the isotropic Heisen- 
berg chain in the field theory regime. The bosonized expression for the 
staggered energy density e„ = (— 1)"S„ • S„+i can be determined from the 
bosonized expressions of the spin operators (3.31) for zero magnetic field. 
The relevant piece arises from the fusion of the smooth with the staggered 
components and is given by [233] 



elx 



cos(^$(x)). (3.104) 



The value of the constant T> is presently not known analytically. In the 
aforementioned fusion procedure subleading (i.e. very irrelevant) terms in 
the bosonized expressions for the spin operators (3.31) contribute to the 
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value of P, so that it cannot be simply calculated from the known expressions 
(3.33), (3.34), (3.35) for a(0), c(0), and ^(0). However, it was found in [241] 
by comparing to numerical results that V is well approximated by 

D.^(f)'. (3.105) 

Bosonizing the 6 = part of (3.103) and then using (3.104) we find that the 
field theory limit of the dimerized chain is 

n= Jdxl^[{d^^f+{d^@f]-ficos{^^) I, (3.106) 



1 



where n = '^'DJ6 and where we have dropped a marginally irrelevant 
interaction of spin currents present in the field theory limit of the isotropic 
Heisenberg chain ^ 



Hmsrg. =jjdx C0S($) + ^ [{d^Qf - {d^^f 



(3.107) 



Hence the dimerized chain at small A is described by a SGM with /3 = ^. It 
is important to note that the field theory respects the spin rotational SU(2) 
symmetry of the original lattice Hamiltonian: the SGM (3.106) is SU(2) 
invariant [135]. This non-obvious fact is most easily understood using non- 
Abelian bosonization, see section VI of Ref. [15]. The spectrum of the SGM 
(3.106) consists of four particles: soliton, antisoliton and a light breather 
forming a SU(2) triplet with gap A and a heavy breather with gap \/3A. 
The latter is a SU(2) singlet. In contrast to the field induced gap problem, 
the nonlinear term is now the cosine of the bosonic field <I> rather than the 
dual field G. Hence the topological charge is the dual of (3.53) for /3 = |, 
which is related to the z-component of the total spin by (3.31) 

Q = — y dxd,^ = ^S]. (3.108) 

The relation of the topological charge to the z-component of the spin agrees 
with our assertion that the breathers have S"^ = and soliton/antisoliton 
spin ^1. A consequence of the SU(2) symmetry is that the following two- 



^ The bare value of 7 has been estimated e.g. in [6] 
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point functions are identical in the SGM (3.106) 

(cos(ie(T, x)) cos(ie(0, 0)) ) = (sin(ie(r, x)) sin(ie(0, 0)) > 

= (sin(i$(r, x)) sin(^$(0, 0)) ) . (3.109) 



Equations (3.109) imply various relations between form factors, which can 
be verified by direct calculation (see e.g. [93]). In order to determine the 
dynamical structure factor it is sufficient to consider the two-point function 
of the operator sin(^$) in the SGM (3.106). Apart from the overall constant 
factor, the latter is the same as the result for the field induced gap problem 
(3.83) with H set to zero ' and an appropriate expression substituted for the 
soliton gap. Using (3.105) and taking the marginally irrelevant interaction of 
spin currents into account in rcnormalization group improved perturbation 
theory, the soliton gap is found to be [241] 



J 



r(i) 



3 r(|) 



(3.110) 



l + §|A|ln|^ 



36125 I 



where A ~ —0.22. The dynamical structure factor is then given by 



ao 



C Re|27rZi ^(s^ - A^) 



^sin(i0).... . ...r, , sin, j; 

^ 2|F+_^^^[g(A,A,g)]p ^ 2|Fi2 



[e(A,V3A,s) 



sVs^ - 4A2 



V(s2 - 4A2)2 - 12A4 



(3.111) 



where Zi, e{a,b,c), Ff'l{9) and Ff^\9) are given by (3.84), (3.75), (3.86) 
and (3.87) with = 5 (C = g) respectively. We have denoted the analog of 
the overall factor C{H) in (3.83) by C. We conjecture that 



C 



- — -ryi > 

(27r)2 2 



(3.112) 



where Qi is given by (3.65) with /3 = ^ and i) = f . We plot the DSF (3.111) 
in Fig. 13. The DSF is dominated by the coherent triplet modes with gap 
A given by (3.110). Starting at twice the triplet gap there is an incoher- 
ent scattering continuum. The latter is singular at its threshold 2 A, as is 
the B1B2 scattering continuum which occurs above A + \/3A. The heavy 



' At if = we have /B = ^, v = ^ as required. 
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Figure 13. Dynamical structure factor of the dimerized chain. Most of the spectral weight is 
carried by the coherent triplet excitations, which are identified as soliton, antisoliton and the first 
breather in the SGM. The delta function corresponding to the coherent triplet excitation has been 
broadened by a Lorentzian to make it visible. 

breather B2 has total spin zero and is invisible to neutron scattering experi- 
ments. However, it contributes to the two-point function of energy densities 
e{x) and hence should be observable in Raman scattering experiments. An 
expression for the two-point function of energy densities was derived in [93]. 

3.6. Quasi-ID Spin-^ Antiferromagnets in the Ordered 
Phase 

A generic feature of quasi-lD spin-^ materials is that at low temperatures 
they develop long range Neel order. This is a consequence of the residual 
coupling between the ID chains as spontaneous symmetry breaking is for- 
bidden in one spatial dimension. At low energies in the ordered phase the 
physics is dominated by Goldstone modes: if the Neel order is along the z- 
axis in spin space, there will be two transversely (in the xy-plane) polarized 
spinwave modes (one mode each for the "broken" symmetry generators 
and S^). Their quantum numbers are = ±1. In addition to the transverse 
Goldstone modes there a priori also must be gapped longitudinal (polarized 
along the z-axis) excitations, which correspond to fluctuations in the size of 
the ordered moment. 

On the other hand it is clear that at temperatures or energies that are 
large compared to the interchain coupling one has to recover the ID spin 
chain physics. 



Integrable QFTs in condensed matter physics 



59 



A long-standing question is how the system crosses over from the 3D 
physics of the ordered phase to the ID physics of the chains, as one varies 
the energy scale at which it is probed. This issue was addressed theoretically 
in Refs. [13, 93, 271, 342] and experimentally in [193] for KCuFa (see also 
[266,297,298]) and in [340,342] for BaCu2Si207. 

Our starting point is the Hamiltonian for an ensemble of weakly coupled 
chains. For dcfinitcness we discuss the form appropriate for BaCu2Si207 
[169,336,337], which reads 

+J3 ^i,j,n • (Si+lj-|-l,n + Si+l,j-l,n) ■ (3.113) 

where = -0.460(7) meV, Jy = 0.200(6) meV, 2J3 = 0.152(7) meV and 
J = 24.1 meV. The Fourier transform of the inter-chain coupling is defined 
as 

J'(p) = Jx cos{px) + Jy cos{py) + J3 [cos{px + Py) + cos{px - Py)] . (3.114) 

As some of the exchange constants are positive and some negative it is con- 
venient to introduce new spin variables 5" : 

^i,j,n = ^i,j,n ' ^i,j,n = i~^y ^i,j,n ' Ot = y^Z . (3.115) 

The transformation (3.115) leaves the Hamiltonians of the ID chains invari- 
ant, but flips the signs of Jy and J3 in the interaction of the y and z compo- 
nents of the spin operators in H' . In the new spin variables the interchain 
coupling is ferromagnetic so that we may take the staggered magnetization 
at T = to be 

{St,j,n)=5a^z {-IT mo. (3.116) 

3.6.1. Mean-Field Theory 

As a first step we take the long range Nccl order into account in a self- 
consistent mean-field approximation [268,271]. Writing 

^t,j,n = ^^t,j,n) + ^^t,j,n 1 (3.117) 
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where dS'j^j ^ denote (smah) fluctuations around the expectation value, and 
then substituting (3.117) in H' we obtain 

-f^MF = X X "^^^'i'" ■ - /l(-l)" Sfj^j^ , 

i,j n 

h = 2{Jy -Jx + 2J3) mo = J'mo . (3.118) 

The Hamiltonian (3.118) describes an ensemble of uncoupled spin-i Heisen- 
berg chains in a staggered magnetic field and is a special case of the field- 
induced gap problem we studied in section 3.3. At low energies, the mean- 
field theory Hamiltonian (3.118) reduces to a SGM for the Bose field $ 
(rather than the dual field 6 as was the case in section 3.3) 

■Hmf = jdxl^ [{dx^? + {d^ef] -^icos(^ $) i, (3.119) 



where fj, = cha^ ^ with some constant c and J 

TrJao 



(3.120) 



2 

The form (3.119) is obtained as follows: bosonizing the Hamiltonian (3.118) 
in the absence of the staggered field gives a Gaussian model of the form 
(3.18), where the field is compactified according to $ = $ + 47r. In addi- 
tion there is a marginally irrelevant interaction of spin currents of the form 
(3.107). Bosonizing the staggered field term by means of (3.31) gives a sine- 
Gordon like Hamiltonian, but with a sine potential rather than a cosine. 
Shifting the Bose field 

$ — (3.121) 

while keeping the dual field unchanged and dropping the marginally irrel- 
evant term then leads to (3.119). The shift (3.121) changes the bosonized 
form of the lattice spin operator S"? to 

SI -'-^d^^ + ca\ {-ly cos . (3.122) 

At the particular value of /3, the spectrum is formed by scattering states 
of four particles: soliton, antisoliton and the first breather Bi all with gap 
A and a second breather with gap -\/3A. As described in Ref. [15], the 



J ao is the period of the spin chains that for BaCu2Si207 is equal to half the lattice constant in 
c-direction. 
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expectation value of the staggered magnetization can be calculated from the 
results of Ref. [211]: 



2 

23 



mo 



£(1) 
r(i) 



£(i) 
r(i) 



In 



(3.123) 



Equation (3.123) is the self-consistency equation of the mean field approxi- 
mation (recall that h = mo J') and has the solution [342] 



mo 



V2 

7 3 

347r2 



rXI) 
r(i) 



£(i) 

r(i) 



J 



In 



2.58495 J y 



J' 



(3.124) 



We note that the constant 2.58495 should not be taken too seriously as we 
have ignored subleading logarithmic corrections. The result (3.124) is found 
to be in good agreement (for small J' / J) with a phenomcnological expression 
obtained from quantum Monte-Carlo simulations in Ref. [263]. The soliton 
gap as a function of the staggered field h has been calculated in Refs. [15,95]. 
Expressing h in terms of mo by (3.118) and then using (3.124) we obtain 



A 

7 



1 



r(|) 



n 2 



£(i) 

r(i) 



J 



In 



2.58495Jj 



J' 



(3.125) 



Let us now turn to the dynamical susceptibilities. A very useful observation 
is that because /3 = ^ in the sine-Gordon description of the mean-field Hamil- 
tonian (3.118), the correlation functions of the staggered magnetizations are 
related by the SU(2) symmetry (3.109). Hence'''' 

{S^+,{t)Stmc^c\o{-iy{sH\^{t,x)) sin(i$(0,0))) , 
(5jViW5f(0)) ^c%(-l)^(sin(l$(t,x))sin(i$(0,0))> , 



^2 
1 



1 



(5|+i(i)Sf(0))^c2ao(-l)^<cos(-$(t,x))cos(-$(0,0))) . (3.126) 

The dynamical susceptibilities are calculated as in section 3.3.4. However, 

now we are interested in both the real and imaginary parts. Taking all con- 
tributions from intermediate states with at most two particles into account 



^ Recall that due to the shift (3.121) the staggered component of S| is now the cosine rather than 

the sine. 

' The same constant enters both transverse and longitudinal correlations because the spin rotational 
symmetry is restored at high energies. 
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one finds [93] 




27rZi 2|Ff°(0)|^ 
' ac'- 



s2 - [2Acosh(^/2)]2 + ie 

+ / I (3.127) 

s2 - 4A2(1 + ^ cosh e) + ie} 

where s'^ = uj"^ — v^q^. The functions Zi, F|'";(^?) and Fl™{9) are given 
by (3.84), (3.86) and (3.87) respectively, where (3 = ^ and hence ^ = |- 
Similarly one finds 



TT 



xfd(a;,- + g)«.C<^-^— + / d9 




27rZ2 p ,^ 2 1 (0) 1 2 + I Fff (e) 1 2 



ao ^ 1 s2 - 3A2 + ie Jq s"^ - [2A cosh(6'/2)]2 + ie 



s2 - [V12Acosh(e/2)]2 + ie 



(3.128) 



Here Zg, F^°!(6I) and Fff (0) arc given by (3.78), (3.80) and (3.82) respec- 
tively, where /3 = ^ and hence ^ = |- The contribution from B2B2 two- 
particle states is given in [93].™ 

3.6.2. Random Phase Approximation for the Interchain Coupling 

The mean-field approach clearly does not give a good description of the mag- 
netic response in the ordered phase: spin excitations have a gap and there 
are no Goldstone modes. In order to reproduce the latter it is necessary to 
go beyond the mean-field approximation. It is possible to develop a pcrtur- 
bative expansion in the interchain coupling along the lines of [44,45,271]. 
This gives 

" l-2J'(p)[xfd"(a;,,) + S-(a;,p,fc)] ' ^^'^^^^ 

where a = x,y, z. In Eq. (3.129) are the self-energies that are expressed 
in terms of integrals involving three-point, four-point etc correlation func- 
tions of spin operators. The analogous expressions in the disordered phase 
were derived in Refs. [44,151]. To date, the relevant multipoint correlation 
function have not been calculated for the SGM. Furthermore, there is no 
small parameter in the expansion (3.129). The reason is that the mean-field 
gap is generated by the interchain coupling itself and as a result the gap 



' Our normalizations are related to those of [93] hy Z = Att^C 
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(3.125) is of the same order as the interchain couphng. Setting this issue 
aside, a simple Random-Phase Approximation in the interchain couphng 
(3.129) was suggested by Schulz in [271]. Its essence to simply neglect the 
self-energies in (3.129). In other words, one sets 

S'*" = 0. (3.130) 

One problem is that in this approximation the transverse susceptibility will 
not have a zero-frequency spin wave pole at the 3D magnetic zone-center, as 
it should, spin waves being the Goldstone modes of the magnetically ordered 
state. In order for the pole to be exactly at a; = the full self-energy 
must be included. A work-around was suggested in [271]. If the RPA is a 
good approximation, there will almost be a zero frequency pole 

l«.2J'(0,7r)xS(0,^). (3.131) 
The idea is to replace (3.131) by an equality 

l = 2J'(0,7r)xS(0,i^), (3.132) 

and then use (3.132) to fix the overall normalization of Xij- Following this 
logic, we may carry out the integral in (3.127) numerically and obtain 

CRi 0.0645 — . (3.133) 

Now it is a simple matter to determine x^^i^^P' ^ + ?) by evaluating the 
ID susceptibilities numerically and then inserting them into (3.129). 

As discussed in Refs. [93,271], the resulting dynamic susceptibility for 
transverse spin fluctuations ™ coupled chains model contains a 

pair of spin wave excitations that disperse perpendicular to the spin chains, 
and are, by design, gapless. At higher energies there is an incoherent scatter- 
ing continuum with threshold 2A, independent of the transverse momentum. 
This momentum independence is likely to be an artefact of the approxima- 
tion made, rather than being a genuine feature of the model. The RPA result 
for the longitudinal susceptibility exhibits a gapped, coherent mode that 
disperses perpendicular to the chains as well. Above a threshold of 2A an in- 
coherent scattering continuum occurs. Like for the transverse susceptibility 
the threshold for the continuum does not depend on the transverse momen- 
tum. Within the RPA the longitudinal mode is sharp, i.e. it has an infinite 
lifetime. This is certainly a feature of the approximation made: a decay 
of the longitudinal mode into a pair of spinwaves is permitted both by the 
quantum numbers of the excitations involved and phase space. Such decay 
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processes are simply not taken into account in the RPA. The physical pic- 
ture that emerges in the RPA has been found to be in good agreement with 
inelastic neutron scattering experiments on KCuFs [193]. There a damped 
longitudinal mode has been observed. It occurs approximately at the energy 
predicted by the RPA and its spectral weight is found to be close to the 
RPA result. On the other hand, the RPA does rather poorly when applied 
to BaCuaSiaOr [340,342]. 

4. 0(3) Non-linear Sigma Model and Gapped, Integer Spin 
Chains 

The existence of a gap in one-dimensional, integer-spin, Heisenberg antifer- 
romagnets was first predicted by Haldane [137]. He found that such spin 
chains can be mapped onto a gapped integrable field theory, the 0(3) non- 
linear sigma model (NLSM), in the large-spin, continuum limit. It is the goal 
of this section to explore the properties of the 0(3) NLSM and its relevant 
predictions for the physics of gapped spin chains. 

While the 0(3) NLSM is derived in a large S limit, a variety of checks 
imply that this behavior persists down to spin iS = 1. A spin-1 chain with 
a specific {S • S)'^ coupling has been rigorously shown to exhibit a spin 
gap [4]. While at a differing value of the {S ■ S)'^ coupling, the spin chain is 
gapless [23,195,295], this critical point is believed to be unstable in the two- 
parameter space of couplings. Gapless behavior thus only arises as a product 
of fine-tuning. Numerous numerical studies carried out on spin-1 chains 
observe a gap [50, 70, 129, 130, 223, 224, 249, 262, 281, 291, 293, 309, 310, 317]. 
Experimentally, inelastic neutron scattering studies on a number of quasi- 
one-dimensional, spin-1 chain materials are consistent with a finite spin gap 
[49, 217, 218, 225, 231, 257, 283, 303, 324] . 

Here we will suppose that the spin chain is described by a minimal 
Heisenberg Hamiltonian and so ignore (for the most part) the affects of 
anisotropics upon the physics. These can take (at least) two forms. Easy- 
axis anisotropics, 

AH = J2(Sff + DyY,{Stf + Y^^Stf, 

i i i 

of varying strengths are often found in spin-1 chain materials. Furthermore, 
actual spin chain materials never take the form of isolated chains. Rather the 
chains exist in three dimensional arrays with weak but non-zero interchain 
couplings, J'. Thus the chains are at best quasi-one-dimensional. With a 
finite J', there will be some correspondingly finite Neel temperature, Tjv. 
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Below r^v the (ultra low-energy) physics, consisting of higher dimensional 
magnon modes, will be dramatically different than that described here. 

CsNiCls was the first material for which evidence of a Haldanc gap was 
found [49,225,303,324]. CsNiCls has only a weak easy axis anisotropy while 
possessing an interchain coupling, J', of strength J' / J ~ .017. This coupling, 
when combined with the large number (6) of nearest neighbour chains, is 
sufficient to induce Neel order at T ~ 5K. Provided however one is interested 
in physics at energies scales around the gap (A •4J), the absence or 
presence of long range order plays an unimportant role. Another material in 
which a Haldane gap is present is AgVP2SQ [231]. It has an extremely small 
interchain coupling, J' /J ~ 10""^ and a similarly small easy axis anisotropy, 
D^/A ~ 10~^. However it possesses a comparatively large gap, A ~ 320K. 

Other Haldane gap integer spin chain materials arc studied experimen- 
tally, but these compounds are characterized by large easy plane anisotropies. 
One such material is Ni(C2H8N2)2N02C104 (NENP) [217,218,256,257,325]. 
It has an easy-plane anisotropy given by Dz/J ~ .16; Dz/A ~ 2/5 
{D^ Dy 0) [325]. For NENP, the ratio J' /J ~ 8 x 10""^ is suffi- 
ciently small that 3D Neel order has not been observed down to temperatures 
~ 1.2K. Related materials Ni(C5Hi4N2)3(PF6) (NDMAP) [131,338,343,344] 
and Ni(C5Hi4N2)2N3(C104) (NDMAZ) [145, 146, 339] share similar easy- 
plane anisotropies. These latter compounds share the additional feature 
of field-induced antiferromagnetism [60,145-148,341]. The Luttinger liq- 
uid that results from magnetic fields large enough to extinguish the Hal- 
dane gap leads to quasi-long range antiferromagnetic correlations. With 
a small finite J', these quasi-long range correlations are promoted to full 
fledged long range order. The corresponding Neel temperature increases 
with applied magnetic fields. In this review we will not explore this phe- 
nomena nor related behavior occurring in the presence of a magnetic field 
(but see [9, 10, 106, 119, 123, 133, 187, 206, 220, 239, 240, 276, 292, 300]). 

The physics underlying the Haldane gap is particularly robust: related 
systems such as two- leg spin- 1/2 or Hubbard ladders also exhibit a gap to 
spin excitations [67,119,123,187,220]. (We will see this gapped behav- 
ior again when we discuss ladder materials in the context of SO (8) Gross- 
Neveu). Roughly speaking, integer-spin composites form across the rungs of 
the ladder making it into an effective integer-spin chain. Both the ability to 
fabricate these materials and their relationship to high Tc cuprate supercon- 
ductors have made them the focus of intense theoretical and experimental 
study [27, 82, 102, 178, 184, 185, 203, 236-238, 272]. 

In this section we first turn to how to derive the 0(3) NLSM. We will 
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then provide a description of its spectrum, scattering matrices, and form 
factors. We then turn to using these form factors to compute correlators 
both at zero and finite temperature. Our goal at finite temperature will be 
to analyze whether transport in the 0(3) NLSM is ballistic or diffusive. 

4.1. From Integer Spin Heisenberg Chains to the 0(3) 
Non-linear Sigma Model 

2i-4 2i-3 2i-2 2i-l 2i 2i+l 2i+2 2i+3 2i+4 2i+5 

f \ i \ i \ f \ i \ 

• • • • • 

k-2 k-1 k k+1 k+2 

Figure 14. The Heisenberg spin chain 

We begin by deriving the 0(3) NLSM from the Hamiltonian of the Heisen- 
berg spin chain 

i/ = jJ^Si-Si+i, (4.1) 

i 

where S is a spin of arbitrary size. We begin the derivation by decomposing 
each spin, S"" (a=l,2,3), in terms of two operators, n" and L"", via 

S2i = •S"-2j+l/2 + ^2i+l/2 ' (4-2) 

where 

^^21+1/2 — (.^2i+l ~ ^2i)/'^^ ; 

^2m/2 = (^2Vl + 52")/2. (4.3) 

Roughly speaking, the and operators govern antiferromagnetic and 
ferromagnetic fluctuations respectively. For convenience we introduce a new 
lattice indexed by k (see Fig. 14) with sites at 2i + 1/2 and so with twice the 
lattice spacing of the original lattice. Indexed in terms of this lattice, the 



° This development of the 0(3) NLSM follows I. Affleck's in Ref. [5]. 
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operators, n" and L", have the following commutations relations 

[^k^^k'] = ^(^"'^''LlSkk' ; 

[nln\,]=ie'''"^^Ll5kk'. (4.4) 

where 5^^/ = 6^1 / 2, 6^1 being the Kronecker delta function on the original 
lattice, the factor of two taking into account the change in lattice spacing. 
The middle commutation relation indicates that in the continuum limit L 
becomes the generator of rotations for the field n. We can also easily show 
that these operators obey the relations 

n-L = 0; 

n • n = 1 + - - . (4.5) 

s 

In the large s limit we see n ■ n = 1. Our assumption that n satisfies this 
relation in general will be the main approximation of this derivation. 
We are now in position to recast the Hamiltonian. Noting that 

S2i • S2i+i = 2L2i+i/2 • L2i+i/2 + const. ; 

S2i • S2i-i = 2L2i+i/2 • L2i+i/2 - 2s^n5^n + 2s(L • d^rv + ^^n • L) 

+ const. , (4.6) 

we can write the original Hamiltonian as 

H = jj y|4L • L + 2s'^{d^rif + s(L • d^^n + d^n ■ L) 

= ^ 1 jff [L + |^5.n]' + ^ (5,n)2| , (4.7) 

where the spin velocity, f , equals 2Js, the coupling constant g is given by 
2/s and the parameter, 9, the theta angle, equals 27rs. In the above two 
equations, we have set the (original) lattice spacing to 1. We henceforth also 

set V = 1. 

This Hamiltonian corresponds to the Lagrangian 

£ = 2^ • a^n - |-e'*'^n • (^^n x 5,n) , (4.8) 

subject to the constraint that n • n = 1. This model is known as the 0(3) 
NLSM with theta angle, 9. The Lagrangian is given solely in terms of the 



68 F.H.L. Essler and R.M. Konik 



field n. L, the generator of rotations for n, is expressible in terms of n and 
it conjugate momentum, p. p is given by 

1 6 

p = - dtn H n X dxn . (4.9) 

9 47r 

L is then readily found to be 

L = n X p. 

We have that L • n = in accordance with Eq. (4.5). Using this form of L 
we verify that H = p ■ n — C is consistent with Eq. (4.7). 

The theta angle plays a crucial role in determining the properties of the 
above Lagrangian. The term to which it serves as a coupling constant, 
e^'^n ■ {df^n x dun)/8iT, counts (when integrated) the number of times the 
field n winds itself about the sphere (i.e. the second homotopy group of the 
sphere). As it is always an integer, this term in the Lagrangian only affects 
the physics if s is a half integer (and so 6 is an odd multiple of tt). In this 
case this Lagrangian describes a massless theory which at low energies is 
equivalent to the SU{2)i Wess-Zumino-Witten theory. If on the other hand 
s is integer, the theta term has no affect on the physics and the theory is one 
of gapped bosons. It will be this latter case in which we will be interested. 

4.2. Basic Description of the 0(3) Non-linear Sigma Model 

With ^ = 0, the low energy theory of integer spin chains is then simply 

3=^1 dtdx (a^n • d„n). (4.10) 

In terms of and n", the original spins of the theory appear as 

qo. _ a I 

•^21+1 — ■S"2i+l/2 ^ -^2i+l/2 

= sn«(x) + |5^n«(x) + V'ix) + ^5^L"(x); 

^21 = ~'^'^2i+l/2 + -^2i+l/2 

= -sn\x) + U)xn^{x) + L\x) - h)^.L\x). 

The staggered component of the spin (governing fluctuations near /c ~ vr) is 
then sn(a;) + 5a;L(x)/2 ~ sn(x) (keeping into mind we are working at large 
s) while the corresponding smooth component (governing fluctuations near 
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A; ~ 0) on the spin (termed M), is given by 

s 1 

M{x) = L(x) + -d^n = -n{x) x dtn{x), 
^ 9 

We thus see M(x) is quadratic in the field n{x). Unhke its ancestral theory, 
the 0(3) NLSM is integrable [313,332,333]. Local conserved charges were 
constructed in Ref. [127, 254] while non-local conserved charges were first 
found by Liischer [208]. 

The low energy excitations in the 0(3) NLSM take the form of a triplet 
of bosons. The bosons have a relativistic dispersion relation given by 

£(p) = (p2 + A2)V2. 

Here A is the energy gap or mass of the bosons related to the bare coupling, 

via A ~ Je~'^'^l^ . The dispersion relations of all three bosons are identical 
as the model has a global SU(2) symmetry. The exact eigenfunctions of the 
0(3) NLSM Hamiltonian are then multi-particle states made up of mixtures 
of the three bosons. Scattering between the bosons is described by the S- 
matrix [332]: 

<^a3 



e{e - m) 

{9 + m){e - i2TT) ' 



2m{m — 9) -.-.^ 

Here again 9 parameterizes a particle's energy /momentum via E = 
Acosh(0), P = Asinh(0). These S-matrices are relativistically invariant. 
While we stress that this relativistic invariance is a natural feature of the 
low energy structure of the spin chain, we do point out for spin-1 chains, 
A ~ AJ. As J serves as the cutoff for the theory, the low energy sector of 
the theory in this case is not unambiguously defined. 



4.2.1. Form Factors for the Staggered Magnetization Operator, n 

We will be interested in computing correlators involving the field, n, de- 
scribing excitations near A: ~ vr. Form factors for the field n(x, t) have been 
computed by both Smirnov [280] and Balog and Niedermaier [31]. However 
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Ref. [31] presents them in a more amenable form, possible in this particu- 
lar case because of the simple structure of the S-matrix of the 0(3) sigma 
model. 

The staggered magnetization operator is the fundamental field in the 
theory. It thus has non-zero overlap with single particle excitations. By the 
Z2 symmetry, n — n, it also only couples to odd numbers of particles. 
Using the axioms of Section 2, Ref. [31] finds for the one and three particle 
form factors: 

= (0|n„(0,0)|^(^)) = Sab ; (4.12) 

fh:di^i,02,03) = {0\naiO,0)\Ad{e3)AMA,iei)) 

where 9ij = 9i — 9j and ip{9) is defined in Eq.(4.14). 



4.2.2. Form Factors for the Magnetization Operator, M 

We will also be interested in computing correlators involving the magneti- 
zation, M. This field, M = Mq, is part of a Lorentz two-current (Mo,Mi). 
As this current is topological, both components of the current can be written 
in terms of a single Lorentz scalar field, m{x,t): 

M^{x,t) = e^^d''m{x,t) . (4.13) 

We note that only even numbers of particles couple to the magnetization 
current and density operators. This is a consequence of M being express- 
ible as a bilinear in the (fundamental) field n, the field which creates the 
excitations. 

Given Eq.(4.13), it is sufficient to compute the form factors of the scalar 
field m: the form factors of M are then given by 

f^--ar,i&l,- ■ ■ ,(^n) = ^ixi'P''{9i)f^...a„i01, - ■ ■ ,&n), 

where P° = J2i^ cosh(6'i) and = Y.i^ sinh(6lj). Ref. [31] then finds the 
following for the two and four particle form factors of m: 
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(4.14) 







— 6'32 

'3i-27rz 

y^32 — ^31 J 



The reader should note however that we use a different particle normalization 
than Ref. [31] and so the above presented results differ by multiplicative 
constants. 



4.3. Zero Temperature Correlators 

In this section we consider the zero temperature dynamic spectral functions 
of the spin operator near A; ~ and A; ~ tt. These spectral functions govern 
the response in inelastic neutron scattering experiments. 



4.3.1. k ^ Magnetization Correlations 

This problem was first studied theoretically in Ref. [11]. Prom this work, we 
expect at A; ~ to see a two-magnon continuum. A high resolution study 
mapping out this continuum has been carried out in Ref. [327]. Observing the 
two-magnon continuum is a difficult task as the scattering function vanishes 
at small momentum (as we will demonstrate). 
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The spectral function for an integer spin chain near A; ~ is given by 

/•oo 



S{(jO, k) = Im 

TT 



x{-(r(M3(x,r)M3(0,0)))}| 



(4.15) 



Here (r(M^M^)) is an imaginary time-ordered correlator. We evaluate this 
correlator using two particle form-factors. As the next form factor that 
contributes is a four-particle one, our result will be exact for energies, lo < 
4 A. Moreover we can expect the contribution to the overall spectral weight 
of the four-particle form factors and beyond to be small. 

The time-ordered correlator evaluated to the two-particle level equals 

G{x,t) = -@{t){M^{x,t)M^{0,0)) - @{-t){M^{0,0)M^{x,t)) 

3 



il,a2=l 
3 



(4.16) 

where 9 = 9 /{2tt), c{9) = cosh(^), and s{9) = sinh(^). If we now take Fourier 
transforms in space and time, make the necessary analytic continuation in 
u, and use the results for the form factors in the previous section, we find 



S{u; >0,k) = k^ — 



1 1 (anh^((9i2/2) e'f^ + tt- 



712 



cosh 



16 - A;2 ^/^:P^^W^^A^ 9^^ 47r2 + 9^^ 



We have only given the spectral function for w > Q. This response function 
of an incoherent two magnon continuum is plotted in Fig. 15. We can ob- 
serve there the lack of spectral weight at small momentum [11]. Ref. [327] 
compares these theoretical predictions with a neutron scattering study of 
CsNiCla. 



4.3.2. « TT Magnetization Correlations 

We now consider the dynamic response function for wavevectors near vr. 
While at A: = the response is governed by the magnetization operator, M, 
here at A; = tt the relevant operator is n. We expect then that the response 
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Figure 15. A plot of two magnon contribution to the spectral weight of the spin-spin correlator 
around fc = in the 0(3) NLSM. This plot is exact (within the confines of this model) for energies 
less than 4 A. 



will be dominated by a coherent single mode (the magnon). Beyond single 
magnon excitations exists a three magnon continuum. The continuum's 
contribution to the response function has been studied in the context of the 
0(3) NLSM in Ref. [96, 149] where it was found to be small. This accords 
with the intuition built up in Section 2.5 of this review that higher particle 
form-factors make negligible contributions to their corresponding spectral 
function. 

The computation of the dynamic response function in terms of the 0(3) 
NLSM amounts to computing 



1 poo poo 

S{u;,k) = — dx dTe'^^-'''^{-{T{n^{x,T)n^{0,0)))%^-i^+,. 

^ J — oo J — oo 

(4.18) 

We compute contributions to this correlator up to and including three 
magnon form- factors. The next contribution comes at five magnons and 
so this result will be exact for frequencies, w, up to 5A. The time-ordered 
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correlator evaluated to the three-particle level is given by 
G{x,t) = -Q{T){n^{x,T)n^{0,0)) - Q{-T){n^{0,0)n^{x,T)) 
= -e(r)| [ deJ2 |/^^(e)|2e-"^^^W+^"^^^W 

ai=l 

01,02,03 = 1 



-{t,x^ -T,-x). (4.19) 

If we again Fourier transform, using the form-factor results of the previous- 
section we find for positive frequencies [96] 

S{u; >0,Tr + k) = Si{u, k) + ^3(0;, A:) 

+ — / (37r2 + + Y^il;{2z)^l^{z + Y)^{z - Y)Y 
3 7o 

X , ^ , (4.20) 

^(^2 _ ^2 _ ^2 _ 4^2 cosh^ (z) )2 - 16 A4 cosh^ {z) 

where 

zo = cosh~^[(x - l)/2] ; 

^ 1 x2 - l-4cosh2(2;) 

r = cosh \ — ; 

4 cosh(z) 

x^ = i^^. (4.21) 

The (5-function term, Si{u!,k), in the above expression for S{u>,k) corre- 
sponds to the coherent one magnon contribution. The next term, 53(0;, fc), 
arises from the incoherent three magnon continuum. Prom Lorentz invari- 
ance, both are a function of s = Vw'^ + k"^. 

We can evaluate the three magnon contribution numerically. The result 
is plotted in Fig. 16 as a function of s. It is zero for s < 3A. At threshold, i.e. 
s — 3A small and positive, S{lo, k + 11) behaves as (,s — 3A)'^. Wc sec that it 
peaks at around s ~ 6A. Furthermore it is clear that its contribution to the 
response function is small. We can make this observation more qualitative. 
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If we define 



r 

Ii{k + 7r) = / 



'30A 



(4.22) 



^0 



we find that 



^3(vr) 



= 0.018. 



(4.23) 




Thus 98% of the total weight at wavevector tt is found in the coherent one 
magnon contribution [96,149]. 

While the three magnon continuum's contribution to the spectral function 
is small, we do note that there are reports of its observation in the Haldane 
gapped material, CsNiCls [168, 170]. These reports suggest that the three 
magnon continuum is larger than computed in the 0(3) NLSM. Ref. [168] 
finds instead that the ratio in Eq. (4.11) is given by /s (7r)//i (tt) ~ 0.1. That 
the 0(3) NLSM might underestimate the three magnon weight is not an 
impossibility. It is a low energy description of an integer spin chain. Given 
the spin gap, A, equals 0.4J where J, the spin-spin interaction strength, is 
an effective UV-cutoff, we would not necessarily expect accurate predictions 
for energies beyond to > J. On the other hand, the NLSM calculations are 
in fair agreement with numerical computations [192,289], which suggests 
that the large continuum reported experimentally has its origin in physical 
mechanisms not included in the description by means of a spin-1 Heisenberg 
model. The effects of one such candidate, a biquadratic interaction of spins 
[302] , were investigated in [96] . A calculation of the corresponding response 
function gives ^ .2. 

4.4. Low Temperature Properties of Correlation Functions 

In this section we show that finite temperature correlation functions admit 

low temperature expansions as discussed in Section 2. We begin by comput- 
ing the finite DC susceptibility of the system. As we can access this quantity 
exactly through other means, we are able to test the methodology. We will 
find that the form factors reproduce precisely the known exact results. We 
then turn to computing the NMR relaxation rate as well as the finite field 
spin conductance. These quantities both probe the nature of transport in 
integer spin chains as described by the 0(3) NLSM. We will thus be able to 
argue that transport here is ballistic not diffusive. 
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Figure 16. A plot of the three magnon contribution to the spectral weight of the spin-spin cor- 
relator near wavevector tt in the 0(3) NLSM as a function of s. 

4.4.1. Magnetization 

The susceptibility, H = can be computed from the magnetization- 

magnetization operator using a Kubo formula: 

^^H = 0)=Ciu; = 0,k = 0), 



dx / ciTe™""e*''^(r(Mo^(x,r)Mo'(0,0))) 
Jo 



To evaluate this correlator we first expand out the thermal trace: 



. (4.24) 



(M3(x,r)M3(0,0)) = 4Tr(e-'3^O(x,i)O(0,0)) 



Z 



Ens„ e-^""'" Sn\0{x, t)O{0, 0)|n, S„) 

En5„e-^^-" {n,Sn\n,Sn) 



(4.25) 



Here \n,Sn) is a state of n excitations with spins described by Sn = 
{si,-- - ,Sn}- In writing the above we have suppressed sums over the en- 
ergy and momenta of the excitations. A term in the thermal trace with n 
excitations is weighted by a factor of e~"'^^. Thus, as discussed in Section 
2, at low temperatures it is a good approximation to truncate this trace. 
For this computation we keep only terms with one and two excitations, i.e. 
n = 1, 2. To evaluate the matrix elements appearing in Eq. (4.25), we insert 
a resolution of the identity in between the two fields. As we only consider 
matrix elements involving one and two excitations from the thermal trace, 
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we thus have 

(si|M3(x,r)M3(0,0)|si) = ^ {si\M^{x,T)\mSm){mSm\Mi{0,0)\s,) 

mSm 

= (si|Mo^(x,r)|4){4|Mo3(0, 0)|si) + • • • ; 



(siS2|M3(a;,r)M3(0,0)|s2Si) = ^ {siS2\M^{x,T)\mSm) 

mSm 

x(m5^|M3(0,0)|s2Si) 
(siS2|Mo^(x,r)|s'i4)(4s'i|M3(0,0)|s2Si)) + • • • . (4.26) 



In the above we have truncated the sum arising from the resohition of the 
identity. With the first matrix element of the thermal trace, we only keep 
terms from the resolution of identity with one excitation. We are interested 
in the behavior of the susceptibility at a; = and this term provides the 
only contribution [188]. Similarly, the only term arising from the second 
matrix element of the thermal trace contributing to the DC susceptibility 
comes from keeping the term from the resolution of the identity involving 
two excitations. We can then evaluate C{uj,k) with the result 

C{co,k) = Ci{L0,k) + C2{u},k). (4.27) 

For Ci{u!,k), we then have the corresponding expression 

C,{x,r) = Y: J ^^e-^/^-M^)(^.,(^)|M3(x,r)|^,(ei)> 

x(^,(^i)|m3(0,0)|^,(^)) 

= E 

= E 



^ ^ g-A^cosh(6l)g-TA(cosh(ei)-cosh(6»))+ia;A(sinh(6»i)-sinh(6l)) 
27r 27r 



^ ^ g-A/3cosh(6l)g-rA(cosh(6>i)-cosh(6/))+ia:A(sinh(ei)-sinh(6l)) 

27r 27r 

SlS2 " 

xf^iie ~ i7r,ei)f^A{ei - i7T,e). (4.28) 

We have used crossing symmetry in the second hne. From Section 4.2.2, the 
form factor faai ^i) is given by 

2 A 

f^a!ie,9i) = i^e^'^Asmh{0) + smhiei))4^ie-ei). (4.29) 



78 F.H.L. Essler and R.M. Konik 



Then upon Fourier transforming the above in x and r and continuing a;„ — ^ 
—ico + S, we obtain 

Ci(a; = 0,k = 0) = ^ I decosh{e)e-'^^'^'^^^^ = ^i^i(/3A), (4.30) 

TT J TT 

where Ki is a modified Bessel function. This has the expected small tem- 
perature behavior, Ci(c<; = 0, fc = 0) ~ T^^/^e"^^. 
On the other hand €2(00 = 0,k = 0) is given by 

d9i d92 dO^ dOi 



27r 27r 27r 27r 

SlS2Sj,S4, •■ 



x(^,(ei)^2(^2)|M3(x,r)|^3(f?3)^4(^4)) 

x(^,(e4)^3(^3)|M3(0, 0)1^,(^2)^1(^1)) 

-Ci(x,r) j ^e-^^-^^W(^„(^)|^„(0)). (4.31) 

The last term is disconnected (proportional to Ci(x,r)) and will ultimately 
cancel (good, as it is proportional to (5(0)). It arises from the expansion 
of the partition function in Eqn. 4.25. The four particle matrix elements 
appearing in the above takes the form (as per Section 2.7.1): 

{As,{eMsM\Ml{x,T)\As,{eMsM)) 

+ S,,^,, 2775(03 - ^2)<^.1(^12)4i(^34)/j^yei -i7r,^4) 
+ <5,, ,,27r,5(^2 - ^4)4.1 (^i2)/i^° 3(^1 - iTT, 63) 

+ 5,^,>27r6{ei - ^3)4.1 (^34)/j^y^2 - i7r,e^) 

+fSus,,ss{02 - i7T,ei - m,^4,^3)c , (4.32) 

where fc refers to a connected form-factor. Substituting this expression into 
the above and using the form factors of Section 4.2.2, we find (see Ref. [188] 
for details) 

C2{uj = 0,k = 0) = i^i(2/3A) 

TT 

+ / (ieid^ne-^^(^°^'^(^i)+'=°'^^('^2))coshf^i) llTT^ + 2gf2 

= Ki(2/3A) + ^Ko(/?A)ifi(/3A) + o(^e-2^^) , (4.33) 
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where O12 = 9i — 62 and Kn are standard modified Bessel functions. The 
first term in C2 is a 'disconnected' contribution related to Ci. The second 
term is a connected contribution and as such is genuinely distinct from Ci. 
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Figure 17. Plots of the zero-field susceptibility computed both from the TBA equations and from 
the form factor expansion. The first of these is an exact numerical solution of the TBA equations 
for the 0(3) non-linear sigma model. The second is arrived at from a small temperature expansion 
in powers of e~l^^ of these same equations. The final plot gives the form factor computation of 
the susceptibility. We have truncated the form factor expansion at the four particle level. 




It is possible in the case of the 0(3) sigma model to arrive at exact 
expressions (in the form of coupled integral equations) for the zero-field sus- 
ceptibility [301,313]. These equations, in their most compact form, appear 
as 



^(^ = °) = -^/^^^°^^(^)-l + e/3^W ' 

€(9) = A cosh{e) -T j dO' log(l + e'^^2{e')^ 



e^{e) = T J de's{e - e')yo^{i + e^^-^^^')) 

+ log(l + e'^^'^+i^^')) + 62n log(l + e^<''^] 
H = hm ^ . (4.34) 

n^oo n 

These equations admit a closed form low temperature expansion. The details 
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of this expansion may be found in Ref. [301]. Here we just quote the results 



TT TT 



Remarkably, we see this expansion agrees exactly with the corresponding ex- 
pression derived with the aid of form factors. Thus the form factor expansion 
at finite temperature meets an important test. 

In Fig. 17 are plotted the susceptibilities computed via an exact numeri- 
cal analysis of the TEA equations, a low temperature expansion of the same 
equations, and a computation based upon the two and four particle form 
factors. We see that as indicated previously that the form factor computa- 
tion and the low temperature expansion match exactly. Moreover these two 
computations track the exact susceptibility over a considerable range of tem- 
peratures despite the fact these computations are truncated low temperature 
expansions. 

4.4.2. NMR Relaxation Rate 

In this section we compute the NMR relaxation rate, l/Ti. We are interested 
in computing this rate in order to compare it to the experimental data found 
in Ref. [285] on the relaxation rate of the quasi one-dimensional spin chain, 
AgVP2SQ. For temperatures in excess of lOOK (the gap. A, in this compound 
is on the order of 320i^), the experimental data [285] shows the relaxation 
rate to have an inverse dependence upon i/ff: 

1/Tioc ^ 



This dependence is nicely reproduced by the semi-classical methodology in 
Refs. [69,259]. Moreover the semi-classical computation reproduces the ac- 
tivated behavior of 1/Ti in this same temperature regime: 

We are interested in determining whether a calculation in the fully quantum 
0(3) NLSM can reproduce these results. To this end we compute l/Ti 
using a form factor expansion. Sagi and AfHeck [261] have already done 
such a computation to lowest order in e~^^. But they do not find the above 
behavior. Rather they see 

l/Ti oc \og{H)- l/Ti oc e-^^. 
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We continue this computation one further step, computing to 0(e~^^^). 
Given the behavior, 1/Ti ~ H^^^"^, appears only as T is increased beyond 
lOOK (i.e. T/A 1/3), it is not unreasonable to suppose higher order terms 
in a low temperature expansion of 1/Ti are needed to see this singularity. 




To proceed with the computation of 1/Ti, we review its constituent el- 
ements. 1/Ti can be expressed in terms of the spin-spin correlation func- 
tion [261]: 

V^i= E /^^a/3(A:)^a7(-fc)(Mo^Mo^)(fc,a;iv), (4.36) 

a = l,2 
/3=1,2,3 

where cun = ^nH is the nuclear Lamour frequency with 77V the nuclear 
gyromagnetic ratio and the A^fj are the hyperfine coupling constants. In 
the above we assume H is aligned in the 3-direction. The above integral is 
dominated by values of k near [261]. Moreover in the relevant experiment, 
the hyperfine couplings are such that only (MqMq) contributes. Hence 

l/Ti oc {mImI){x = 0, wjv ~ 0) . (4.37) 

We now proceed to compute {MqMq). 

To compute {M^M^), we again employ a form factor expansion. Akin to 
the computation of the susceptibility and the spin conductance, this compu- 
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tation amounts to a low temperature expansion of (MJMq), 

where we are able to compute ai and 02- For the details of this computation 
we refer the reader to Ref. [188]. We find there 

{MM)i^ = 0,0; = 0) = e-/^^(log(f ) - 7) - ^ e-^/^^(log(|^) - 7) 
+Ae-^/^^(log(^) - 7)y||(247r + ^J)) (l + 0{H/T) + 0{T/A)) , 

(4.38) 

where 7 = .577 ... is Euler's constant. We are interested in the regime 
^ T <C A (the regime where it is expected spin diffusion produces sin- 
gular behavior in l/Ti). The terms that we have dropped do not affect this 
behavior. In principle there is no difficulty in writing down the exact ex- 
pression (to C'(e^^'^^)); it is merely unwieldy. This expression for 1/Ti is 
plotted in Fig. 18 for a variety of values of the ratio T/A. 

We see that we do not obtain the same behavior as found in Refs. [69, 
259]. Going to the next order in 0{e~'^^'^) produces a behavior in 1/Ti as 
H ^ identical to the lower order computation of 0{e^^^): we again find 
a logarithmic behavior consistent with ballistic transport. An alternative 
comparison we might make to the results of Refs. [69, 259] is to perform a 
low temperature expansion (in 0{e~^^)) of the semi-classical computation 
of {M^M^){x = Q,u = 0). Doing so by treating Te'l^^/H as a small 
parameter, we find 

{mMKx = 0,0; = 0) oc Ae-/^^(log(f ) - 7 + (J - i) -g, e-'^^ 

+0(6-^^^)) . (4.39) 

We see that the low temperature expansion of the semi-classical result agrees 
to leading order with our computation but afterward differs. It possesses no 
term of 0{e'~'^^^). The next term rather appears at 0{e~^^^) and possesses 
a 1 /H'^ divergence. That the small H behavior should be 1/ ^/H does suggest 
the importance of summing up terms. But the lack of a term of 0{e^^^^) in 
the semi-classical result nonetheless hints that the two results are genuinely 
different. 
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4.4.3. Drude Weight of Spin Conductance at Finite Field 

In this section we compute the spin conductivity, Ug- The spin conductivity 
gives the response of the spin chain to a spatiahy varying magnetic field. It 
is defined via 

ji{x,t)=asVH, (4.40) 
and so can be expressed in terms of a Kubo formula, 

Reas{k,u) = -\j da;dte^^-+-*Im(jo(x,t)ii(0,0)),,tarded ■ (4-41) 

In the notation used here the spin current ji is synonymous with Mf of Sec- 
tion 4.2.2. (for a field in the 3-direction), the Lorentz current counterpart of 
the uniform magnetization, Mq = jq. We wih focus primarily on computing 
the Drude weight, D, of Recr^, i.e. computing the term in Us{k,(jo) of the 
form 

as{k = {),io) = D5{u). (4.42) 

However we are able to compute Gs for general k, ui. We find that for uj <C 2A, 
A; = 0, the spin conductivity is described solely by the Drude weight. In 
particular, we find no indication of a regular contribution to as{k = 0,u;). 

To evaluate o"s, we employ the identical form factor expansion to that 
used in computing the susceptibility. And like the susceptibility, our result 
is an exact low temperature expansion of D, 

n 

Here we will compute Di and D2 exactly. As the details of the computation 
are nearly identical to that of the susceptibility, we merely write down the 
results: 

D(H = 0) = pA [ dee-P^'-''^('^^-^^^^(l - Se-Z^^-^i^C^)) 
J cosh(0) 

+23 A [ de^de,c-^^^''''^^^^^^+^°^^^^^^^ sinh^(g2) Utt^ + 291^ 
J M,M2e ^^^^^^^^ ^ ^^2^2^ ^ 

+0(e-^^^) 
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This expression, like the susceptibihty, involves only the two and four particle 
form factors. We plot this result in Fig. 4.4.3 as a function of T/A. 
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Figure 19. In this plot we present the form factor computation of the zero field {H = 0) Drude 
weight, D, of the spin conductance. 

We first observe that D{H = 0) 7^ 0. This is in accordance with Ref. [117] 
where D is computed using an argument involving the finite size scaling of the 
thermodynamic Bethe ansatz equations. (We do note that the computation 
of D at = in Ref. [117] appears only as a note added in proof and 
so is decidedly sketchy. However the equations governing D developed in 
Ref. [117] are manifestly positive with the consequence D cannot vanish.) 
But our results do differ from the semi-classical computation of Ref. [260] 
where it was found that D vanishes at H = 0. We find as well no additional 
regular contributions to as{uj,k = 0) near uj = ~ only the Drude term 
is present in contrast to Refs. [69,259]. (There will, however, be regular 
contributions at higher frequencies, in particular for uj > 2A, which persist 
even in the zero temperature limit). 

We have only given the spin conductivity at H = 0. However it is ex- 
tremely straightforward to generalize the form factor computation to finite 
H. As H couples to the total spin, a conserved quantity, the form factors, 
f^{x,t), of an operator, 0{x,t), carrying spin s, are altered via the rule 

/^(t)^e^^^V^(t). 

(In the case of the spin conductance, the spin currents, = M^, carry no 
spin and so are not altered at all.) The only remaining change induced by a 
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finite field is to the Boltzmann factor appearing in the thermal trace. If an 
excitation with rapidity, 9, carries spin s, its Boltzmann factor becomes 

g-^(Acosh(e)-sif) 
For example we find D clS Si function of H (to ©(e"^^)) to be 

DiH) = ^Acosh(^i7) / d^e-^^-^^W ^^^^ . (4.44) 
J cosh(t^) 

Again this in agreement with Ref. [117]. Indeed Ref. [117] computes D(H) 
at large H/T (but ii" < A) to be 

D = ^^PH (' ^m^M e-/3Acosh(e) + c»(e-2^^) . (4.45) 
47r J cosh(6') 

Up to a factor of 27r, this expression is in exact agreement with 4.44. In 
this particular case our derivation of D{H) agrees with the semi-classical 
computation [260] (provided T <C -ff <C A). The symmetries in the semi- 
classical model that lead D{H = 0) to vanish are broken for finite H. 



4.5. Transport: Ballistic or Diffusive 

We have compared our transport calculations to the semi-classical compu- 
tations in Refs. [69,259]. The essence of this method lies in treating the 
spin-chain as a Maxwell-Boltzmann gas of spins which interact with one 
another through the low energy limit of the scattering of the 0(3) NLSM, 

S^t{e = 0) = -SadScb- (4.46) 

In contradistinction to the semi-classical computation, we have found that 
the Drudc weight of the spin conductance is finite in the limit of zero external 
field. Our results for the NMR relaxation rate, 1/Ti, indicate a similar 
discrepancy. We, like Ref. [261], find that 1/Ti is characterized by ballistic 
logarithms. These logarithms are relatively robust: they continue to appear 
at higher orders in the low temperature expansion. We do not, however, 
see diffusive behavior in the relaxation rate, i.e. 1/Ti ~ nor does 

our low temperature expansion match the low temperature expansion of the 
semi-classical computation of the correlator. 

To come to some sort of judgment between the form-factor and the semi- 
classical approaches, an understanding is needed of the differences between 
our computations of the spin conductance and the NMR relaxation rate. In 
the case of the first quantity, it is likely this difference is real and not an 
artefact of our methodology. The data that goes into the spin conductance 
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is identical to that needed to compute the susceptibihty and we know that 
we can match the low temperature expansion of the susceptibility with a 
similar expansion coming from the exact free energy. Moreover we know 
that the Drude weight of crs{H = 0) has been found to be finite from an 
approach [117] independent of ours. 

In generic systems, the Drude weight, D, at finite temperatures will be 
zero. It is then the integrability of the 0(3) NLSM and the attendant ex- 
istence of an infinite number of conserved quantities that leads to a finite 
weight. The existence of these quantities can be directly related to a finite 
D. As discussed in Refs. [56,345,346], D is bounded from below via an 
inequality developed by Mazur: 

D>cY.^^, (4.47) 



iQl) ' 



where J is the relevant current operator, Q„ are a set of orthogonal conserved 
quantities, i.e. {QnQm) = ^nmiQu), and c is some constant. For a finite 
Drude weight, we then require that at least one matrix element, {JQn), does 
not vanish. While we do no direct computations, we can obtain an indication 
of whether the matrix elements vanish by examining the symmetries of the 
model. Under the discrete {Z2) symmetries of the 0(3) NLSM, the spin 
current, J, transforms via 

Z2iJ)^±J. 

In order that the matrix element, (JQn), not vanish we require that 

The Z2 symmetries in the 0(3) NLSM include Ua —>■ —Ua, a = 1,2, 3, parity, 
and time reversal. The spin current we are interested in transforms under 
rotations as a pseudo- vector. Thus any charge, Qn, coupling to the current 
must also transform as such. Prom the work by Liischcr [208], there is 
at least one conserved pseudo- vectorial quantity such that (JQn) does not 
vanish due to the action of one of the above Z2 symmetries. For the sake of 
completeness we exhibit it. Rewriting the magnetization and spin current, 
Mo,i, explicitly as antisymmetric tensors, 

Mf = n'^d^n^ - n^d^n'' , /i = 0, 1 ; 

the conserved quantity takes the form, 

Q''\t) = Y, J dxidx2 sgn{xi - X2)MS''it,xi)M§\t,X2) - J dxMf{t,x). 
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While the first term of Q does not contribute to the matrix element, 
(JQ"^), as it is bihnear in the currents, M"^, the second term does. We 
point out that Q"'^ is an exotic object inasmuch as it is a non-local conserved 
quantity. As pointed out in Ref. [208], it is the first in a series of non-local 
charges. 

While the structure of the conserved quantities in the 0(3) NLSM seem 
to be consistent with the existence of a finite Drude weight, this is not the 
case in the semi-classical approach. The dynamics of the semi-classical ap- 
proximation used in Refs. [69,259] are also governed by an infinite number 
of conserved quantities. Importantly however, these are different than those 
appearing in the fully quantum model. In particular, the semi-classical ap- 
proximation does not admit non-local conserved quantities. As shown in 
Refs. [69,259], the structure of the Z2 symmetries in the semi-classical ap- 
proach is such that all matrix elements, {JQn), vanish. It would thus seem 
the absence of a Drude weight in the semi-classical case is a consequence of 
differences in the symmetries between the semi-classical and fully quantum 
models. 

To understand the discrepancies in the case of the NMR relaxation rate, 
1/Ti, is not as simple. However if we believe that the spin conductance 
demonstrates finite temperature ballistic behavior, it is hardly surprising to 
find the NMR relaxation rate characterized by ballistic logarithms. Again 
the difference between the fully quantum treatment and the semi-classical 
approach will lie in the dift'ercnces between the models' conserved quantities. 
Nonetheless one possibility that we must consider is that merely going to 
0{e~^^^) in the computation of 1/Ti is insufficient. It is possible that we 
need to perform some resummation of contributions from all orders to see 
the desired singular behavior, 1/Ti ~ I/^/H. While this would belie our 
experience with computing the susceptibility and the spin conductance via 
the correlators, the data that goes into the two computations is not exactly 
identical. Thus the possibility that the low temperature expansion of l/Ti 
is not well controlled cannot be entirely ruled out. 

The differences in the nature of the conserved quantities between the 
0(3) NLSM and the semi-classical model of Refs. [69,259] suggest the latter 
is not equivalent to the 0(3) NLSM, even at low energies. An indication 
of this lack of equivalency may lie in the universal nature of the ultra low 
energy S-matrix. This quantity is the primary input of the semi-classical 
model. The semi-classical model imagines a set of classical spins interacting 
via 



5S:(^ = o) = -<5aA6, 
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i.e. in the scattering of two spins, the spins exchange their quantum num- 
bers. However this specification may be insufficient to adequately describe 
the 0(3) NLSM. Even beyond the quantum interference effects which are 
neglected by the semi-classical treatment, it is not clear that the zero- 
momentum S-matrix is enough to determine the model. 

In this light it is instructive to consider the sine-Gordon model in its 
repulsive regime. The sine-Gordon model is given by the action, 

S = -l- j dxdt {d^<i>di^^ + A cos(/5$)) , (4.48) 

where f3 falls in the range < /3 < 1. The model is generically gapped. Its 
repulsive regime occurs in the range, 1/2 < < 1. The model's spectrum 
in this same range then consists solely of a doublet of solitons carrying U(l) 
charge. It is repulsive in the sense that the solitons have no bound states. 
The sine-Gordon model has a similar low energy S-matrix to the 0(3) NLSM, 

-5^6(^ = 0) = 

where here the particle indices range over it, the two solitons in the theory. 
Thus we might expect that sine-Gordon model to possess identical low energy 
behavior over its entire repulsive regime. 

This is likely to be in general untrue. For example we might consider 
the behavior of the single particle spectral function. We might thus want to 
compute a correlator of the form 

(V'+(x,t)V'-(o,o)), 

where ip± are Mandelstam fermions given by 
il^±{x,t) = exp ^(-^ + V2/3)0i(x,t) T |(^ - V2(3)Mx,t)^ ; 

<t>L/R = \ {^{x, t)±ij dydtHy, t)^ . (4.49) 

Now consider the finite temperature spectral function corresponding to this 
correlator with energy, uj, fixed at just above the one-soliton gap. This 
spectral function has a contribution from a one particle form factor. This 
contribution's /3-dependence, while present, is trivial in that it only appears 
in the overall normalization. But because we are working at finite temper- 
ature, higher particle form factors will also contribute. Their contribution 
will depend upon /3 in a more complex fashion than an overall normalization. 
We thus expect the low energy properties of this spectral function to have 
an overall non-trivial dependence upon (5. This accords with intuition. /3 
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determines the compactification radius of the boson in the model and so is 
related in a fundamental way to the model's properties. 

It is useful to point out that Mandelstam fermions are the unique fields 
that create/destroy solitons that carry Lorentz spin 1/2, i.e. a spin that 
is independent of /3. They would then be the only fields with a chance of 
matching any semi-classical computation. However there are other soliton 
creation fields, for example, 

g±i0i,fl/(V2/3)^ 

for which one could determine the corresponding spectral density. As these 
fields carry spin that varies as a function of /3, their spectral functions will 
depend upon more than the ultra low energy soliton S-matrix. In general, 
the semi-classical treatment of the sine-Gordon model cannot capture its full 
quantum field content. 

As with the 0(3) NLSM, the conductance of the fully quantum model 
differs from that of the semi-classical treatment. If one were to compute 
the conductance at finite temperature in the sine-Gordon model one would 
again find a finite Drude weight, D, while the semi-classical approach yields 
D = [120]. The notion of under-specificity appears here again. The semi- 
classical approach for the sine-Gordon model equally well describes the Hub- 
bard model at half-filling (the solitons are replaced by particle/hole excita- 
tions in the half-filled band). But it fails to give the correct Drude weight. 
An analysis of finite size corrections to the free energy in the presence of an 
Aharonov-Bohm flux [116] again flnds a finite Drude weight in the half-filled 
Hubbard model at finite temperature. 

Interestingly however, there are certain properties at low energies that 
seem to be independent of /3. For example, if one were to compute the 
low temperature static charge susceptibility, the term of 0{e^^^) would 
be independent of p. However at the next order, 0(e~^^^), this would 
almost certainly cease to be true. And the energy/temperature ranges we 
are interested in exploring do not permit dropping terms of 0(e~^^^). 

It is important to stress we do not question the agreement between the 
semi-classical model and experiment. What we do question is whether the 
fully quantum 0(3) NLSM exhibits spin diffusivity. If we are then to un- 
derstand spin diffusion in terms of the 0(3) NLSM, it is possible we need 
to include additional physics explicitly such as an easy axis spin anisotropy 
(weakly present in the experimental system, AgVP2SQ), inter-chain cou- 
plings, or a spin-phonon coupling (as done in [117]). 

Beyond these, another mechanism that might lead to diffusive behavior 
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are small integrable breaking perturbations of the 0(3) NLSM. Generically 
any physical realization of a spin chain will possess such perturbations, even 

if arbitrarily small. Such perturbations may introduce the necessary ergodic- 
ity into the system, ergodicity that is absent in the integrable model because 
of the presence of non-trivial conserved charges, and so lead to diffusive be- 
havior. As discussed in the semi-classical context by Garst and Rosch [120], 
such perturbations introduce an additional time scale, r, governing the de- 
cay of conserved quantities in the problem. For times, t < t, the behavior 
of the system is ballistic and the original conserved quantities do not decay. 
For times, t > t, the behavior is then diffusive. Consequently the Drude 
weight in the purely integrable model is transformed into a peak in a{uj) at 
Lo ~ l/r. 

Now the difference in the physics between the 0(3) NLSM and its 
semi-classical variant is not that of integrable breaking perturbations. As 
demonstrated in Refs. [69,259], their semi-classical model is classically in- 
tegrable. However as discussed above the models do possess different con- 
served charges. It might then seem for certain transport quantities, the 
semi-classical model cures the lack of ergodicity present in its quantum coun- 
terpart. 

5. U(l) Thirring Model and Quasi One Dimensional Mott 
Insulators 

The Mott metal-insulator transition [121, 227, 228] is a zero temperature 
( "quantum" ) phase transition between a gapless metallic phase and a gapped 
insulating one. It occurs at some critical ratio of the strength of the electron- 
electron interaction to the kinetic energy, which is usually measured in 
terms of the tunneling amplitude t of an appropriate tight-binding model. 
In practice U/t can be varied by applying pressure to a crystal, which results 
in a better overlap between the orbitals of the conduction band and leads to 
an increase int. In high-energy physics the Mott transition corresponds to 
the phenomenon of dynamical mass generation. 

The intrinsic difficulty in describing the Mott transition quantitatively 
in the general case is that it occurs when the kinetic and potential ener- 
gies are of the same order of magnitude. This regime is difficult to access 
from either the "band" limit, where one diagonalizes the kinetic energy first 
and then takes electron-electron interactions into account perturbatively, or 
the "atomic" limit, in which the electron-electron interaction is diagonalized 
first and the electron hopping is taken into account perturbatively. Inter- 
estingly, the Mott transition on lattices where tunneling along one direction 
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is much larger than along all others can be understood in some detail by 
employing methods of integrable quantum field theory. The systems under 
consideration can be thought of in terms of weakly coupled chains of elec- 
trons and will be referred to as quasi one dimensional Mott insulators. When 
the band is half-filled and the interchain tunneling is switched off, Umklapp 
processes dynamically generate a spectral gap M and we are dealing with 
an ensemble of uncoupled ID Mott insulating chains. The same Umklapp 
scattering mechanism can generate gaps at any commensurate filling e.g. 
quarter filling, but only if the interactions are sufficiently strong. There are 
two questions we want to address: 

(i) What is the dynamical response of the uncoupled Mott-insulating 
chains system? 

(ii) What are the effects of a weak interchain tunneling? 

As far as point (i) is concerned we will concentrate on the optical conductivity 
and the single particle spectral function. Results are also available for the 
density-density response function [64] and the dynamical structure factor 
[41]. Examples of materials that are believed to fall into the general category 
of quasi- ID Mott insulators are the Bechgaard salts [48] and chain cuprates 
like SrCu02, Sr2Cu03 ° or PrBa2Cu307. They have been found to exhibit 
very rich and unusual physical properties such as spin-charge separation 
[118,173,180,222,326,328]. 



5.1. Lattice Models of Correlated Electrons 

The simplest models used in the description of one dimensional Mott insu- 
lators are "extended" Hubbard models of the form 

H = -tY^ [4,aCn+i,a + h.c. nk,fnk,i + rikUk+j , (5.1) 

n,a k j>l k 

where nj-^^ = c|. ^Cfe,a- and = rifc^-f -|- n^^ j are electron number opera- 
tors. The electron-electron interaction terms mimic the effects of a screened 
Coulomb interaction. Two cases are of particular interest from the point of 
view of application to, for example, the Bechgaard salts [48]: 

(1) half filling (one electron per site); 

(2) quarter filling (one electron per two sites). 



° More precisely, these compounds are considered to be charge-transfer insulators. 
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We will discuss both these cases and emphasize similarities and differences 
in their dynamical response. 

5.2. Field Theory Description of the Low Energy Limit 

The field theory limit is constructed by keeping only the low-energy modes 
in the vicinity of the Fermi points ztkp. We may express the lattice electron 
annihilation operator in terms of the slowly varying (on the scale of the 
lattice spacing ao) right and left moving electron fields R{x) and L{x), 

ci^a — ^ [exp(iA;i?x) Ra{x) + exp(— ifc^a;) i^o-(x)] . (5-2) 

Here kp = 7r/2ao for the half-filled band and x = la^. The resulting fermion 
Hamiltonian can then be bosonized by standard methods [128]. In order to 
make our presentation reasonably self-contained we review the relevant steps 
in the half-filled case next. 

5.2.1. Half Filled Band 

Using (5.2) in (5.1) we arrive at the following low-energy effective theory 
"H = ^ vpjdx: idxL^ — i?^ idxRa ■ J dx : [J • J -|- J • J] : 

-gsjdx :3.3:+^thA^ J :[l ■ I + 1 .l\ : 

+ Jdx^ -.{I+I- +I-I+):+gl -.PP: . (5.3) 

Here vp = 2tao is the Fermi velocity, 

go = 2Uao, Qs = g^ = 2aoiU-2Vi+2V2), ffjJ = 2aoiU +6V1+2V2) , (5.4) 
and J, I are the chiral components of SU(2) spin and pseudospin currents, 

P = \: + : , /+ = (I-)t = , 

P = ^: + : , /+ = (/-)t = r\r\ 

P = i : (lJLt - lJl^) : , J+ = (J-)t = , 

P = ^:(^r\r^-R\Ri): , J+ = (J-)t = . (5.5) 

Here ":" denotes normal ordering of point-split expressions [3]. The "kinetic" 
terms in the Hamiltonian (5.3) can be expressed as normal ordered bilinears 
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of currents as well [3,65, 128] 



27r 



J dx : [: 



I I + J J] : 
i I + J Jl : 



dx 



dx ^ : -RJ- idx Ra ■ 

. cr 

^ : lI idx La : 



(5.6) 



Using (5.6) the Hamiltonian (5.3) can now be split into two parts, corre- 
sponding to the spin and charge sectors respectively 

Ti- = 'Hc + 'Hs , 

U,= '^ j dx : [I-I + I-I] : 



3 

+ j dx 

2ttVs 



i^:[I+r+ri+]:+gl:P-P: 



J dx : [j • J + J • j] : —gs J : J • J : . 



(5.7) 



Heref^ = VF — Uao/2'K andt;^ = VF+{U+4Vi+4V2)ao/2Tr. The Hamiltonian 
(5.7) with Vs = v'f, and gs = is known as the U(l) Thirring model in the 
literature. As we will see the difference in velocities can be accommodated 
and as will show now the current-current interaction in the spin sector is 
marginally irrelevant and will only lead to logarithmic corrections at low 
energies. The 1-loop renormalization group equation for gs is 



dgs 
dr 



9s 



(5.8) 



Hence gg diminishes under renormalization and the current-current inter- 
action in the spin sector is marginally irrelevant. In order to keep things 
simple we will drop it from now on. Keeping it would lead to logarithmic 
corrections in many of the formulas below. The Hamiltonian (5.7) can now 
be bosonized using 



4(3 



4(a 



J]£_JU7r/4 



2n 



exp I -2*?^^ ) exp 



-(ps , 



27r 



exp ( -(pc I exp I 



(5.9) 
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where rja are Klein factors with {r/a, 77^} = 2Sab and where = 1, /j^ = — 1 . p 
The canonical Bose fields ^s,c and their respective dual fields &s,c are given 

by 

^a = (l>a + 4>a , &a = (l^a - 4>a , a = S, C , (5.10) 

where the chiral boson fields (f)a and 4>a fulfill the following commutation 
relations 

[(/)a(x),(^a(y)] = 27rz , a = c,s. (5.11) 

We choose a normalization such that for |x — 2/| — > the following operator 
product expansion holds (a = s, c) 

exp (ia$„(x)) exp (ifi^aiy)) — ' \x - exp {ia^a{x)+iP^a{y)) ■ (5.12) 

Applying the bosonization identities we obtain the following bosonic form of 
the low energy effective Hamiltonian (we recall that we have set Qs = 0) 

~ (i)2 / "^"'^^^^ + I ~ idxQcf] , 

'^s = ^jdx [{d^^s? + {d^Qs?] . (5.13) 
Finally we carry out a canonical transformation on the charge boson 

$c P^c , Qc ^Qc , (5.14) 

where 

r i-gl/47r< 

The transformation property of the dual field follows from the fact that 
—dx@c is the momentum conjugate to $c- In terms of the rescaled fields the 
charge sector takes the form of a sine-Gordon model 

'^c = ^jdx [{d^Qcf + [dx'^cf] j dx cos^$e , (5.16) 



4 

(5.15) 



P The phase factors in (5.17) have been introduced in order to ensure the standard bosonization 
formulas for the staggered magnetizations. 
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where Vc = {v'^ + gc/^Tr)P'^. The entire procedure we have been foUowing can 
be summarized as follows. We first project the lattice Hamiltonian to the 
low-energy degrees of freedom using (5.2) and then bosonize the resulting 
fermion Hamiltonian by means of the identities 



rUx) = 



V<7 



1 



Or 



exp 



exp 



4 



+ e,] . (5.17) 



The Hamiltonian (5.16) exhibits spin-charge separation: Hc,s describe 
collective charge and spin degrees of freedom respectively, which are inde- 
pendent of one another. The pure Hubbard model corresponds to the limit 
f3 ^ 1 and the effect of Vj is to decrease the value of f3. From the form 
(5.16) we can deduce a number of important properties. Firstly, the spin 
sector is gapless and is described by a free bosonic theory. Hence correlation 
functions involving (vertex operators of) the spin boson $5 and its dual field 
@s can be calculated by standard methods [128]. Excitations in the spin sec- 
tor are scattering states of gapless, chargeless spin ^ objects called spinons. 
The charge sector of (5.16) is a SGM. Here excitations in the regime P > 
are scattering states of gapped charge =be soliton and antisoliton excitations 
respectively. In the context of the half-filled Mott insulator these are also 
known as the holon and antiholon. The soliton and antisoliton have massive 
relativistic dispersions with velocity Vc and single-particle gap A, 



E{P) = + v^P^ . 



(5.18) 



For Vj = 0, j > S the gap can be determined by renormalization group 
methods , , , 

- (5.19) 



A 



where we have fixed the constant factor by comparing to the exact result for 
the Hubbard model [202], and where 



X = 



1 - 



U-2V1 + 2V2 
U + 6V1 + 2V2 



1/2 



g = (U + 6Vl + 2V2)/2TTt . 



(5.20) 



We note that the gap vanishes on the critical surface U — 2V\ + 2V2 = 
separating the Mott-insulating phase with gapless spin excitations from 
another phase with a spin gap. 
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In the regime < /3 < 1 /\/2, soliton and antisoliton attract and can form 
bound states. In the SGM these are known as "breathers" and correspond 
to excitons in the underlying extended Hubbard lattice model. There are 

■1-/32- 



N 



(5.21) 



different types of excitons, where [x] in (5.21) denotes the integer part of x. 
The exciton gaps are given by 

A„ = 2A sin(n7r^/2) , n = l,...,N , (5.22) 

where we have defined 

We note that in the weak coupling analysis summarized above f3 cannot be 
appreciably smaller than 1. However, numerical Dynamical Density Matrix 
Renormalization Group (DDMRG) computations [97, 155] indicate that the 
U(l) Thirring model description remains valid in an extended region of pa- 
rameter space where U,Vi,V2 are not small compared to t, but the model 
is still in a Mott insulating phase. A rough criterion for the applicability of 
the U(l) Thirring model to the description of the low-energy physics is that 
the charge gap should be small compared to the electronic band width At. 
When applying the U(l) Thirring model in the extended region of parameter 
space, the gap A and the parameter /3 have to be determined numerically. 



5.2.2. Quarter Filled Band 

In the quarter filled case there are no simple Umklapp processes that can 
open a gap in the charge sector. As a result the quarter-filled extended Hub- 
bard model is metallic in the weak coupling regime, i.e. both charge and spin 
sectors are gapless. However, integrating out the high energy degrees or free- 
dom in a path-integral formulation generates "double Umklapp" processes 
involving four electron creation and annihilation operators each [122,322]. 
For small U, Vj these processes are irrelevant, but increasing U, Vj decreases 
their scaling dimension. This suggests that for sufficiently large U, Vj the 
double Umklapp terms eventually become relevant and open up a Mott gap 
in the charge sector. Such a scenario is indeed supported by numerical com- 
putations [250]. Assuming that a bosonized description remains valid beyond 
the weak coupling region, the low-energy effective Hamiltonian is identical 
to (5.16). However, the relations between the Fermi operators and the Bose 
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fields are different 



RUx) 




V2^ 



exp 



exp 



-i[f*'=-|H)''''("xi'^--®-i)' 



Although the low-energy effective Hamiltonian is the same as for the half- 
filled Mott insulator, the physical properties in the quarter-filled case are 
rather different. Firstly, the insulating state emerging for sufRcicntly large 
U, Vj at quarter filling is generated by a different physical mechanism (double 
Umklapp scattering) as compared to half filling (Umklapp scattering) and 
concomitantly is referred to as a Akp charge-density wave insulator in the 
literature [232]. We adopt this terminology here. Secondly, the quantum 
numbers of elementary excitations in the charge sector are different. Like 
for the half-filled case the elementary excitations in the charge sector are 
a soliton/antisoliton doublet, but now they carry fractional charge ib|. A 
simple way to see this is to recall that the conserved topological charge in 
the SGM is defined as 



The soliton has topological charge —1 and the antisoliton +1. A simple 
calculation shows that the right moving fermion creates two solitons 



This implies that fermion number 1 corresponds to topological charge 2 and 
hence solitons have fractional charge. The elementary excitations in the spin 
sector are again a pair of gapless, chargeless spin zt^ spinous. 

5.3. Correlation Functions 

Due to spin charge separation a general local operator 0{t,x) can be rep- 
resented as a product of a charge and a spin piece O = OcOg- As a result 
correlation functions factorize as well 



(0|Ot(r,a;)O(0)|0) =e(0|Oj(r,,x)a(0)|0), s{0\OliT,x) OM\0)s , (5.27) 




(5.25) 



QRl{x)\0) = -2Rl{x)\0) . 



(5.26) 



where \0)s,c are the vacua in the spin and charge sectors respectively. Corre- 
lation functions in the spin sector are easily evaluated as we are dealing with 
a free theory. In what follows we only need correlators of vertex operators. 
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which in our normahzations are given by 

1 



5( exp {ia(f)s{T, x)) exp {-ia(j)s{0)))s = — ^ , 

[vsT — ixy^ 

i(exp (ia^s(r,x)) exp (-^^^(O)))^ = • (5.28) 



5.3.1. Correlation Functions in the Charge Sector 

As the charge sector of the U(l) Thirring model is equal to a SGM, we 
may use the form factor approach to determine the charge pieces of cor- 
relation functions. In the parameter regime of interest the spectrum of the 
SGM consists of soliton, antisoliton and several breather bound states, which 
we denote by labels s,s,Bi,B2,...Bi^. Energy and momentum are then 
parametrized in terms of the rapidity variable 9 as 

Ea{e) = Aa cosh e , Pa{e) = —smhe , (5.29) 

where 

A, = A,- = A , Ab„ = An , (5.30) 
and the breather gaps, A„, are given by (5.22). 



5.4. Optical Conductivity 

The optical conductivity was calculated in Refs. [62,63,97, 153]. In the field- 
theory limit, the electrical current operator is related to the fermion current 
J by 

Jei(T, x) = -eaoJ (r, x) = -i — — 5^$c , (5.31) 

where ^ is a non-universal constant. For U,Vj t we have ^ ~ 1. The 
expression for the current operator is the same for the half-filled and the 
quarter filled bands. As seen from Eq. (5.31), the current operator does 
not couple to the spin sector. This shows that spinons do not contribute to 
the optical conductivity in the field theory limit. Hence, the calculation of 
the optical conductivity has been reduced to the evaluation of the retarded 
current-current correlation function in the charge sector. The real part of 
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the optical conductivity (lo > 0) has the following spectral representation 

In Refs. [24,115,167,209,210,280] integral representations for the form factors 
of the current operator J(t, x) in the SGM were derived. Using these results 
we can determine the first few terms of the expansion (5.32). From (5.32) it 
is easy to see for any given frequency lo only a finite number of intermediate 
states will contribute (as per Section 2): the delta function forces the sum of 
single-particle gaps Y2j t° ^^^^ than lo. Expansions of the form (5.32) 
are usually found to exhibit a rapid convergence (Section 2.5), which can be 
understood in terms of phase space arguments [54,230]. Therefore we expect 
that summing the first few terms in (5.32) will give us good results over a 
large frequency range. 

Using the transformation property of the current operator under charge 
conjugation one finds that many of the form factors in (5.32) actually vanish. 
In particular, only the "odd" breathers Bi, B3, . . . (assuming they exist, i.e., 
/3 is sufficiently small) couple to the current operator. The first few non- 
vanishing terms of the spectral representation (5.32) are 



Re (t{uj) 



\2aoPj 



[(1+0/2?] 
n=l 



(5.33) 

Here C7s„ (lo), agsiLo) and 03^82 {^) a-'^s the contributions of the odd breathers, 
the soliton-antisoliton continuum and the B1B2 breather-breather contin- 
uum respectively. The latter of course exists only if ^ < ^. We find 



tan(IIfS) ^ 2 ^ 



X exp 



-'I 

Jo 



dt sinh(t(l -g)) sinh^(mtO" 
t sinh(t^) cosh(i) sinh2f 



. (5.34) 
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The soliton-antisoliton contribution is [62] 



4Vw2 _ 4^2 _ 2A) 



X exp 



[cos(|) +cosh(|)] 

•oo sinh(t(l 1 - cos(^) cosh2t 

'0 't~ 



sinh(t^) cosh(t) sinh 2t 



(5.35) 



where = 2arccosh(a;/2A). The result for the B1B2 breather-breather 
continuum is given in Ref. [97]. As a function of the parameter /3, the 
optical conductivity behaves as follows. 

• 1 > > 1/2: 

In this regime the optical spectrum consists of a single "band" 
corresponding to (multi) soliton-antisoliton states above a threshold 
2A. The absorption band increases smoothly above the threshold 2A 
in a universal square root fashion 



a{uj) oc \J Lo — 2 A for 



2A+ . 



(5.36) 



In Fig. 20 we plot the leading contributions for the case 0^ = 0.9. 
Clearly, the four-particle contribution is negligible at low frequencies. 



10 




Figure 20. Optical conductivity for 0^ = 0.9. Shown arc the dominant contributions at low 
frequencies: the soliton-antisoliton part <Tss{lo) and the two soliton - two antisoliton contribution 
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• 1/2 > > 1/3: 

Here the optical spectrum contains one band and one excitonic 

breather peak below the optical gap 2 A at the energy ujbi = ^i- 
The optical weight is progressively transferred from the band to the 
breather as decreases down to 1/3. The absorption band again 
increases in a square-root fashion (5.36) above the threshold for all 
values of except = 1/2, where the breather peak merges with 
the band. In this case (t{lo) shows a square-root divergence at the 
absorption threshold 

a{uj) DC ^ ^ for uj 2A+ (/^^ = 1/2) . (5.37) 

Vlj — 2A 

The field theory results discussed here have been compared to DDMRG com- 
putations (see e.g. Ref. [154] and references therein) computations of a{uj) 
for extended Hubbard models in Refs. [97, 153, 155] and good agreement 
has been found in the appropriate regime of parameters. We note that the 
DDMRG method can also deal with parameter regimes in the underlying 
lattice model, to which field theory does not apply. Let us discuss the above 
results from the point of view of an application to optical condTictivity mea- 
surements in the Bechgaard salts [139,274,305]. There it is found that up to 
99% of the total spectral weight is concentrated in a finite-frequency feature, 
which has been attributed to Mott physics of the type discussed here [274]. 
A comparison of (5.33) to the experimental data gives satisfactory agree- 
ment at high frequencies, but the detailed peak structure at low frequencies 
is not reproduced [62] . A likely source for these differences is the interchain 
tunneling. 



5.5. Spectral Function 
5.5.1. Half-Filled Mott Insulator 

The zero-temperature spectral function of the half-filled Mott insulator has 
been studied in many previous works. There have been extensive numerical 
studies on finite size t-J and Hubbard models e.g. Refs. [32,173,275,323]. 
The limit where the single-particle gap is much larger than the bandwidth 
was treated in Refs. [247,248,251]. This regime is complementary to the case 
we address here. The weak-coupling limit we are interested in was studied 
in Refs. [282,306,312], where a conjecture for the spectral function was put 
forward. Here we derive these results by means of an exact, systematic 
method. In what follows we will for simplicity fix /? = 1, i.e. deal with the 
Hubbard model only. 
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5.5.2. Zero Temperature 

The single particle Green's function is calculated by following the steps out- 
lined above [98]. The creation and annihilation operators for right and left 

moving fermions factorizc into spin and charge pieces upon bosonization 
(5.17). The spin part is easily calculated: in imaginary time we have 



ifty 

e 4 



\/vsT — ix 



(5.38) 



The correlation function in the charge sector is calculated by means of the 
form factor bootstrap approach. Taking into account only processes involv- 
ing one soliton we obtain 



g-|[<i>c+ecl g|[*c+ec 



~ = exp ( — — a/x^ + vIt"^ , (5.39) 

^/vcT-ix \ Vc J 



where Zq ~ 0.9218 [213]. The form of (5.39) is fixed by the transformation 
properties of the operator exp(— | [$c + ©c]) under Lorentz transformations. 
The corrections to (5.39) involve intermediate states with three particles 

and are negligible at long distances/low energies [80]. The imaginary time 
Green's functions of left and right moving electrons are then given by 



i?^(r,x)4(0,0)) ~- 



A + x'^vc^ 



-{Tr L,(r,x)4(0,0)) 



27r y/VsT — iX\/VcT — ix 



A^/T^+xhj^ 



27r ^/v^T^flx^/v^T^ 



IX 



(5.40) 



(5.41) 



Fourier transforming and analytically continuing to real frequencies (we sup- 
press the spin index a in the formulas below) we arrive at the following result 
for the retarded single-particle Green's function 



LO + VcQ 



\ Vc + Vs ^A2 -I- v^q^ - a;2 
(a + VA^T^^IT^)' - ^^^^(o; + i;cg) 



(5.42) 



We note that the charge velocity Vc is larger than the spin velocity Vg. The 
spectral function is obtained from the imaginary part of the single particle 
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Green's function (5.42). In the case Vs = Vc = v it takes the simple form 



1 TRW X Zo A e(U -^A2 + u2g2) 

A{u, kF + q) = - -ImG(«) a;, kp + q) = r ^' ' . 

TT 7r|a; - vq'l ^co^ - - ?;2g2 

(5.43) 

In Fig. 21 we plot the spectral function in the case Vs = 0.851 If c and 
/3 = 1 in a series of constant q scans. There is a continuum of states above the 
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Figure 21. Spectral function for a half-filled Mott insulator with Vs = 0.8511iic and fS = 1, 
corresponding to the half-filled Hubbard model with U = t. 



Mott gap, which is smallest sd, kp- The most striking feature is the presence 
of two distinct "peaks" dispersing with velocities Vg and Vc respectively. 
Most of the spectral weight is concentrated in these features, which are a 
direct manifestation of spin-charge separation. The higher (lower) energy 
feature corresponds to the situation where the momentum carried by the 
spin (charge) sector is held constant. Concomitantly the high/low energy 
feature is referred to as antiholon/spinon peak. Neither peak is sharp but has 
intrinsic width (more precisely, they correspond to square root divergences) . 
The threshold at a; > in the vicinity of A;ir is [98] 

^thres(A:F + g) = <^^ 'J (5.44) 

[vsq + AVI- iiq>Q, 
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where 

Vc VcV'"i - 

We see that for g < Q the threshold follows the antiholon dispersion, whereas 
for g > Q it follows the (linear) spinon dispersion (shifted by a constant). 



5.5.3. Finite Tem,perature T <C A 

It is possible to extend the results for the spectral function to small temper- 
atures r <C A [100] by the method described in Ref. [188] and summarized 
in section 2.7. As T ^ A the effects of temperature on the charge piece of 
the correlation function are small, but the spin piece can be affected quite 
strongly because the spinous are gapless. The spectral function in the vicin- 
ity oi kp can be represented as 



/oo 
-oo 



+ e 



-Ac{z)/T 



gs{uj - Ac{z),q s{z)) 

Vc 

A 



gs{uj + Ac{z),q-\ s{z)) + 

Vc 



(5.46) 



where A 



Vc (27r) 



c{z) 



coshz, s{z) = sinh^; and 



gs{u},q) = i/^^Re 



' 4 AttT ' 2 



5{lo - Vsq) . (5.47) 



In Fig. 22 we plot the spectral function for Vs = 0-5vc and /3 = 1 at a 
temperature of 0.05 times the single particle gap. A significant temperature 
dependence of the "charge" peak is apparent. It gets damped rather strongly 
at temperatures that are still small compared to the gap. This may make 
an unambiguous detection of SC-separation by ARPES more difficult as 
the experiments are done at elevated temperatures (room temperature for 
Sr2Cu03) in order to avoid charging effects. 



5.5.4. Quarter- Filled CDW Insulator 

The single-particle Green's function in the quarter-filled case can be deter- 
mined by the same method [99] employed in the half-filled case. The spin 
sector is again gapless and the spin part of the Green's function is easily de- 
termined. The charge piece is analyzed by means of the form factor bootstrap 
approach. The difference to the half-filled case is that the charge piece of 
the single-electron operator (5.24) now couples to at least two (anti)solitons. 
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Figure 22. Spectral function for a half- filled Mott insulator with Vg = O.Sfc, T = (dotted) and 
T = 0.05A (solid). The different curves are constant-q scans with Vcq/^ = —1.5, —1, . . . 2.5 (from 
bottom to top), which have been offset. 



Neglecting contributions of four or more particles in the charge sector we 
obtain the following result for the single-particle Green's function in the 
vicinity oi kp [99] 



G^^\oj,kF + q) = -C 



l + a 



°° ,^ E{2e)smh^{9) uj + Vcq 
do 



1 — a 
l + a 



'u + Vcqf 



(5.48) 



where E{6) is given in (3.81), = uj"^ — v^q^, c{0) = 2Acosh^, and a = 
Vs/vc- C is a numerical constant. In Fig. 23 we plot the spectral function 
Aiuj^kp + g) as a function of uj for Vg = 0.8fc and different values of q. 
We observe that the spectral function is rather featureless and there are 
no singularities. Furthermore, in contrast to the half-filled Mott insulator 
discussed above, there are no dispersing features associated separately with 
Vc and Vs- The absence of any distinct features is clearly related to the fact 
that the electron has "fallen apart" into at least three pieces. 

Just above the threshold at = — v^q^ = 4A^ one finds for Vc = Vg 



Aiuj^ q) oc 



1 



{— 

vq\ V2A 



(5.49) 
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Thus the spectral weight increases hnearly with s — 2 A above the threshold. 



+ 
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CO/A 



Figure 23. A{(jj, kp + q) as a function of w/ A for Vs = O.Svc- The curves for different q are offset. 



The tunneling density of states is 



2C 



arccosh(u;/2A) E (20) Sinh'^ (9) 

dd . 

V cosh 6 



, , , , (5.50) 

i^-\/VcVs Jo Vcosh 9 

and displays a roughly linear increase after an initial (a; — 2A)^/^ behavior 
just above the threshold at uj = 2A. 

Let us now turn to a comparison with experiments. PrBa2Cu307 (P123) 
is a quarter-filled quasi-lD cuprate, to which the theory presented here may 
have some relevance. The ARPES data for P123 (Fig. 3 of [222]) show a 
single, very broad, dispersing feature that is asymmetric around kp. If we 
interpret the underlying increase in intensity in the data as background, the 
signal has a form similar to Fig. 23. In order to assess whether the theory 
presented here is indeed relevant to P123, it would be interesting to extract 
(-) 

a value Ap^ for the gap from the ARPES data and compare it to gaps seen 

(-) (-) 
in optical measurements A^p^ and the thermal activation gap Ar^ extracted 

e.g. from the dc conductivity. The theory presented here predicts 



A, 



(1) 

opt 



_ OA^^) _ /\(4) 

.ii-lrp PE 



(5.51) 
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This is in contrast to the case of the half-filled Mott insulator, where one 
has 

Aij; = 2A^^^ = 2Ag . (5.52) 
5.6. A Remark on Luttinger's Theorem 

In all cases we have discussed, the Green's functions have branch cuts but 
no poles. In particular there are no poles at zero frequency and hence no 
Fermi surface. Nevertheless Luttinger's theorem holds as we will now show 
following [100]. Luttinger's theorem reads [2,85] 

where N is the total number of electrons, V is the volume and the integration 
is over the interior of the region defined by either singularities or zeroes of 
the single-particle Green's function. The former is the case for a Fermi liquid 
whereas the latter is the case at hand. The interest in the equality (5.53) 
derives from the fact that it implies that the integral on the right hand 
side is independent of electron-electron interactions. In a Fermi liquid this 
means that the volume of the Fermi surface is unaffected by interactions. 
The Green's functions we have discussed above all fulfill (5.53) by virtue of 
them having zeroes at the position of the non-interacting Fermi surface. 

The mechanism underlying this fact can be understood as follows. For 
simplicity let us concentrate on the Lorentz invariant case, where Vg = Vc- 
Then 

• As the spin sector is gapless, we have (i?o-(T, x)lJ-(O)) = 0; 

• Lorentz invariance of the low-energy effective theory implies that 

(^'^(r, x)^liO)) ~ exp(±i(/))7^(r); ^ = R,L. (5.54) 

Here r and (p are polar coordinates and TZ denotes the radial part of 
the correlation function. 

As we are dealing with an insulating state we have TZ{r) oc exp(— Ar) at 
large distances and hence J dr TZ{r)r is finite. Thus 

G(0, ±kF) = J exp(±z(/)) j dr n{r)r = . (5.55) 

For a metallic state the r integral would diverge and the Green's function 
would have a singularity rather than a zero. 
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5.7. Interchain Tunneling 

Let us now consider a quasi one dimensional situation of Hubbard chains 
weakly coupled by a long-ranged interchain tunneling 

I ly^m;n;a 

= E + h-c. + uY: -SnS • (5.56) 

n,a n 

We allow the interchain tunneling to be long-ranged in order to be able to 
carry out a controlled expansion in the case where the Fourier transform 
t±{k) of the interchain tunneling becomes of the same order as the ID Mott 
gap A (see the discussion below) . For simplicity we take t± long-ranged only 
in the direction perpendicular to the chains, although it is straightforward 
to generalize all formulas to an interchain tunneling of the form 

E ML'^cgtcM+h.c. . (5.57) 
5.7.1. Expansion Around Uncoupled Chains 

Following the analogous treatment for the case of coupled Luttinger liquids 
[18, 45, 308] we take the interchain tunneling into account in a perturbative 
expansion. The building blocks of this expansion are the n-point functions 
of fermion operators for uncoupled chains, which are represented pictorially 
in Fig. 24. The full single-particle Green's function is given in terms of 




Figure 24. Elements of the diagrammatic perturbation theory in the interchain tunneling. 

a diagrammatic expansion, the first few terms of which are shown in Fig. 25 
Unlike for Luttinger liquids, it is extremely difficult to determine four- 
point functions for ID Mott insulators. On the other hand, it is straightfor- 
ward to sum all diagrams involving only two-point functions of uncoupled 
chains. This approximation is known as the random phase approximation 
(RPA) and has a long history [268]. The small parameter making RPA a 
controlled approximation for any form of the interchain tunneling is the ra- 
tio of the interchain tunneling to the Mott gap (at zero temperature) [45]. 
This is because the only energy scale entering the n-point functions is the 
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+ 

Figure 25. Diagrammatic expansion for the single-particle Green's function of coupled chains. 

dynamically generated gap A. Dimensional analysis then shows that the 
diagrams neglected in RPA are suppressed by powers of t±/A. Within RPA 
the single-particle Green's function is given by 

GsD (W, g, k) = -^r— ; , 

l-t^{k) GiD{uJ,q) 

t±{k) = Y,t'reMik-[Ri-Rm]) ■ (5.58) 

m 

The RPA Green's function G3d(w,(/, /c) has the interesting property that it 
has a pole, which corresponds to a bound state of an antiholon and a spinon 
with the quantum numbers of an electron. For small interchain tunneling 
this bound state still has a gap. Fig. 26 is a density plot of G3D(w,g,A;) 
as a function of uj and q for a fixed value of k and hence a fixed value of 
t±{k). At energies above the Mott gap there is a continuum of states, which 
is similar in nature to the result for uncoupled chains. Below the Mott gap 
the coherent electronic mode is visible. 



5.7.2. Formation of a Fermi surface 

As long as the "binding energy" of the electronic bound state is small, RPA 
is a controlled approximation for any form of the interchain tunneling [45]. 
However, in the most interesting situation when the gap of the bound states 
becomes very small, RPA becomes uncontrolled; there is no longer any small 
expansion parameter for a generic t±(k). An exception is the case of a long- 
ranged interchain tunneling: here the support of t±{k) becomes very small, 
so that any integration over the transverse momentum generates a small 
volume factor proportional to the inverse range of the hopping. Recalling 
that RPA takes into account all terms not involving any integrations over the 
transverse hopping (i.e. "loops"), we conclude that RPA is the leading term 
in a controlled loop expansion. Increasing t±{k) in RPA reduces the gap 
of the electronic bound state, until it eventually vanishes and electron and 
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Figure 26. The spectral function for fixed momentum perpendicular to the chains k. The elec- 
tronic bound state is visible as a delta-function peak below the antiholon-spinon continuum 



hole pockets are formed: we have crossed over from weakly coupled ID Mott 
insulators to an anisotropic Fermi liquid. As an example, let us consider a 2D 
square lattice with interchain tunneling between nearest neighbor chain only. 
Here RPA is an uncontrolled approximation, but we still find it instructive 
to discuss its predictions. In this case electron pockets are formed in the 
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Figure 27. Fermi surface predicted by the RPA for a 2D square lattice with nearest-neighbor 
interchain tunneling. 
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vicinity of the points (±^^,0) and hole pockets form around (±fcF,zb7r). 
In the particle-hole symmetric case the volume of the electron pockets is 
precisely the same as the volume of the hole pockets. 

6. Hubbard Ladders and Carbon Nanotubes: SO(8) 
Gross Neveu model 

In this section we address another set of examples of quasi one-dimensional 
interacting electronic systems, two-leg Hubbard ladders and single- walled 
carbon nanotubes. 

Two-leg ladders have been the focus of much theoretical (i.e. [22, 27- 
29, 33, 67, 68, 102, 123, 156, 157, 184, 185, 203, 204, 234, 258, 270, 273, 307, 311]) 
and experimental activity (i.e. [142-144,157,196,219,267]). At half-filling 
they are Mott insulators, exhibiting gaps to all excitations, and in par- 
ticular, a spin gap. They are typical examples of 'spin-liquids'. Upon 
doping, the gaps to all excitations except for those with charge-two sur- 
vive [27,102,184,203,204,273]. The gapless charge modes induce quasi 
long-range superconducting pairing correlations, with approximate d-wscve 
symmetry, reminiscent of underdoped cuprate superconductors. 

Carbon nanotubes are novel materials whose mechanical and electronic 
properties promise potential for new technological applications [86]. They 
are formed by wrapping graphite sheets into cylinders of nanoscale dimen- 
sions. They support electronic excitations, which, for a prominent mem- 
ber of the nanotubc family, the armchair (n,n) type, can be described by 
the same theoretical model as that used for the two-leg Hubbard ladders 
[28,189,203,204,235]. Even though these systems would be one-dimensional 
band metals in the absence of interactions, they become Mott insulators at 
half-filling due to the presence of short-ranged electronic interactions, which 
play an important role due to their onc-dimcnsional nature. ^ It is these 
interaction effects that we analyze exactly in this section using form factors. 

After experimental techniques had been developed to fabricate long 
single-walled nanotubes with high yields in the laboratory, this field of ma- 
terial science has seen a explosive development [42,81,255,314]. Electronic 
properties can be measured relatively easily by attaching metallic leads [296] 
or by tunneling into these materials with scanning tunneling microscope 
(STM) tips. The practical feasibility of such tunneling experiments from 
an STM tip into an individual single-walled nanotube placed upon a gold 



1 1f forward scattering interactions axe assumed predominant, the metallic nanotube is a Luttinger 
liquid, not a Mott insulator [88,89,158,200,204,273]. 
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substrate (screening long-range Coulomb forces) has been demonstrated in 
a number of experimental reports [159, 160,304,314]. 

As said, both systems, two-leg ladders and single-walled armchair nan- 
otubes, would be one-dimcnsional metals in the absence of electron inter- 
actions. These are described theoretically (on scales much smaller than the 
non-interacting band width) by two species of spinful massless Dirac fermions 
in (1 -|- 1) dimensions. (These two species have equal Fermi velocities due to 
particle- hole symmetry present at half- filling.) For the ladder compounds, 
the two species arise in an obvious way from the two legs of the ladder, 
whereas for the nanotubes they arise from the particular band structure of 
the underlying hexagonal graphite lattice, characterized by two Fermi points 
in the Brillouin zone [28,189]. These massless Fermi surface excitations in- 
teract with short-range interactions whose detailed nature is determined by 
non-universal microscopic considerations. 

A remarkable low energy feature of these system was observed by Lin, 
Balents, and Fisher [204]. These authors argued that within an 1-loop RG 
any such model with generic, non-chiral, short range interactions flows at 
half- filling into a theory with an immense symmetry, namely the S0{8) sym- 
metric Gross-Neveu model. This model not only has a large SO (8) global 
symmetry, which encompasses an onc-dimensional version of 50(5) advo- 
cated by S. C. Zhang [335], but in addition has an infinite number of hidden 
conservation laws, which are a consequence of the model's integrability. 

Motivated by this work, we study in this section the low energy properties 
of nanotubes and two-leg ladders through computing correlation functions 
via form factors in SO (8) Gross-Neveu. In order to be pedagogical, we ini- 
tiate our discussion by explaining why the low energy excitation spectrum 
of armchair carbon nanotubes and Hubbard ladders (four Dirac fermions) 
are identical. We then move on by reviewing the RG argument (and its lim- 
itations) by which these systems flow onto the S0(8) Gross-Neveu model. 
Having done this, wc consider the excitation spectrum of SO (8) Gross-Neveu 
together with the connection between the fields in the theory and the orig- 
inal lattice operators. Finally we compute a number of physically relevant 
correlation functions for these systems. 

6.1. Armchair Carbon Nanotubes and Hubbard Ladders: 
Identical Low Energy Behavior 

We begin by showing that the set of low energy excitations in a non- 
interacting Hubbard ladder is equivalent to the low energy sector in an 
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Figure 28. The bonding and anti-bonding bands of a Hubbard ladder. 

armchair carbon nanotube [204,205]. The Hubbard ladder Hamiltonian is 

Ho = ta[ai^ + + tal^{x + l)a2a{x) 

+ t±a\^{x)a2a{x) + h.c.). (6.1) 

Here the ai/a\ are the electron annihilation/creation operators for the elec- 
trons on rung / of the ladder, x is a discrete coordinate along the ladder, and 
a =t,i describes electron spin, t and t± describe respectively hopping be- 
tween and along the ladder's rung. The first step in the map is to reexpress 
the a's of Hq in terms of bonding/ anti-bonding pairs: 

Cja = + (-lF«2a) ■ (6.2) 

With this transformation, the Hamiltonian can be diagonalized in momen- 
tum space in terms of two bands, the bonding (B) and anti-bonding bands 
(AB). Their dispersions are given by 

Es/ABikx) = -2tcos{kx) =F i± ■ (6.3) 

These bands are pictured as half-filled in Fig. 28. Working at this filling, 
particle-hole symmetry dictates that the Fermi velocities, vpj, of the two 
bands, j = 1,2, are equal. As we are interested in the low energy behavior 
of the theory, the Cja's are linearized about the Fermi surface, kpj: 



Cja^CRjaC *J +CLjae ^ , 



(6.4) 
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where L, R correspond to right and left moving modes. With this Hq be- 
comes, 

Ho = VF j dx^^ [cRjJdxCRja - C^^JdxCLja] ■ (6.5) 
ja 

Thus the lower energy description of a non-interacting Hubbard ladder is 
that of four Dirac fermions (once spin degeneracy is taken into account). 
We now show that the low energy spectrum of a non-interacting armchair 
carbon nanotube is exactly the same. 




Figure 29. The graphite lattice is composed of two interpenetrating triangular lattices. 



A single- walled armchair carbon nanotube is made by rolling up a sheet of 
graphite [132,204,221,264]. As such we first focus upon the graphite lattice, 
picture in Fig. 29. This lattice consist of two interpenetrating triangular 
lattices marked as blue circles and violet squares in Fig. 29. A tight-binding 
hopping model of electrons on this lattice is described by the Hamiltonian 

H= ^ \ -taL(r)a2a(r + a_+d)-t4^(r)a2a(r + a++d) 

-txaL(r)a2a(r + d) + h.c.j . (6.6) 

Here R is the set of lattice vectors generated by the bases a_ and a+ where 
a± = a(±l/2, -v/3/2). The two triangular lattices are displaced by a vector 
d = a(0, 1/^/3). Here again Or/aJ create and destroy an electron at site r. 
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Figure 30. The dispersion of the bonding band of graphite (with a = 1). 



In the absence of interactions, it is possible to diagonalize this Hamilto- 
nian. One looks for eigenfunctions of the form 

|V'k.) = E(4<x + /3(k)aL)e^'^"- (6.7) 

reR 

Solving Hlipka) = -E'(k)|'0fc(T), one then finds the following solution 

r(k) = 2tcos{k^a/2)e'''y^/^ + t^e'^'^^"/^ ; 
i?(k) = ±|r(k)|; 

/32(k) = r(k)r*(k). (6.8) 

The two signs in the above equation for the energy, E{k), correspond to the 
bonding and antibonding bands. We plot the dispersion of the bonding band 
in Fig. 30. 

We see that the dispersion is characterized by a series of cones whose tips 
extend to zero energy. The positions of these tips are known as Dirac points. 
The Dirac points are located at k = (ib47r/3a,0) (and equivalent reciprocal 
lattice vectors). In Fig. 31 we mark with violet squares the position of the 
Dirac points in reciprocal space provided t = t±. At half-filling the bonding 
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band is completely filled and there will be a set of low energy excitations. 
These excitations, however, have a vanishing density of states at the Fermi 
surface. Consequently graphite is known as a semi-metal. 

At this point we want to consider the effects of rolling up the sheet of 
graphite on the low energy spectrum. These effects will depend on how 
precisely this is done. We will only be interested when the sheet is rolled up 
such that an (A'^, N) armchair carbon nanotubc is obtained. This tube has a 
circumference of \/^aN. Referring to Fig. 29, a (2, 2) tube can be obtained 
by rolling the sheet so that the two lattice points labelled (0, 0) and (2, 2) 
are identified. When the sheet is rolled up, momentum transverse to the 
direction of the resulting tube will be quantized. For an armchair tube, 
ky momentum is quantized. For an (N, N) tube, there will be N distinct 
allowed values of ky given by 

Ott r? (N — 1) 

ky = -^-, n = 0,±l,±2,... ,±^^-^, Nodd; 

fcy = ^|:, n = 0,±l,±2,--- ,± ^^~'^\ n/2, Neven. (6.9) 

The allowed values of discrete fcj^-momentum are shown in Fig. 31 for a (3, 3) 
tube. 

When the tube is rolled up, one generically expects the hopping param- 
eters to change. If once t = t±, we then expect t ^t±. With such a change, 
the Dirac points are shifted. How these Dirac points are shifted is pictured 
in Fig. 31 for a (3, 3) tube supposing that t±_ is now less than t. We see the 
shift involves only k^ and not ky. Thus rolling up an armchair nanotubc 
leaves the Dirac points (and so the presence of low energy excitations) in the 
allowed spectrum. The linearization of the spectrum about each of the two 
Dirac points gives a single Dirac fermion. Thus with spin degeneracy, the 
low energy spectrum of an armchair carbon nanotube is four Dirac fermions, 
identical to that of a Hubbard ladder. 



6.2. Weak Coupling Flow Onto SO(8) Gross Neveu 

There exist a number of treatments of interactions in carbon nanotubes. 
As carbon nanotubes are quasi-one dimensional, a number of authors have 
employed bosonization, in combination with the renormalization group, to 
understand the effects of correlations [88,89,102,158,204,273]. Wc will be 
interested in treatments of half-filled nanotubes with short range interactions 
where all excitation are gapped as opposed to nanotubes with long range, 
unscreened Coulomb interactions. In the latter case, the nanotubes are 
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Figure 31. The Brillouin zone for the graphite lattice. Pictured are Dirac points for both the 
case i = (violet squares) and t > tj_ (green circles). The blue dashed lines mark the allowed 
values of quantized fcy-momenta for a (3, 3) armchair carbon nanotube. For the armchair tubes, 
the Dirac points always lie at an allowed k value. 



Luttinger liquids [88,89,158,200,204,273]. In particular, we will focus upon 
Ref. [204] where it was shown that weak interactions lead the interacting 
electrons to flow in an RG sense onto the SO (8) Gross-Neveu model. 

To set up this analysis, we consider the low-energy, non-interacting 
Hamiltonian for the ladders /tubes: 

Ho = VF i dx ^ ['^^RjJd^r.^Rja " c\jJdxCLja\ ■ (6.10) 
a j = l,2 

We now consider adding all possible allowed interactions to this Hamiltonian. 
To organize this addition, we introduce various SU (2) scalar and vector cur- 
rents ^ 

^ij = 4a^cra'Cjcr' hj = -jc\^{e<T)au' Cja' ■ (6.11) 

Here a are the Pauli matrices and e^^' is the e-tensor. The crystal- 
momentum conserving interactions divide themselves into forward and back- 
ward scattering terms: 

Hp = ^ fijJuJjj - fij^fi ■ ^ij ■ (6.12) 

Here / and b are the forward and backward scattering amplitudes. From 
hermiticity and parity, we have 612 = 621 and /12 = /2i- Furthermore at 
half- filling, 611 = 622- Thus we obtain six independent amplitudes. The 
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Umklapp interactions take the form 

Hu= J: <M^- E + (6.13) 

ij=l,2 i7^j=l,2 

where 1 = 2, 2 = 1. As lij, lij are anti-symmetric, one may choose u\2 = U21 
leading to three additional independent couplings. For this set of generic 
interactions, the one-loop RG equations were derived in Ref. [204]. They 
show that Ho + i^int flow onto the S0(8) Gross-Neveu model. 

To express this remarkable result compactly we must invoke a change to 
variables. We begin by bosonizing the c's: 

Cpja = KjaC'^''^" , P = +, - = R, L . (6.14) 

Here Klein factors satisfying 

{i^joc, i^ip} = '^^ij^ap ■ (6.15) 

In terms of these four Bose fields, four new Bose fields are defined (effectively 
separating charge and spin): 

^pi = -j{^p\-\ + ^pii + 4>P2] + ^P2i) ; 

0P2 = -ji^pi-i - 4>pn + ^P2T - ^P2i) ; 

= -j^i^p^ - (t>pn ~ ^p'^i + ^p'^i) '■> 
p 

(j)p4 = — (^Pl-f + (t)pil - (t)p2\ - 4>P2l)- (6.16) 

Note that 0P4 has a relative sign between the right and left movers. This sign 
effectively masks the SO{8) symmetry in the original Hamiltonian. From 
these chiral bosons, one can define pairs of conjugate bosons in the standard 
fashion 

fi = (t>Ri + (l>Li ; 

Oi = (t>Ri - (l>Li , (6.17) 

which obey the commutation relations, 

[ip{x),e{x')] = -iA7:e{x - x') . (6.18) 

where Q{x — x') is the Hcavisidc step function. 

In terms of these variables the free part of the Hamiltonian can be written 

i^o = g E ((^-^«)' + i^-^af] . (6.19) 

a 
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The interacting momentum conserving interacting term can be written as 

1 4 
a=l 

-26^2 cos(6'4) cos(6'2) + cos(6'2)(26^i cos(6'3) + 2/^2 cos(<^3)) 
- COS(04)((6?2 + 46?2) cos(03) + {b^l2 ' 46?2) cos((^3)) , (6.20) 

where the coefHcients Aa are equal to = 2(0^*^ it/fg) and A2/3 = — (cjj^zt 
/f2)/2. The Umklapp interactions in these variables takes the form 

Hu = -16<i cos(6'i) cos(6'4) - 4<2 cos(6li) cos(6l2) 

- cos(6'i)(2«2 + 4<2) cos(6'3) + 2«2 - M2) cos(<^3)) . (6.21) 

As demonstrated in Ref. [204], the Klein factors combine in such a fashion 
so as to be equal to 1. 

For gcnerically repulsive scattering amplitudes, the various couplings /, b, 
and u flow to fixed ratios under the RG [204] : 

K2 = \K2 = /f2 = = 2<i = 2<2 = \^l2=9>^- (6.22) 

With these values, the interaction Hamiltonian dramatically simplifies to 

-f^int = Hb + Hp + Hu 

4 4 

= -^Xl^^'^-R«^^'^^»-^5 cos(^„) cos(06) . (6.23) 

a=l a^b=l 

We now write this Hamiltonian in fermionic form. We thus refermionize the 
bosons ^pa, a = 1, • • • ,4, i.e., 

*Pa = Kae''^", a =!,•••, 4; 

^PA = Pi^Ae'^\ (6.24) 
where the Klein factors are given by 

«1 = K2T> K2 = «1T> «3='«U, Ki = K2l. (6.25) 

We then find the free Hamiltonian can be written as 

Ho= f dxY, [^li^^x^aL - '^'ipid^'^aR] , (6.26) 
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where the Fermi velocity, vp, has been set to 1, while the interaction Hamil- 
tonian can be written as 

4 

Hint = - (6-27) 

a=l 

This is, of course, Hint for the 50(8) Gross-Neveu model. 

It will sometimes prove convenient to recast the theory in terms of Ma- 
jorana fermions, ipaP- In terms of the Dirac fermions, *aP) they are given 

by 

^aP = ^(^2a,P+#2a-l,p), (a=l,...,4). (6.28) 

In this basis, Hint can be recast as 

Hint = gGtGt , (a > 6 = 1, 8) , (6.29) 
where Gp = iijjaP'4'bP is one of the 28 SO{8) Gross-Neveu currents. 

6.3. Limitations of the RG Analysis 

In the previous section we have sketched the argument by which Ref. [204] 
demonstrated that generically interacting Hubbard ladders/armchair carbon 
nanotubes behave at low energies according to the SO (8) Gross-Neveu. In 
this section we want to consider how precise a statement this is. 

The analysis in [204] is based on an 1-loop RG. As such the initial micro- 
scopic (bare) interactions must be small enough so that the integrable S0{8) 
invariant RG trajectory is approached sufficiently closely after a number of 
RG steps, before leaving the range of validity of the 1-loop RG equations. 
Whenever this is the case, it can be argued that the integrable model is 
approached independently of the (sufficiently weak) values of the bare inter- 
actions. The situation for the 2-leg ladder is thus similar in spirit to that 
of the point contact device encountered in [107-109], where only a single 
operator was relevant, and this relevant operator was integrable. In the lat- 
ter case all other interactions were irrelevant in the RG sense, and could in 
principle be treated perturbatively. 

The requirement of the RG that the interactions be short-ranged is nat- 
ural in the case of the Hubbard ladders. However it may not seem so in 
the case of the carbon nanotubes. Both theory [88, 89, 158] and experi- 
ment [43,61] have discussed the case where long-ranged Coulomb forces drive 
Luttinger liquid behavior in single-walled carbon nanotubes. However we do 
not have such situations in mind for the paper at hand. Rather we want to 
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consider situations such as those found in the experiments [314] where the 
long range forces are screened. 

Although the restriction to such experiments in the case of the carbon 
nanotubes places us upon safe ground, it is not inconceivable that experi- 
ments where the long-ranged forces are present would nevertheless see be- 
havior indicative of the SO (8) symmetry. An unscreened force translates 
into an unusually large bare coupling (in comparison with other bare cou- 
plings) in the forward scattering direction. However this docs not mean the 
RG is inapplicable. The RG still indicates a potential enhancement in the 
symmetry. Because of the large bare coupling, the RG must be run a longer 
time before any enhancement would be seen but nevertheless an enhance- 
ment may well occur at some low energy scale. In terms of the experiments 
in [43,61], this would mean that at medium energy scales, Luttinger liquid 
behavior would predominate, while at much lower energy scales, SO (8) be- 
havior would be expected. However at current standing, the material science 
is not advanced to the point where it is possible to accurately probe the very 
low energy behavior. But the potential for advancement in this area is ever 
present. 

The RG analysis further requires the bare couplings to be weak. With 
Hubbard ladder compounds, this condition will not be generically met, al- 
though it certainly will not be universally violated. However with (N,N)- 
armchair carbon nanotubes, the bare couplings are naturally weak. It is 
one of the hallmarks of the physics of the (N,N) armchair carbon nanotubes 
that the electrons are delocalized around the circumference of the tube. This 
in turn leads to a scaling of the effective short-ranged interaction by 1/N, 
making it naturally small [28]. 

It can, however, be questioned on a more fundamental level whether an 
1-loop RG adequately describes the system's behavior. Difficulties with the 
analysis in [204] take two forms. As a first objection, the authors of Ref. [21] 
point out that an RG flow can imply a symmetry restoration which in fact 
does not occur. As an example they consider a U(l) symmetric Thirring 
model, 

C = ^a^d^^a + IgwUz? + \9±[{j.f + Uyf] , (6.30) 

where ja = Although the 1-loop RG equations for this model seems 

to indicate a generic symmetry restoration to a more symmetric SU(2) case 
(i.e. g± = g-y), this in fact only occurs in a certain region of coupling space. 
For TT — [gixl > — > \g±\ > 0, the U(l) model maps onto the sine-Gordon 
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model with interaction cos{(3(p) [161], where /3 is given by 

Stt/?^ = Stt — 8jU; fi = cos~^[cos{g\\)/ cos{g±)] . (6.31) 

The value of /? completely characterizes the model. While g± and flow 
under the RG, the particular combination of these parameters forming (3 
does not. Thus for this particular region of parameter space the model 
moves no closer to the SU(2) symmetric point under an RG flow. In other 
regions however (for example \g±\ > \g\\\), the situation is better; the effect 
of the anisotropy in the couplings is exponentially suppressed. 

However it is reasonably clear that such pessimism is not warranted in 
the analysis of the RG of Ref. [204]. A salient criticism of Ref. [21] is that 
in considering the action of the renormalization group, they fail to consider 
the consequences of working in the scaling limit. The scaling limit is exactly 
the limit in which a field theory becomes available. In turn, the scaling 
limit places constraints upon the possible range of bare couplings consistent 
with a field theory. In the case of sine-Gordon, the underlying integrabil- 
ity /solvability of the theory allows explicit investigation of this question. 
It is found that the allowed range is such that even moderate anisotropic 
deviations are forbidden [186]. The scaling limit, in other words, enforces 
isotropy. In cases where there are RG flows indicating an enhancement in 
symmetry (including the case at hand), this enforcement turns out to be a 
general phenomena and it leads to an expanded notion of symmetry restora- 
tion [186]. 

Wc can see however that this conclusion is a double edged sword. Strongly 
anisotropic models will not have sensible continuum limits and so will not 
be able to be described by a field theory. If the Hubbard ladder/carbon 
nanotube is strongly anisotropic, the continuum theory is not a good start- 
ing point rendering the RG analysis irrelevant. However if interactions are 
only weakly anisotropic, the continuum theory can be used and the weakly 
anisotropic Hamiltonian is better able to flow under the RG to the SO (8) 
symmetric theory. 

On a more concrete level, the breaking of the SU(2) symmetry considered 
in Ref. [21] is a rather special case. Ultimately, the parameter, (3, in the sine- 
Gordon model is protected under an RG flow by the presence of a quantum 
group symmetry arrived at by deforming a Yangian symmetry present at the 
SU(2) point. There is, however, no such known way to deform the Yangian 
in 50(8) Gross-Neveu. Indeed the natural generalization of the sine-Gordon 
model to 50(8) is not to 50(8) Gross-Neveu but to an affine toda 50(8) 
theory where such a deformation of the Yangian symmetry is possible [36]. 
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Another question that one must ask in looking at the analysis in Ref. [21] 
is how the choice of the symmetry breaking terms affects the symmetry 
restoration. The sine-Gordon model still possesses a U(l) symmetry. How- 
ever it is certainly possible to consider perturbations that break this U(l). 
Such perturbations would destroy the quantum group symmetry of sine- 
Gordon and thus might lead to symmetry restoration. This would be per- 
haps closer to the RG analysis of Ref. [204] where a large number (nine) of 
marginal perturbations were included. We in fact consider exactly such a 
situation in Ref. [186] and find that indeed there is symmetry restoration. 
We also note in passing that the authors of Ref. [21] consider an anisotropic 
Gross-Neveu model, a model of direct relevance to the situation at hand. 
They conclude through a mean-field/largc N limit computation that the 
model is intrinsically anisotropic thus throwing doubt upon the analysis in 
Lin et al. [204]. However we would point to how the bare couplings are 
scaled in the large N limit. The anisotropic model they consider has three 
bare couplings. One is chosen to not scale at all, one scales as 1/A^, and the 
last scales as d <1. With this scaling [21] the model is intrinsically 

anisotropic. However this model possesses a diverging bare anisotropy in the 
large N limit while the RG is not allowed to act. As such, we believe this 
example is not so directly telling. 

The second objection to the analysis of Ref. [204] is its omission of chiral 
interactions that alter the Fermi velocities [87]. Such interactions, although 
they are absent from the 1-loop RG, likely play a role at higher order. How- 
ever their effect is less drastic than envisioned in Ref. [87] . There a scenario 
was considered where the invariant RG trajectory of higher symmetry was in- 
herently unstable to perturbations. However the 50(8) RG ray in Ref. [204] 
has a basin of attraction of finite measure. The effect of chiral interactions 
is to then slightly alter the direction of the ray. In turn the ratio of masses 
of the various fundamental excitations will be slightly perturbed away from 
one. 

In taking into account of these objections, prudence suggests a modifi- 
cation in the understanding of the RG analysis of Ref. [204] . This analysis 
tells us that while the RG flow does not restore an exact symmetry, it leaves 
us close to the symmetrical situation. In particular, it indicates that while 
the masses in the actual system may differ from their 50(8) values, they do 
not wildly diverge. One then understands the 50(8) Gross-Neveu theory, 
not as precisely representative of the actual system, but in near perturbative 
vicinity of it, that is, as an excellent starting point about which to perform 
perturbation theory in the non-integrable interactions breaking 50(8) in 
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much the same spirit as done for a non-critical Ising model in the presence 
of a magnetic field [71] . 

We can state with some confidence that there are scenarios where the 
effects of these perturbations will be small. A related Hubbard lattice model 
possessing an SU{4) symmetry has been studied in Ref. [22]. Provided that 
the interaction strength was not overly large {U/t < 3), it was found using 
quantum Monte Carlo that the low energy sector of the S'C7(4) theory was 
enhanced to 5*0 (8) Gross-Neveu. In particular it was found that the charge 
gap, the spin gap, and the single particle gap satisfy the corresponding 50(8) 
ratios. 

6.4. Excitations and Physical Fields in SO (8) Gross Neveu 

Wc now consider the basic excitations of SO (8) Gross-Ncvcu together with 
identifying the fields in the theory that correspond to operators of physical 
interest. 

6.4.1. Excitation Spectrum 

The Gross-Neveu 50(8) model has an exceedingly rich spectrum. There are 
24 fermionic particles of mass m organized into one eight dimensional vector 
representation and two eight dimensional spinor representations. We denote 
the particles of the vector representation by A^, a = 1, . . . , 8. The AaS are 
the Majorana fermions of Eq. (6.28). We will often refer to these particles 
as Gross-Neveu fermions. The kink particles, in turn, will be denoted by 
Aa- Here a is of the form a = (±1/2, ±1/2, ±1/2, ±1/2) and so takes on 16 
values. These 16 particles decompose into the two eight-dimensional spinor 
representations. The division is affected by the chirality (cither even or 
odd) of the kinks. An even chirality kink has an even number of +l/2's in 
its corresponding a while an odd chirality kink has a correspondingly odd 
number. 

Beyond the eight dimensional representations, there are 29 bosonic par- 
ticle states of mass \/^m, transforming as a rank-two tensor of dimension 
28 and a singlet. Together they form a representation of the SO(8) Yangian 
symmetry. These particles can be thought of as bound states of either two 
kinks or two fundamental fermions. 

As 50(8) is a rank 4 algebra, the 50(8) Gross-Neveu model has four 
Cartan bosons (i.e. the (pPa , o = 1, ... 4) and so its excitations are char- 
acterized by four quantum numbers, Ni, i = 1, ... ,4. With the Majorana 
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fermions, the combination 

A2a ± iA2a-l , (6.32) 

carries quantum number A^^ = il, = 0,b ^ a. The quantum numbers 
carried by the kinks are directly encoded in a. If a = (01,02,03,04), 
Oj = ±1/2, the Aa carries the quantum numbers, iVj = Oj. The quantum 
numbers carried by the rank two tensor states can be directly deduced from 
the particles forming the bound state. As we will always think of the bound 
states in this way, we will not list their quantum numbers directly. 

The last thing needed in the section is to identify the relationship between 
the quantum numbers, Ni, and the physical quantum numbers of the system: 
the z-component of spin, Sz, the charge, Q, the difference in z-component 
of spin between the two bands, 5*12, and the "relative band chirality", P12, 
defined as Pu = Njn — Nn — A'^r2 + -^L2, where Npj is the number electrons 
in band j with chirality P. We have 



(A^i = 1,0,0,0) ^ 


(Q 


= 2,S, 


= 0,Su 


= 0,Pi2 


= 0); 


(0,7V2 = 1,0,0)^ 


iQ 


= 0,5, 


= 1,^12 


= 0,Pi2 


= 0); 


(0,0,iV3 = l,0)^ 


iQ 


= o,s. 


= 0,^12 


= l,Pl2 


= 0); 


(0, 0, 0, Ni = l)^ 


iQ 


= o,s. 


= 0,^12 


= 0,Pi2 


= 2). 



(6.33) 



With this assignment, we can see that the vector representation of funda- 
mental fermions corresponds to states of two electrons in the original formu- 
lation. For example, the fermion ^2±*^i carries charge ±2 and no spin (the 
cooperons), and the fermion ^4 it iA^ carries spin, Sz = 1, and no charge 
(the magnons). The spinor representations, the kinks, in turn correspond 
to single particle excitations as their quantum numbers are combinations of 
Ni/2. 



6.4.2. Relationship between Lattice Operators and Gross-Neveu Fields 

In this section we make contact between the fields of the SO (8) Gross-Neveu 
model and the original fields of the Hubbard ladders. This identification is 
crucial if we are to compute physically relevant correlation functions. 

As we have already discussed, the fundamental (Dirac) fermions of the 
vector representation are given by 

^aP = Pnae"''^'' , (6.34) 
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and carry quantum numbers corresponding to two electronic excitations. 
However the ^aP are fermionic, whereas such excitations are bosonic. As 
such, ^aP are not simply related to a fermionic bilinear of the original elec- 
trons but must be a fermion bilinear multiplying some non-local field (a 
Jordan-Wigner string). As we will not compute correlators involving such 
fields in this review, we will not elaborate upon this. 

As discussed previously, the kinks correspond to single particle excita- 
tions. Thus we expect to find that the kink fields are related to the original 
electron operators. This is true in part. There are 32 kinks in total (counting 
both left and right movers), but only sixteen electron operators, the c's and 
c^'s (four for each of the four Fermi points). So we expect only 1/2 of the 
kinks to correspond to actual electron operators. 

We represented the fundamental fermions in terms of the four Cartan 
bosons. There is a corresponding representation for the kink fields: 

^aP-e*"-^^, (6.35) 

where ^ = (^1,^2,^3,^4)- The kink fields that then correspond to the 
electron operators c's are as follows: 



CRII ~ e^('f'iR+'t'2R+'p3R+^4R)/^ 

CR2'\ ~ g'^(.'t>lR+<t>2R-4>3R-<p4R)/2 

CR21 ~ Q^'^'t>iR-'t''2R+'t>3R-'t>iR)/'^ 

CR\1 ~ QK4>iR-'l>2R-4'3R+4>AR)/'2 

(even chirality) 



CLIT e*(<Ali'+"?^2L+<i^>3L-<?>4L)/2 

CL2T ~ Q^i't>iL+4>2L-<t>3L+<t>4L)/'^ 

CL21 ~ Q^'^'t>lL-'t'2L+'t>3L+<t>AL)/'i 

CLli ~ e'^('t'iL~'l'2L-'t>3L-<p4L)/'^ 



(odd chirality). (6.36) 

With hermitian conjugates, this totals to sixteen fields. The ~ sign is meant 
to indicate that these equivalences hold up to Klein factors. The cpjaS, of 
course, are fermionic. However the kink fields as defined are not. 
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The last set of fields that are of concern to us are the currents. The 
electric current of the ladder has the lattice representation 



J ^ -i^ al^{x)aia{x + 1) - al^{x + l)aia{x) 



la 



(6.37) 



where we have summed over the contribution coming from each spin (a) 
and each leg {I) of the ladder. Taking the continuum limit, J equals, in 
Gross-Neveu language. 



J ~ z sin kpi dt(f)i ~ G12, 
where G12 is one of the SO{8) currents discussed in Eq. (6.29). 



(6.38) 



6.5. Form Factors in SO (8) Gross-Neveu 

Having identified the fields that we need to compute physical correlators, we 
now turn to the corresponding form factors of these fields. Here we only state 
the needed form factors. Their derivation may be found in [185] and [167]. 

6.5.1. One and Two Particle Current Form-Factors 

At the two particle level, two kinks of the same chirality or two GN fermions 
can couple to one of the Gross-Neveu current operators (see Eq. (6.29)). The 
corresponding form-factors are 

^f^{e^M = {Gf{o)A,ie2)Aa{e,)) = iAoidacSM - SadSbc)U{ei,e2) , 



^'(,{01,92) = {G';f{0)A^{92)Aa{9i)) = i ^ (Ca"^)a/3/M(^i, ^2) , (6.39) 
where 

s{en/2) 



Uei,e2) = (e(^i+^^)/2 - (-l)''e-(^^+^^)/2) 



ciOu) - 1/2 



X exp 



dx Grixj . / X , , „ , , 



GA^) = 2<^''>-'^f (6.40) 

c{x/2) 

Here the current G^'' is a particular combination, G"^ — {—1)^G'^'', of the 
chiral currents of Eq. (6.29). Aq is some arbitrary (real) normalization con- 
stant. 
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The current operators also possess one particle form factors. The currents 
are able to couple to one of the 29 bound state excitations. The correspond- 
ing form factor can be computed following the discussion in Section 2.6 of 
this review. We thus have 



r(2/3) \ 



-1/2 



X exp 



(6.41) 



Notice that the normalization constant Aq is the same as that which appears 
in the two-particle current form factor. 



6.5.2. One and Two Particle Kink Form-Factors 



We recall that the single electron excitations of the theory are represented 
by kink fields. A kink field, tp^, of a given chirality couples to a combination 
of excitations consisting of a kink particle of opposite chirality and a Gross- 
Neveu fermion. The corresponding form factor is given by 



= _^^e±-/^(C7")a/3/±(^i,^2), 



(6.42) 



where 

/±(^l,^2) 

Gf{x) 



; exp 



c{ei2) + 1/2 

2c(a;/6) +e-^^/6 



dx Gfix) 2 
^sm 

X S{X) 



(6.43) 



c(x/2) 

Here is an arbitrary (real) normalization parameter. 

The kink field of course also couples to a single kink. The associated form 
factor can be computed again using the bound state relations of Section 2.5 
(a kink of a given chirality is a bound state of a kink of opposite chirality 
together with a Gross-Neveu fermion). The result of the calculation is 

±/| = {mO)Aa{0)) = c±e±^/2C„^ , (6.44) 

where the constant c± is 

-1/2 

exp 



c± = —= e 



rf 



,^-(5/3)Y 



'dx Gf{x) 2 



s{x) 



s\xl^) 



(6.45) 
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Again Ap is the same normalization constant appearing in the above two- 
particle form factor. 



6.6. Exact Low Energy Correlation Functions in SO (8) 
Gross-Neveu 

In this section we study various aspects of carbon nanotubes and Hubbard 

ladders using the previously stated form factors. In particular we study the 
optical conductivity of a ladder system, the single particle spectral function 
of a ladder/nanotube, and finally the associated differential conductance 
arising from tunneling into a ladder/nanotube via a scanning tunneling mi- 
croscope (STM). We again emphasize that although we calculate these quan- 
tities with a small, finite number of form factors, the results are exact up to 
some energy scale. 



6.6.1. Behavior of Optical Conductivity in a Hubbard Ladder 

We first consider the response of the ladder system to an electric field polar- 
ized along the legs. Apart from the meanfield treatment in Ref. [204], this 
problem has been examined previously, both theoretically [138] and exper- 
imentally [143]. However these two latter papers did not consider undoped 
ladders at zero temperature. 

In linear response, the optical conductivity is given by 

Re [a{LO, k)] = Im [^i^^^l , (6.46) 



CO 



where A is the current-current correlator 

A{u;,k)= J dxdTe'^^e'='''{T{J{x,T)J{0,0)))l^-i^+s. (6.47) 
J is given by Eq. (6.38) 

J~Gi^ (6.48) 

As explained in Section 2, to compute the correlator, {T{G\'^{x, t)G\^{0, 0))), 
we insert a resolution of the identity between the two J's, turning the cor- 
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relator into a form factor sum. We then have 

{T{G\\x,t)G\\0,Q))) 

n=0 ai,--- ,an 

x{A,M---AaAOi)\Gl'm 

n n 

X exp ( - |t| ^ ruai cosh(6'i) + isign(T)x ^ sinh(6'i)) , (6.49) 
i=l i=l 

where the first sum runs over the number of particles in the form factor 
expansion and the second sum runs over the different particle types. 

We have also extracted the spacetime dependence of each term. 

We do not compute this sum in its entirety but truncate at the two 
particle level: 

{T{G\'{x,r)G\\0,0)))= J ^ {G\\0)\AUe^)) {A,2{9i)\G\\0)) 

X exp ( — |r|\/3mcosh(6'i)) + isign(r)a;\/3m sinh(0i)) 
+^ /^i^^exp ( ^ ( — |r|mcosh(^j)+zsign(r)xmsinh(^j))) 

x(Y.iG?mAl{9i)Ali92)){A,i92)A,{9,)\G\\0)^ 

ab 

+ Y,{G\\0)\Ai{9,)Al{92)){Afs{92)Aa{e,)\G\\0))) . (6.50) 

a/3 

The first term gives the single particle contribution to the correlation func- 
tion. The only particle that contributes here is ^412, denoting one of the 
particles belonging to the rank 2 tensor multiplet. At the two particle level 
a variety of contributions are non-zero. The second term in 6.50 gives the 
contribution of two Majorana fermions while the third term gives the con- 
tribution of kinks with the same chirality. 

As discussed in Section 2, this truncation of the form factor sum is better 
than it may at first seem. Because the correlator is evaluated at zero temper- 
ature in a massive system, the higher order terms make contributions only at 
higher energies, to. That is, the massiveness of the system leads to particle 
thresholds. The next contribution comes from a three particle combination 
of even kink/fermion/odd kink that carries mass 3m. Thus for lu < 3m, this 
term gives no contribution to Re[a{(jj,k)] for arbitrary k. Hence our result 
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for Re[a{u), k)] is exact for uj < 3m. 

In the case when lv does exceed 3m, we expect the higher particle form 
factors to make only a small contribution to a(uj). As we saw both in Section 
2.5 and Section 4.3.2, terms in the form factor sum involving higher numbers 
of particles make only negligible contributions to the spectral function at any 
given energy. 

Using the results for the form factors of Section 6.5, we can put everything 
together and write down an expression for Re[cr(u;)]: 



Re[cr(u;)] = 5{w 



X exp 



X exp 



X s{x) 

iO\/ — k"^ — Am? 



2 _ _ 377j,2)2 

dx Gc{x) 



lo X s{x) 
where s{x) = sinh(x), c{x) = cosh(x), and 

■2 k^ 



(u;2 - A:2)3/2 

, , ,x0i2.. 

c[x) cos( j) 

TT 



(6.51) 



712 



cosh 



-1 



UJ 



2m^ 



2m^ 



(6.52) 



As indicated in Section 6.5, Aq is an arbitrary constant normalizing all 
current form-factors while Gc{x) can be found in Eq. (6.40). 



0« 0.2 




2 2.5 

CO 

Figure 32. Plot of the optical conductivity at wavevector k = 0. 
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In Fig. 32 we plot the regular real part of the optical conductivity for 
wavevector k = 0/ We see there is an exciton type peak at w = VSm cor- 
responding to the single particle form factor contribution. The first vertical 
dashed line marks out the beginning of the two particle form factor contribu- 
tion to the conductivity. The onset of the two particle contribution behaves 
as s/uj — 2m and not as l/-\/uj — 2m as would be expected in a free theory 
due to the divergence in the density of states, the van-Hove singularity, that 
occurs in one dimensional systems. This singularity is removed by the corre- 
sponding current matrix element which behaves as (u; — 2m) with to 2m"'" 
as the low energy behavior becomes strongly renormalized in the presence 
of even weak interactions. 

The optical conductivity was computed in Rcf. [204] using the large 
limit of SO{2N), or in an alternate language, an RPA approximation. In 
such an approximation, the model becomes equivalent to a theory of four 
massive, non-interacting Dirac fermions. Hence Ref. [204] finds that the 
van-Hove singularity is present. 

The second vertical dashed line in Fig. 32 at 3m marks the point where 
the three particle form factors would begin to make a contribution. Up to 
this point, the result is exact. We note that the three particle contribution 
is strictly a consequence of interactions. In SO {8) language, two kinks of 
opposite chirality together with a fcrmion will couple to the current operator. 
In a free theory these different particles would not all exist and there would 
be no three particle contribution. 

If we were to compute the three particle contribution, three possibili- 
ties present themselves. The three particle density of states approaches a 
constant as to ^ 3m,^. If the corresponding matrix clement vanishes as 
LJ 3m"'", the contribution will open up gradually, leaving cr(u;) continuous 
at a; = 3m. If the three particle matrix element also approaches a constant 
value as uj ^ 3m"'", the conductivity will be marked by a jump at a; = 3m. 
But if the matrix element diverges in this limit, we expect to find a corre- 
sponding divergence in the conductivity at a; = 3m. Of these scenarios, it 
is our belief that the first is most likely. Moreover we expect that the to- 
tal spectral weight in such a contribution in comparison to the two-particle 
contribution will be extremely small. 



' We will not comment here whether a finite Drude weight exists beyond saying that it is a possi- 
bility as the model is integrable. 
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6.6.2. Single Particle Spectral Function 

In this section we compute the single-particle spectral function of the elec- 
trons of the ladder/nanotube. To do so we first consider the correlator, 

/oo 
dx e-'^y'-'^--{T{aUx,T)a\^{Q,m ■ (6-53) 

Here ky takes on the values 0, tt. We then define the particle/hole spectral 
functions, ^p/^, as follows: 

/oo 
dr e-'^^G{k^,ky, T)l_,_.^+g . (6.54) 
-oo 

We note that we have not explicitly summed over spin, a. 

As described in Section 6.4.2, electronic excitations around the Fermi 
point correspond in the Gross-Neveu language to low energy excitations of 
kinks. We thus expect to recast the Green's function, G, above in terms of 
kink correlators. This in fact can be done with the result, 

/oo 
dxe'''^T{cpia{x, T)4i„(0, 0))) , (6.55) 
-oo 

where i = 1,2 and kyi = (2— i)7r. The c's, the bonding-anti-bonding electrons 
are in turn related to the various kinks via Eq. (6.36). The Greens function 
on the r.h.s. of Eq. (6.55) is thus equal to 

{T{cpia{x, r)4,(0, 0))) = (T(k,V±(x, r)KamO, 0))) , (6.56) 

where a (a being its charge conjugate) is the particular kink corresponding 
to the Fermi point (kpijkyi). The Klein factors included to ensure 

the ip" are anti-commuting. Because of the 5*0(8) symmetry together with 
its associated triality symmetry, {T(Katp^ix,T)Kai^'^{Q,0))) turns out to be 
independent of the type a of kink. It is only sensitive to whether the kink 
field is right (-I-) or left (— ) moving. 

To compute this correlator, we again expand to the two lowest contribu- 
tions: 

/oo J/Q 
-^{^P%{x,T)AU9l)){A^{9,)^|;^{0)) 

1 /•oo JO jg 

The first contribution, the one particle contribution, comes from the kink 
excitation, A^, destroyed by the field, V'a- The second contribution, a two 
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particle contribution, arises from kinks, Ap, of opposite chirality to A^, 
and Major ana fermions, Aa- (This reflects the group theoretical fact that 
the tensor product of an 5*0(8) spinor representation with an SO (8) vector 
representation gives the other SO{8) spinor representation [277].) The first 
contribution not included is a bound state-kink pair. It begins to contribute 
at w = (1 + V3) m. 

From the form factor expressions from Section 6.5, we can then write 
down the expression for the spectral functions, ^p//j(a;, k), 

Ap{uj,PkFi + k,kyi) 



Ah{LU,-PkFi + k,kyi) 
|2 uj + Pk 



TT Cp 



m 



S{u) — s/k"^ + vn?) 



1 



X exp 



Pk (u;2 _ ^2 _ ^2)2 ^ 

dx Gf{x) ( _^ , \ ^ ( xd\2 



X 



- 



TT 



/j;2 _ 4^2 

, (6.58) 



where Ap the (unspecified) normalization of the two particle kink form- 
factor, Gf{x) is given in Eq.(6.43), and c± is found in Eq.(6.45). For P = 
R = + (i.e. right- moving electrons/kinks), this function is plotted in Fig. 33. 






k/m 



Figure 33. Plot of the single particle spectral function for right moving kinks. The more darkly 
shaded region corresponds to greater spectral weight. 



The parabolic line in Fig. 33 arises from the single particle form factor 
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contribution, and represents the standard dispersion relation of a particle of 
mass, m. Above this curve comes the two particle form factor contribution 
to the spectral function. This contribution is bounded by the curve, uj = 
Vk"^ + 'im?, and so the single particle states do not cross into the two particle 
region. As can be seen from Eq. (6.58), the two particle contribution opens 
up at threshold with a square-root singularity, indicative of the van-Hove 
singularity in the density of states. 

The plot is manifestly chiral with weighting greater for A; > than for 
A; < 0. This is to be expected as we are plotting the excitations linearized 
about the Fermi momentum, +kFi- The heavier weighting for > indicates 
that is easier to create excitations above the Fermi sea than below it. It is 
interesting indeed that excitations below the Fermi surface can be created 
at all and is a mark that interactions are at play. 

6.6.3. STM Tunneling Current 

In this section we study the tunneling between a metallic lead and the carbon 
nanotube/Hubbard ladder through a point contact. Our starting point is 
a Lagrangian describing the nanotube/ladder, the metallic lead, and the 
tunneling interaction: 

^ = ^SO{8) + ^lead + ^tun ■ (6.59) 

Cso{8) is the Lagrangian of the 5*0(8) Gross-Neveu model. 

The electron gas in the lead is, in general, three dimensional. However, 
in the context of tunneling through a point contact, the electron gas can be 
mapped onto an one dimensional chiral fermion (see for example [8,12,57]). 
The general idea is well illustrated by its application to the Kondo problem. 
There an electron scatters off a spin impurity at x = 0. The scattering 
is determined by the electron operator, '>p{x = 0). As tpiO) only depends 
on its spherically symmetric, L = 0, mode, one can consider the scattering 
electron in terms of an ingoing and outgoing radial model defined on the 
half-line, r G [0, oo]. Unfolding the system onto the full line leaves one with 
a chiral fermion. We emphasize however that the map requires no special 
symmetry; the result is exact regardless of particular anisotropics [12]. As a 
consequence, we write Ciead ^ 

Aead = ^^'^5.-*, (6.60) 

where ^' is a massless, left moving fermion, and z = (t + ix)/2. 

It remains to specify Ctun- In order to preserve charge, the electrons must 
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couple to the kinks of the SO (8) Gross-Neveu model, the excitations with 
the quantum numbers of the electron. Thus 

+gR[^\T)^l{T) + ^Pl{T)^{T)]6{x) . (6.61) 

Here we have coupled the lead electrons to both the right and left moving 
fields creating the kink, a, and have allowed the two couplings, gi and gn, 
to be unequal. However as we will work to lowest non-vanishing order in 
the tunneling matrix elements gi/B., the tunneling current will depend upon 
the sum, gf^ + gj^, that is, the contribution of the left and right channels to 
tunneling will add linearly. Similarly, permitting other kinks to couple to 
the lead electrons will give lowest order contributions which simply add. 

15 



10 



5 





1 2 3 4 5 

V/m 

Figure 34. Plot of the tunneling current as a function of applied voltage. The dashed curve 
describes tunneling into a non-interacting fermionic system of mass, m. The vertical dashed line 
marks where a second set of two particle states begins to make a contribution. 

The calculation of the current to lowest non- vanishing order in g^^ jji fol- 
lows the standard route. The tunneling current operator is given by 

I{T)=igR{r+{T)^{T)-^!\T)r+{T)) 

+igL(V'5(T)^'(r) - ^'t(r)V^^(r)) . (6.62) 

In order to induce current flow, {I{t)), through the point contact, one biases 
the lead with a voltage, V. This bias can be taken into account via a gauge 
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transformation, 

*(t) ^ e'^'^^ir) . (6.63) 

In effect we have shifted the energy levels of the electrons. Treating the 
couplings, gh/gRi with linear response theory, we find 

-e^^^{V'5(r)V'^(0))(*(T)^'t(0))l|__,^+, 

V^-iV 

+ (L-^i?,C/r->^?/V'?). (6.64) 
The lead electron correlator (^'^^(r)^'(O)) is well known: 

(*t(T,x)*(0)) = {^(t,xWUO)) = . (6.65) 

T + ix 

With this it is straightforward to express the dc current, (/(w = 0)), in terms 
of the single particle kink spectral function, 

{I{lo = 0)) = —J^ckoj dk[glA-iio,k)+glA+{io,k)], (6.66) 

where A±{u},k) = Ap{ui,±kFi + k,kyi) + Ah{-ui,^kFi - k,kyi), and Ap/^ 
are the spectral functions given in Eq. (6.48). We note that as a technical 
point, in deriving the above equation we have displaced, ^, the lead electron 
operator, slightly from x = 0. In this way we cure the UV divergence 
attendant as r — > 0. At the end of the calculation we then take x to 0. 

In the previous section we have computed A±{u!,k) exactly for energies 
uj < {V3 + l)m. Inserting Eq. (6.58) into Eq. (6.66), we find (7(0)) takes the 
form 

m) = ^-^ial + glW - m'y/'e{\V\ - m)sgn(F) 

+0(|y| — 2m)sgn(y) X two particle contribution. (6.67) 

We sec that for \V\ < 2m, the system behaves as a gapped free fermion. The 
first sign that there is any interaction comes for |V^| > 2m where the voltage 
begins to probe the two particle states, a signature of interacting fermions. 

We explicitly plot (/(O)) in Fig. 34. The square root behavior near V/m = 
1 and subsequent linear form is typical of a gapped fermion. At V/m = 2, 
the two particle states begin to contribute leading to a small change in the 
slope of the I — V curve. At V/m = \/3 + 1 (marked by the the vertical 
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dashed line), a second set of two particle states (a mass VSm bound state 
together with a kink) begin to contribute and at this point the result ceases 
to be exact. However as with the current correlators, we expect this higher 
energy contribution to be small. 

The change in slope in the I — V curve at V/m = 2 can be explicitly 
computed. To do so we consider dv{I{0))- This quantity is given by 

dv{im = -ig% + gl) / dk A±iV, k) . (6.68) 

We can thus see 5y(/(0)) directly measures the local density of states at 
X = of the nanotube/laddcr system. 

We plot dv{I{^)) in Fig. 35. The square root singularity at V = m 
signals the singularity of the density of states in an one dimensional system. 
AtV = 2m we see a sudden jump, indicative of the onset of the two particle 
contribution. The height of the jump can be determined exactly: 



dv{IiO))iV/m = 2+) - dv{I{0))(V/m = 2') 

(9% + d) exp [T^^i^- c{x))] . (6.69) 



16^1 
9m 



The region m < V < 2m of 5y(/(0)) completely determines m (by the 
location of the jump), as well as an overall scale (the product of {gi+g%) and 
the constant Ap, normalizing the spectral function). Dividing out these non- 
universal quantities leaves a universal number, characterizing the magnitude 
of the jump: 

16 



exp 



/•~ dx Gf{x) , , , „ , 

/ tV • 6.70 

7o X s{x) J 



This number represents a definite prediction based upon the integrability of 
the model. 



6.7. Effect of Integrable Breaking Perturbations 

It can now be asked how perturbations to SO (8) Gross-Neveu will affect the 
various computations discussed in the previous section. We consider this in 
the broadest terms by focusing upon how the spectrum of SO (8) Gross- 
Neveu is changed under a perturbing term. We do so through straightfor- 
ward stationary perturbation theory, in the same spirit that Ref. [71] treated 
the off-critical Ising model in a magnetic field. The most general possible 
perturbation takes the form 



-ffpert — ^GabGcd , 



(6.71) 
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Figure 35. Plot of the diflferential conductance as a function of applied voltage. The dashed curve 
marks the single particle contribution to this quantity while the solid curves give both the single 
and two particle contribution. The latter plot is exact up to V/m = \/3 + 1 (indicated by the 
dashed vertical line) where a bound state-kink pair begins to make a contribution. 

where Gab, Gcd are 50(8) currents (of unspecified chirality). For such a per- 
turbation it is necessary to consider degenerate perturbation theory. Thus 
in a given particle multiplet (for example, the fundamental fermions in the 
vector representation), the perturbed energies arise through diagonalizing 
the matrix 

Mij = {Me)H^,nA]{e)), (6.72) 

where here the index i,j indicates the particles Ai,Aj belong to the multiplet 
of concern. In the case that Gab = Gcd, Mij is necessarily diagonal, i.e. 
nondegenerate perturbation theory is sufficient. 

It is important to emphasize that this procedure can be handled in the 
context of integrability. The expression in Eq. (6.72) is no more than a form- 
factor which can readily be computed. Moreover as the theory is massive, 
perturbation theory is well controlled. We expect the unperturbed theory 
to describe all qualitative features of the model while the perturbations to 
only introduce small quantitative changes. 

We first consider the consequences of such a perturbation on the optical 
conductivity. They are two-fold. We expect the exciton peak (found, for ex- 
ample, in Fig. 32 at A; = and uj = \/3) to split. In the unperturbed model 
the peak results from a single rank-two bosonic bound state coupling to the 
current operator. When the matrix Mab above is diagonalized, this partic- 
ular state should be mixed into many others resulting in several states that 
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couple to the current operator. However we do not expect the functional 
forms of the exciton peaks to change: they should remain delta functions. 

They must do so provided the perturbation is not so large as to push the 
exciton peak past the threshold of two particle states where it then conceiv- 
ably could decay. As there is a gap between the excitonic peak and the two 
particle threshold, this will not happen for small perturbations. Experimen- 
tally changes to the exciton peak may not be detectable. Given that any 
experiment will be conducted at finite temperature, the excitonic peak will 
be thermally broadened, perhaps washing out any splitting of the original 
zero temperature peak. 

Wc also expect the perturbation to affect the onset of the two particle 
threshold, although in a less dramatic fashion. Like the unperturbed case, 
there will be several two particle contributions to the optical conductivity. 
However unlike the unperturbed case, the thresholds of the two particle con- 
tributions will not all occur at a; = 2m but be distributed about this energy. 
Thus the two particle contribution is arrived at (approximately) by super- 
imposing several slightly shifted two-particle contributions similar to that 
found in Fig. 32. But given the optical conductivity vanishes at threshold, 
the qualitative picture remains effectively unchanged (i.e. the superimposed 
contributions will appear nearly identical to the original picture). That the 
optical conductivity vanishes at the two particle threshold is a result of the 
vanishing of the relevant matrix element at threshold. This should be robust 
under perturbation as it is ultimately a consequence of the mere presence of 
interactions and not some particular type of interactions. 

We can also consider in a similar spirit the approximate effect of perturba- 
tions breaking intcgrability on the tunneling conductance. As the tunneling 
conductance is determined directly from the single particle spectral function, 
we can deduce how the former is affected from how the latter is changed. 
At a given momentum, the single particle contribution to the single parti- 
cle spectral function under perturbation comes at a discrete set of energies. 
In terms of the tunneling conductance, we expect a series of closely spaced 
square root divergences (a sawtooth behavior) about V = m indicative of a 
series of van-Hove singularities. As the perturbation is removed these singu- 
larities would collapse on top of one another leaving the original picture in 
Fig. 35. 

In the unperturbed case the two-particle threshold is characterized by a 
jump in the differential conductance. Under a perturbation, this jump would 
become a staircase or a series of smaller, closely spaced jumps. This is a 
reflection of the series of van-Hove singularities found about u) = Vfe^ + ^rri^ 
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in the two particle contribution to the single particle spectral function. 
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